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New Frontiers 


By CARL 
Chrysler Corporation, 


ROGRESS in any form of human endeavor 
depends largely on finding new things to do 
and the discovery of better ways of doing old 


things. The very word “research” implies not 
only a search for new things and new ways of 
doing them, but also a re-searching of old things 
and finding new uses and applications for them. 
re-examination of the 
old ways of doing the things we now do. 

All progress 


This entails a constant 


individual, corporate or national 

is based upon the understanding of problems 
id, secondly, the solution of those problems. 
\mong savage peoples and, indeed, in all simple 
societies, men have been content because their 
problems were few and elementary—and their 
progress was correspondingly slow. They used 
the ol “trial and error’’ methods 


he old “‘hit or miss,” 
of problem solution. 
It has been the same in the business world. 
\Vhen business was simple there was not so 
much need of research, but today we see a 
society that has become complex—complex in 


> Tequ 


ements, its mode of living, its needs 
ind aspirations. So it is natural, and inescapable, 
that business, too, should have become compli- 
cated for it lives only by satisfying the multi- 
\udinous needs of today’s complex society. 

: lo meet the new demands made upon it, 
business must find new things to do and better 
loing old things. It is senseless to 


“ontinue doing old things in old ways in a 
| from an address by Dr. Carl Breer at the 
ociation of Manufacturers’ Congress of 
ustry, December 8, 1938. 
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for Industry” 


BREER 
Detroit, Michigan 


world that demands new and better things. And 
to find the ways and means of doing these new 
things in the best possible manner is the essential 
function of research. 

I have devoted my life to research. To me it 
is high adventure. All true research men find in 
their work the call, the lure and fascination to 
which, in other days, the pioneer responded. 
They must be dreamers, but their dreams must 
come true. To them the future and far-off 
horizons are still entrancing. Today many miss 
the thrill of the frontiersman because they 
mistake the frontiers. Geographical frontiers 
have almost vanished. You may regret that 
there are few new lands awaiting discovery, few 
shores unexplored, few hidden bays in which to 
drop anchor. But why worry! The early pioneer 
split rails. The modern pioneer splits molecules! 

These new modern frontiers are not reached 
on foot, on horseback or in creaky old ships. 
They are sought by active, imaginative, in- 
quiring minds, but, once found, they yield the 
same satisfying returns to our modern pioneer 
as the discovery of new and fertile lands gave to 
the geographical pioneer of former days. 

No modern industry, or business, can survive 
or function efficiently—and certainly cannot 
hope to progress—without well organized and 
intelligently directed research. Examples by the 
score are all around you of the discovery of just 
one new law of physical science leading to the 
creation of a new industry or the revolutionizing 
of an old industry or product. And thousands 
or, perhaps, millions of others as yet unknown— 












because they are beyond our present horizons 


await our coming. 





A common expression among executives is, 





“If we only had a financial endowment what a 





wonderful job we could do in the discovery and 





development of new things.”’ 





There was a time in our industry when we 





had that very viewpoint, but we soon found it 





was not money we needed primarily. In reality 





it was new ideas and men with ideas. An idea 





does not cost much to explore, and a valuable 





idea is not difficult to sell for development costs. 





When bankers or investors look at industry, 





they are attracted by past and present perform- 
ance, but in addition they want to know, and 






have a right to know, about futures. No pro- 





gressive manufacturing industry can succeed 





unless it -has futures in progress. The best 





investment in the future any enterprise can have 





is a good research organization. 





In our yearly-cycle industry, success certainly 





depends upon futures. In your industry, in any 





industry, success depends upon futures. Why 





not stimulate futures? Open up the fields to 





men of imagination, guide them, become inter- 





ested in what they have to offer, give them 





facilities to proceed and they will leave a wake 





behind them by which all mankind will benefit. 





By coordinated and organized effort new in- 





dustries that we dare not even dream of today 





will be opened up. New industries cannot start 





without creating happiness, through the added 





employment and enjoyment they give to human 





beings. 





With proper emphasis upon research the next 





century will be outstanding in its accomplish- 





ments. We know much about science, but what 





we now know compared to what remains to be 





learned is nothing. Research is not a complex 





thing. It is simply recognizing the fact that the 





universe itself is, so to speak, an engineered 





project and then attempting to understand and 








utilize as many of its fundamental laws as 
possible. 

The laws are the same today as they were 
10,000 years ago. We must simply try to learn 
more ways to use them. We must try to use 
some of the same ingenuity as the man who 
first made the wheel. All he did was to use 
existing materials and put them to use. We can 
still do that. 

If you have a research program in your 
business, cooperate fully with the men who are 
carrying it on. If you do not have such a program, 
may I recommend that you move in that direc- 
tion soon. Organize your research wisely, then 
stand back of it sympathetically, and if you are 
not in too much of a hurry for results, it will 
repay you handsomely in satisfaction, and in 
dollars and cents. 

But lest I be misunderstood, let me add that 
there are two important, two essential factors 
necessary. These are interest and guidance on 
the part of management, and a_ reasonable 
investment. It will take some money, of course, 
yet more adventures in research have come to 
grief because of lack of interest on the part of 
management than because money was lacking 

But more important still it will take interest 
on the part of management to secure the right 
men. The genuine type of research mind is not 
common, but such minds do exist and can be 
found. It is management's responsibility. 

In conclusion, research in industry is simply 
search, and then search some more, and, having 
done that, start searching all over again. In a 
changing world, a business cannot stand still. 
It either moves forward or backward. If we 
move by accident, by outside pressure, or b) 
guess, our business is probably doomed to 
failure. But if by intelligent research we direct 
our products and our policies toward new 
horizons to meet the new services needed, our 
progress will be sure and sound. 
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A Review of the Photoelastic Method 
of Stress Analysis. II 


By RAYMOND D. MINDLIN 


Department of Civil Engineering, 
Columbia University, New York, New York 


Part IV. Independent Determination of the 
Principal Stresses 


20. METHODS FOR MEASURING 
PRINCIPAL STRESSES 


HI: techniques described in Section 15 make 
possible the determination of the directions 
of and the difference between the two principal 
stresses in the plane of the plate. As we have 
seen, at an unloaded boundary of the plate, the 
principal stresses normal to the boundary vanish, 


so. that obtained is 


the information already 
sufficient to determine the principal stresses at 
the boundary independently. In many cases of 
engineering interest, the greatest stresses in the 
model and prototype occur at a free boundary 
ud, since these stresses are the most important 
i wide range of design problems, it is often 
unnecessary to proceed further with the inves- 
ition. However, there are also problems in 
knowledge is required of the complete 

of stress at points in the interior of the 

ile or at points along a loaded boundary. 
lhese requirements have led to the development 
idaptation of the following methods for the 
ndependent the principal 
which will be described in Sections 21 


determination of 


graphical integration method, 
’s graphical method. 
extensometer method. 
interference method. 
mcentration method. 
nterferometer method 
nterferometer method. 

’s converging light method. 
ne method, 

| method, 

| method, 


n assisted by the Ernest Kempton Adams 
sical Research of Columbia University. 

last half of this paper. The first half ap- 
\pril issue. 
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21. FILON’s METHOD 


A consideration of the equilibrium of a small 
curvilinear rectangle bounded by four nearby 
isostatics (Fig. 21.1) leads to the following two 
equations of equilibrium* which govern the 
variations of the principal stresses along the 
isostatics : 

op p-4q 


Os 1 pe2 


} 
\ 
( 


0q p-q| 


OS2 pi j 


in which s; and se are orthogonal curvilinear coor- 
dinates measured along the p and gq isostatics, 
respectively, while p; and ps are the respective 
radii of curvature of these isostatics (Fig. 21.1). 


35% 


2, 
qt $2 ds, pr 5s, 





ssoclinic 


. 21.1, Curvilinear rectangle bounded by 
four near isostatics. 


Maxwell pointed out!®* that these equations, 
combined with a knowledge of the isochromatics 
and isostatics, are sufficient 


to determine the 


* Coker and Filon, reference 5, p. 142. 





values of p and gq throughout the plate. The 
numerators of the right-hand sides of the equa- 
tions are known from the measurement of relative 
phase retardation and the denominators can be 
measured from the isostatics. The equations can 
therefore be integrated graphically, starting from 
a point on the boundary where the stress is 
known. However, the form of Eqs. (21.1) is not 
conducive to accurate graphical manipulation 
on account of the difficulty of measuring p; and 
p2 With any degree of precision. 

Filon®® has made the method practical by 
showing that the equations may be written in the 
form 


? ) 
p=po+ | (p—g) cot ydo| 
J 6 

\, (21.2) 


g=qo- | (p—q) cot y dd | 


; 


in which y is the angle between the isoclinic and 


the isostatic p, as shown in Fig. 21.1; @ is, as 


iv 
> 
Y) 

‘2 











Fic, 21.2. Filon’s integration method along an isostatic p. 


before, the angle between the isostatic p and an 
arbitrary x direction; py and qo are the known 
values of p and g at the starting point of the 
integration. The procedure is as follows: At a 


A 
wt a 


Az 


yy 
yy 














Ls 


Hh % p 


. . . . > 
Fic. 21.3. Graphical integration of j" (p—q) cot yd¢. 
“ do 





- along an isostatic 
(say p, Fig. 21.2), the values of (p—q), y and 
@ are determined by one of the methods de- 
15. (p—q) cot y is then 
plotted against @¢ as in Fig. 21.3. The area under 


series of points A», A, As, 


scribed in Section 


the curve from @» to, say, 2 is 


(p—q) cot ¥ do 


and this, added to the known value pp at Apo, is 
the stress p at A». Since p—q is known at A, ¢ 
is thus also determined. A similar procedure is 
followed for succeeding points and also for the 
isostatic q. 

The inaccuracies of the method arise chietly 
from the difficulty of precisely determining the 
angles @ and y. These depend upon the isoclini 
lines which are usually observed in the polari- 
scope as broad, poorly defined bands. 

Variations of the method suitable for special 
types of stress distributions have been devised by 
Filon,? Féppl,® Baud,?* and Frocht.® 


22. NEUBER’S GRAPHICAL METHOD 


Still another variation of the form of Eqs. 
(21.1) has led Neuber'® to a method for deter 
mining p and g by a graphical construction not 
This method yields 4 
family of curves, called isopachics, which 


involving integration. 
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ent contours of constant p+q. The sig- that AP=AQ. For the isoclinics, the successive 
ce of the term isopachic will appear in the parameters will be taken to differ by Ag. Let a 
section. Having p+q and p—q at any be the distance between two successive iso- 
the values of p and gq are obtained by pachics, 6 that between two successive isochro- 
on and subtraction. matics, and c that between two successive iso- 
vy, 22.1, O1 and O2 represent the directions clinics, so that a, b and c are measured in the 
sostatics p and q, respectively, through a_ directions 4, 6 and 8, respectively. Then Eqs. 
O. These are the same as the directions (22.1) may be written in the approximate form 
Jong which ds; and dse are measured in Fig. 
21.1. The direction of the isopachic line through AP aon AQ Pp Rat = 
0 is given by O3 and its normal by O04. The a —a j an 9+20 i rsaali 
tion of the isochromatic through O is O5 


ind its normal is O06. Finally, the direction O7 AP ae on aad ei _— 
l 


that of the isoclinic through O and O08 is its a , 2 


il. As before, @ is the angle which the p 
tatic (O1) makes with the x axis and ¥ is the All the quantities on the right of Eqs. (22.2) 
ingle between the isoclinic and the isostatic ~. can be measured from the isoclinics, isochromatics 
In addition, € and » are the angles which the and isostatics, and on the left AP is taken equal 
isopachic and the isochromatic, respectively, to AQ. There are thus only two unknown quan- 
make with the isostatic p. tities, € and a, and two equations to determine 
From Figs. 21.1 and 22.1 we find the following 


tities letting p+q=P, and p—q=) 


aren oP oP 
sin ¥ : =-—siné : 
OSs OSs Os , ea 
' 7 ssoclhnic 


dQ dQ 5 ssochromatic 
—= —sin 7 
OS OS¢ 


3 isopachic 


/ fsostatic 





oP oP 
=cos £ : 


OSe OSs O x 





aQ dQ Fic. 22.1. Diagram used in explanation of Neuber's method. 
—=cos 7—. 
OSe OS¢ 
them. The angle & gives the direction of the 
Substituting these in Eqs. (21.1) there result isopachic at any point. Having these directions 
for many points throughout the plate, the iso- 


90) : 
at 09 pachics may be constructed from them in the 


sin 7+2Q0— cos y ; 

as, OSs same manner as the isostatics are constructed 
from the isoclinics. 

aQ 42 Si ; At a free boundary, P=Q, if the boundary is a 

ase Con gH Ose =e line 1, and P= —Q if the boundary is a line 2. 

Hence the order of an isopachic is determined by 

Neuber’s forms of the equilibrium the order of the isochromatic which it meets at a 

free boundary. Neuber gives an alternative 

stress increments of the successive procedure when the isopachic does not intersect 

the isochromatics be AQ so that the aboundary. 
s of the isochromatics are Q=0, AQ, As in Filon’s method, the accuracy of the con- 
The same stress increments will be struction depends upon the precision with which 





the successive orders of isopachics, so the angles ¢, y and 7 can be measured. 
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23. LATERAL EXTENSOMETER METHOD 


6 


in the 
direction of the normal to the plate is propor- 


Mesnager'’® observed that the strain 
tional to the sum of the principal stresses so that 
by measuring the change in thickness of the 
plate at any point, we may determine p+q at 
that point. Since p+q is constant along contours 
of constant thickness, the origin of Filon’s term 
isopachic, meaning same thickness, is apparent. 

If the plane of the plate is the x,y plane, o, 
vanishes and the strain in the direction of the 
plate normal becomes, from Kqs. (7.1), (10.7) 
and (10.8), 


€-=—(v E)(p+q). (23.1) 


Taking d as the thickness of the plate and Ad as 
the change in thickness we have, therefore, 


p+q= —AdE/ vd. 


(23.2) 


To obtain an idea of 


the magnitude of Ad, 
suppose the model is of Bakelite, for which we 
may take E= 600,000 Ib. per sq. in. and v=0.24. 
Further, let p+q=1000 Ib. per sq. in. and 
d=0.25 in. Then d=0.0001 in. Hence we require 
an extensometer which will measure at least to 
a hundred-thousandth of an inch. 

\Mesnager designed an interferometer caliper 
for the measurement of Ad. It consists of a pair 
of bars, AB and CD, Fig. 23.1, hinged together 





A 


— 





E 


Fic. 23.1. Mesnager’s interferometer caliper. 





























L_| 


Cc 














at the ends A and C with a flexible steel plate. 
A pair of contact points at the intermediate 
points E and F engage the model so that as the 
thickness of the model changes there is a some- 
what magnified change in distance between the 
ends B and D of the bars. The relative motion is 
measured by observing the movement of inter- 
ference fringes formed in the air film between a 
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pair of optical flats fixed at B and D. If the 
5461A line of the mercury spectrum is used as a 
light source, a movement of one fringe corre- 
sponds to a displacement of 0.00001075 in. 











An instrument based on the same principle as 
Mesnager’s but of improved design has been used 
by Vose.'” 

The Westinghouse Research Laboratory has 
adapted the Huggenberger tensometer to the 
measurement of Ad as shown in Fig. 23.2. 

In Coker’s® lateral extensometer, Fig. 23.3, the 
gage points actuate a mechanical lever which, in 
turn, rotates a small concave mirror. The sweep 
of a beam of light reflected from the mirror is a 
measure of Ad. The instrument is mounted on a 
movable platform so that the changes in thick- 
ness all along a line in the plane of the model 
may be determined by traversing the line before 
and after loading. The three instruments previ- 
ously mentioned have an advantage over Coker’s 
extensometer in that their gage points move with 
the model as the latter deflects under load. 

The Ad method for determining p+ is one of 
the best available in that (1) the measurements 
are independent of p—gq and ¢, (2) the measure- 
ment at one point of the model is made inde- 
pendently of the measurements at other points, 
(3) the apparatus is inexpensive. 


24. Arr-F1LM INTERFERENCE METHOD 


A procedure applied by Maris' involves the 
determination of Ad by observing interference 
fringes in the thin air film formed between the 
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of the model itself and an auxiliary 
flat. The method has been further de- 
d by Tesaf®? and Frocht®* 7° so that the 


field of the isopachic contours may be 


viewed or photographed. 


Since the interference fringes are loci of points 
nstant air film thickness, they will also be 
f points of constant model thickness if the 


iry optical flat is properly oriented with 














23.3. Coker's lateral extensometer. 


to the model. Successive fringe orders 


hen represent isopachic contours of equal 


ments of p+q. 


rding the proper 


orientation 


of the 


flat with respect to the model, Tesa#!®? 


ribed the following procedure. Three 


ot on the same straight line, are chosen 


boundaries of the 


() 


The 


model. 


values 


of p—q at these points are known 


isochromatics so that the values (P;, Pe, 
+g are also known at the same points. 
iary flat is then adjusted so that the 
he isopachics at the three points are in 


P, m P» 


: Ps=Q1 : Qe : Qs, proper con- 


being given to algebraic sign as ex- 
Section 15. This procedure can be used 
surface of the model is so well polished 
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that, when the model is not loaded, the surface 
is flat to a degree commensurate with the least 
error desired in the Ad measurement. For models 
not so well polished, a preliminary set of fringes 
may be the unloaded model to 
determine surface irregularities so that appro- 
priate corrections may be applied to the isopachic 
fringes. Tesa¥ also suggests making observations 
on both faces of the model. For each face the 
fringes determine the ordinates to the model 
surface with respect to a plane determined by 
three points on the edge. The relative orientation 
of the two reference planes is then determined 
from the known thickness changes at the three 
pairs of points. 


observed on 


The method involves an elaborate and delicate 
technique but it has the advantage that it gives 
a complete set of isopachic curves rather than 
point by point information. 


25. STRESS CONCENTRATION METHOD 


This method, attributed to Tesa¥!® and Baud! 
makes use of the Kirsch solution* for the stress 
concentration at a small hole in a uniform ten- 
sion field. Referring to Fig. 25.1(a), if p is the 
uniform tension in a plate far from the hole, the 
stresses at A and A’ are 3p and at B and B’ are 
—p. If on this stress field there is superposed a 
stress g at right angles (Fig. 25.1(b)) the stresses 
at A and A’ will be Py =3p—q and P,, = — p+3q. 
Py and P,, will be, respectively, the maximum 
and minimum stresses around the hole. Then 


p=5(Pn+3Pm) 
g=4(Put+3Pn) , 
p 


tt tt 


“Pp 
8 


p 
ttt t 


“p+3q 


SpiA A'\3p 


tht 


(a) 


pai. 
Hid | 


(b) 


Fic. 25.1. Diagram used in explanation of stress 
concentration method. 


* Coker and Filon, reference p. 5, p. 482. 





To determine p and g at a point in a photo- 
elastic plate it is proposed to drill a small hole at 
that point, determine Py and P,, from the result- 
ing isochromatics on the boundary of the hole 
(Fig. 16.1) and calculate p and g from Eqs. 25.1. 

There are at least two difficulties involved in 
the method: (1) It does not seem likely that a 
small hole can be drilled without producing 
“edge effect ;’’ (2) in the light of Vose’s experi- 
ments (see Section 13) the values of P,, and Ps 
might not be reliable since the diameter of the 


hole 


would have to be of the same order of 


magnitude as the thickness of the plate. 


26. FAVRE’S INTERFEROMETER METHOD 


As we have seen in Section 11, measurement 
of the absolute phase retardations of the waves 
polarized in the directions of p and gq will give p 
and g directly. Favre: **:°* has designed an 
the Mach-Zehnder 
effect the required measurements. In this instru- 


interferometer of type to 
ment (Fig. 26.1) light from a source S is polarized 
by a prism P and is split into two rays, 7 and R, 
by the half-silvered 1/,. The transmitted ray (7) 
is reflected at the full-silvered mirror M2, passes 
through the model .M, and is reflected into a 
telescope by the half-silvered mirror, M3. The 
other ray (R) is reflected at 4, passes through 
M;, and rejoins the first ray. A half-wave plate 
II, is used to rotate the plane of polarization so 
the the 

a point in the 
determined previously). 


as to coincide with one or other of 
principal stress directions at 
can be 


model (which 


Msg c 


~- 
-. 


/ 


M2 


Fic. 26.1. Favre's interferometer-polariscope. 


The second half-wave plate (H2) rotates the plane 
of polarization back to its original position. A 
glass plate C is rotated to bring the phases of the 
rays 7 and R, traversing the separate paths, into 
coincidence when the model is not loaded. With 
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retardation, 


the plane of polarization in the direction of, say. 
p, the model is loaded and the phase of the ray 
T, is altered with respect to that of the ray R 
for two reasons. First, the index of refraction 
for a wave polarized in the direction of p changes 
from mp to nm, and, second, the thickness of the 
model changes by an amount Ad. The phase 
difference between 7, and R is then 


A, =(22/d)[d(mi— my) +(mo—1)Ad], (26.1 


the notation being the same as in the preceding 
sections. A measure of A, is obtained by recording 
the rotation of the glass plate C required to 


bring 7°, and R back into phase. In the same 








Fic. 26.2. Favre apparatus. (J. H. A. Brahtz, U.S. Bureau 
of Reclamation). 


manner, the plane of polarization of T is rotated 
so as to coincide with g and we find another phase 
this T, and R ol 


time between 


magnitude 


As=(2r/d)[d(n2—mo)+(mo—1)Ad]. (26.2 


Substituting into Eqs. (26.1) and (26.2) the 
values of (m;—m9), (m2—mo) and Ad given in Eqs. 
(11.1) and (23.2), we find 


A, =adp+bhdg | 
A» = bdp+adq }" 


2a , . 
. Ci- (mj—1) 
4 E , 


, Bs 
lc, (no—1) ] | 
: E 


where 








/ 


rr ] 
The constants a and 6 depend upon the model 
material and the wave-length of the light source, 
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iy be determined in terms of the angle of 
of the glass plate C by a simple tension 
exp lent. 

Hence, knowing 4), As, d, a, b in Eqs. (26.3), 
we may solve for p and g. This procedure is 
ried out point by point over the entire plate. 
Very accurate measurements can be obtained by 
this method if (1) the model surfaces are accu- 
rately plane and parallel (glass is usually used), 
?) the interferometer-polariscope and the load- 
device are well designed and constructed, 
proper precautions being taken especially to 
maintain constant temperature. A photograph 
of the Favre type instrument built at the photo- 
the U. S. 
Reclamation is shown in Fig. 26.2. 
lank 
Favre's interferometer may be used to measure 

g) and (p+q) point by point. 


laboratory of Bureau of 


elastu 


has suggested a procedure whereby 


27. FABRY'’S INTERFEROMETER METHOD 


In this method,®” the two faces of the model are 


half-silvered and are used as _ interferometer 


surfaces. Parallel light is passed through the 
it normal incidence and the interference 


model 


hich takes place between rays transmitted 
directly through the plate and rays undergoing 
27.1) at the 
served through an analyzer. The light entering 


he loaded model is split into two waves, p and q, 


iteemen 
+ ——— senor 
95 ane 
A 


% 
2 


two reflections (Fig. surfaces is 




















hic. 27.1. Fabry’s method. 


wized at right angles to each other. 
hese waves is partly transmitted (p, 
nd partly reflected (p2 and gz) at the 
face. The reflected waves (2 and qo) 
nother reflection (at the first surface) 
pass out of the model in their original 
nd through the analyzer with p, and 
\ing the analyzer, one pair, say g; and 
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F1G. 28.1 Index ellipsoid showing circular sections. 


gz, may be extinguished, since they are both 
polarized at right angles to p; and pe, and inter- 
ference will be observed between p; and p2 due 
to their path difference through the plate. By 
rotating the analyzer through 90°, interference 
will be observed between gq; and qe. 

Again, this is necessarily a point by point 
method since the polarization directions of p 
and g vary over the model, and also the médel 
surfaces must be accurately polished plane and 
parallel. 


28. HILTSCHER’S CONVERGING LIGHT METHOD 


Returning to our consideration of the index 
ellipsoid (Section 8 and Fig. 8.3), we observe that 
there are two planes BOD and BOD’ in Fig. 28.1) 
containing the intermediate principal axis, OB, 
which intersect the ellipsoid in circles. The 
normals, OQ and OQ’, to these planes are called 
the optic axes and the included angle, 2Q, is the 
opticaxial angle. When the differences between the 
principal refractive indices m., m», m,. are small, 
as they are in photoelasticity, it may be shown, 
from the geometry of the ellipsoid that 

Na— Ny 
sin? Q=——— 
Na—N- 


‘ Nb — Ne 
cos* (.}=————_-. 
Na— Ne 


(28.1) 


’ 
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We also recall (Section 8) that any radius ON 
(Fig. 8.3) represents a possible direction for the 
that 
perpendicular to ON intersects the ellipsoid in 


wave normal and the plane through O 
an ellipse DE whose semi-axes OD and OE are 
proportional to the indices of refraction nm, and 
no of the two waves traveling through the doubly 
refracting medium. Now, if 6; (Fig. 8.3) is the 


angle between the wave normal and the optic 


Sx 














s 4 


p p A’ 





x 


Fic. 29.1. Strain and displacement along the x axis. 
axis OQ, and if 62 is the angle between the wave 
normal and the optic axis OQ’, it may be shown, 
from the geometry of the ellipsoid that 


N;—Ne=(nN,—N,) sin 0; sin Oe (28.2) 


the differences between the principal 


indices are small. 


when 


Now, from Eq. (9.1) we find 
N.—Na= C(o3—0)1) 
Na —N,= C(a,;—<a2) }, 
n, —n.= C(o2—@3) 


so that, combining Eqs. (28.1), | 


and (8.2), there results 


1 

* ain 6, sin Oo 

sin? 2 

sin 6, sin 

cos? Q 

2ndC sin 6; sin 

We have already seen how to 


the quantities on the right of Eqs. (28.4) except- 


ing Q, 6, and 6. Hence, if we can measure these 


three angles, we can determine the three principal 
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measure all of 


stress differences. Further, since at least one of 
the three principal stresses is zero at any point 
in a two-dimensional stress system, the other two 
can then be determined directly from Eqs. (28.4). 
Hiltscher*! has measured &, @; and 62 by means 
of converging polarized light, a technique widely 
used by petrologists* for determining the cor. 
responding information for natural crystals. 
Another technique might 
adapted involves the use of the universal rotating 


which also be 


stage.** 


29. MEMBRANE METHOD 


Consider a small strain ¢, in a line segment PA 
(Fig. 29.1) of length 6x parallel to the x axis. The 
displacement of P to P’ is u and the displacement 
of A to A’ is u+(du/dx)bx. Hence the strain of 


PA is 
P'A'—PA 
a | 


€, 


z 


[(éx—u)+(u+(du/dx)ix)|—bx du 
=— =—, (29.1 
6x Ox 


Similarly, the strain ¢, parallel to the y axis is 


y 


4 
T 


P A 


FG. 29.2. Shearing strain in the x,y plane. 








expressed in terms of the displacement 9 in the 
y direction by 


y= ov Oy. (29.2 


To determine the corresponding relations for 


*A. Johannsen, Manual of Petrographic Methods, 
Chapters 24 to 34. 
** A. Johannsen, 


Chapters 35, 36. 


Manual of Petrographic Method 
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aring strain y,, in the x,y plane, consider 
inge in angle between adjacent edges PA 
B of an element in this plane (Fig. 29.2). 
lisplaced to P’, A to A’, and B to B’. If 
perpendicular distance from P’ to PA is 1, 


to a first approximation, the perpendicular 


the plate to the end that the material remains 
geometrically continuous. 

Consider, now, the laws which govern the 
variations of the stress components so that the 
forces on an element in the x,y plane shall be in 
static equilibrium. Such an element, with the 


e from A’ to PA is v+(dv/dx) dx. Hence, 
ne of the angle between P’A’ and PA is 


stresses acting on it, is shown in Fig. 29.3. 
Summing the forces in the x and y directions, we 
7 : find the following equilibrium equations : 
(v+(dv/dx)éx)—v dv 


bx ax Oc; IX + OT zy dy=0 


as \ : (29.5) 
, a 00,/dy+07,,/d0x=0 
since the angle is small, dv/dx represents 
between P’A’ and PA. Similarly, 
is the angle between P’B’ and PB. The 
small shearing strain y,, is the total change in 
vle between the two edges. Hence 


angle Differentiating the first of Eqs. (29.5) with 





Yry = 0v/dx+0u/dy. (29.3) 


It is apparent that, when the strain varies 





throughout the x,y plane, the three strain com- 
ponents €,, €y and y,, are not independent quan- 
tities but are related through the two components 
of displacement u and v. To find the relation, 
differentiate Eq. (29.1) twice with respect to x, 
Iq. (29.2) twice with respect to y, Eq. (29.3) 








ce with respect to x and once with réspect to y, 


ind add the results; obtaining F 1G. 29.4. Models used in membrane analogy 
(E. E. Weibel). 


Pe,/ OV? +07 €,/0x* =A y2,/dxdy. (29.4) 


(his is known as a compatibility equation and respect to x and the second with respect to y, and 
expresses the law which governs the permissible adding, we obtain ' 


ions of the strain components throughout 


Gyr ia by 


o,/dy?+0?0,/dx°+20*7,,/dxdy=0. (29.6) 


Now, the strains are related to the stresses 
through Hooke’s law (Eq. 7.1). In the case of 


Tayo gee Sy plane stress in the x,y plane, we have 





Taye Sy 5x 





Ee,.= 
Ee, =0,—vez 
Eyry=2(14+ 7) rey 


7} 
Oz— VO, 


r (29.7) 
) 





Substituting Eqs. (29.7) into (29.4) and using 
(29.6), we find 




















(0?/dx?+0?/dy?)(o,+0,) =0. (29.8) 


Since, from Eqs. (10.7) and (10.8), (¢:+¢,) 


7" =(p+q), we have, finally, Laplace’s equation 


resses acting on a rectangular element in the 


x,y plane. (0?/dx?+0?/dy")(p+q) =0 (29.9) 
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as the equation governing the variation of the 
Den 


ol served 


Hartog,*” 
that 
29.9) is the same as the equation governing 
the 


sum of the principal stresses. 
Biot®! 
lq. 
the 


and Biezeno and Koch*® 


ordinates to surface of a uniformly 
stretched membrane which is deflected a small 
amount by forces applied only at its boundaries. 
To make use of this analogy, a hole of the same 
shape as the boundary of the model is cut in a 


plate, and around the edge of the hole ordinates 








Fic. 29.5. Membrane analogy apparatus (University of 


Michigan ° 


are erected normal to the plate and proportional 
to p+q 
membrane is stretched across the opening and is 


at each point of the boundary. The 


constrained to follow the boundary ordinates. 
The ordinate to any point of the membrane is 
then proportional to the magnitude of p+q at 
that point. 

The method is very easy to apply and is 
applicable to all regions in the neighborhood of 
boundaries along which p+q can be previously 
determined. It has been applied with success by 


Weibel!77: '7> (see 


film 


Figs. 29.4 and 29.5) using a 
Biot 
and MeGivern and Supper"! ' with a 


soap as the membrane, and by and 
Smits,” 
thin rubber sheet as the membrane. The method 
is necessarily limited to regions in the neighbor- 
hood of boundaries along which the stresses are 
completely known. 


30. ELECTRICAL METHOD 


If the boundary values of p—gq are known, Eq. 
(29.9) may be solved in other manners than the 
membrane method. 
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It is known, for example, that potential of a 
two-dimensional electric field is governed by the 
same equation. If, then, a metal plate is cut to 
the shape of the model and an electrical potential 
proportional to p+q at each point is impressed 
along the boundary, the potential at any point 
in the interior of the plate will be proportional 
to p+q at that point. 

The 
Koch,*® 


Tank.'' 


method is described by Biezeno and 


Biot,*" Malavard,'®? and Meyer and 


31. NUMERICAL METHOD 


Liebmann® has shown that Laplace’s equation 
may be solved for given boundary values by a 
process of successive numerical approximations 
In this method, a network of rectangular lines js 
superposed on the region, and approximate 
values of p+q are assigned to each pe int of inter- 
section. Liebmann pointed out that the value of 
the function at any point should be equal to the 
average of the values at the four neighboring 
points. Hence, by repeatedly traversing the 
network, adjusting the value at each point with 
regard to the values at the four neighboring 
points, a set of values will be reached which fits 
the boundary values. 

Shortley and Weller'® have devised methods 
for greatly accelerating the convergence of the 


pre CESS. 


Part V. Special Developments and Applications 
32. GoopiER’s HypRODYNAMIC ANALOGY 


Goodier™ has described a very complete 


analogy between the quantities observed in 


photoelasticity and those involved in the slow 
motion of a viscous fluid in two dimensions. The 
isopachic lines correspond to ‘lines of constant 
fluid vorticity; the lines orthogonal to the 1so- 
pachics correspond to lines of constant fluid 
pressure; the isochromatic lines correspond to 
lines of constant maximum rate of shearing 
strain (distortion) in the fluid; the isostatic 
network corresponds to lines showing the orien- 
tation of the elements having maximum distor- 


tion at points in the fluid. It seems possible that 


this analogy might form the basis of a method for 
solving fluid flow problems by the photoelastic 


method. 
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Biot’s GRAVITATIONAL AND THERMAL 
ANALOGIES 


‘ has shown how the stresses due to 
tional loading or temperature variation 
inalyzed photoelastically without actually 


ducing the gravitational force or heating 
odel. 


sses due to the weight of a large body may 


be calculated from the stresses produced in a 
small model of the same shape loaded on its 
boundaries by a normal pressure distribution 
varving linearly with the depth. 

lhe stresses, due to a steady state of heat 
flow in a hollow cylinder, may be found from the 
stresses in a plate model by slitting the plate so 
is to connect the inner and outer boundaries 
nd displacing or rotating the two edges of the 

dislocation) by an amount which may be 
culated. Weibel'”® has applied the method to 
irious cross-sectional shapes (see Fig. 33.1). 











plication of Biot’s thermal stress 


lalogy (E. E. Weibel). 
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Stresses in a solid cylinder due to a transient 
thermal state may be calculated from a plane 
model of the same shape as the cylinder cross 
section and loaded on the boundary. Biot and 
Smits® have applied this to the investigation of 
thermal and shrinkage stresses in a dam. 

Stresses in a hollow cylinder due to a transient 
thermal state require a dislocation in addition to 
boundary loading. 


34. THREE-DIMENSIONAL PHOTOELASTICITY 


In the preceding sections, photoelastic methods 
for analyzing two-dimensional stress systems 
have been described. The various techniques are 
not directly applicable to the analysis of three- 
dimensional distributions of stress because, in 
these, both the magnitudes and directions of 
the three principal stresses vary along the line 
of propagation of the light. Only integrated 
effects over the total path may be observed, 
whereas the optical properties of each element 
along the path are required. No analytical 
procedure for relating the integrated effects to 
the state of stress at an interior point has yet 
been evolved. 

Favre® suggested observing the stress-optical 
effect on a small cube of photoelastic material 
model made of a 


embedded in a transparent 


material whose stress-optical coefficients are 
zero; but the experimental difficulties involved 
in the procedure seem insurmountable. 
Mesnager'®* proposed a method applicable to a 
limited class of three-dimensional problems. The 
successful application of the method would lead 


‘ 


to a determination of the ‘‘extreme fiber’’ stress 
produced in a plate by forces applied normal to 
the plane of the plate. The procedure begins 
with the coating of the polished surface of a 
metal plate with a thin layer of transparent 
material which becomes doubly refracting when 
strained. Bending of the plate out of its plane 
strains the coating so as to produce an essentially 
two-dimensional stress system, and the resulting 
double refraction is measured by analyzing light 
reflected from the polished metal surface. No 
successful application of the suggestion seems to 


have been made. 
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Timby and Hedrick'™* have proposed another 
technique for analyzing extreme fiber stresses. 
They suggest cementing to the surface of a 
of Polaroid. 
Over this is cemented a sheet of photoelastic 


transparent model a thin sheet 
material. Light emerging from the model passes 
through the Polaroid so that its ellipticity is 
annulled. Hence only plane polarized light passes 
through the outer layer of doubly refracting 
material. Observations on the emerging light then 
yield information on the two-dimensional state 
of stress in the outer layer. 

In 1851 Maxwell'** performed an experiment 
which has been developed, in the last few years, 
into a very promising three-dimensional method. 
Maxwell took a hollow cylinder of hot isinglass 
and rotated the inner surface of the cylinder, 
with respect to the outer surface, through a 
small angle about the axis of the cylinder. He 
allowed the isinglass to cool with the twisting 
still that the 


resulting isinglass 


observed 
the 
represented an elastic stress distribution. He was 


couple applied and 


double refraction in 
unable to explain the phenomenon and _ the 
matter was apparently ignored until 1935 when 
Solakian' heated a solid cylindrical bar of Mar- 
blette, applied a twisting couple about the axis of 
the bar, and allowed the material to cool with the 
couple still applied. He then cut out a slice per- 
pendicular to the axis of the bar and examined 
the plate in a polariscope. However, the resulting 
fringe pattern did not correspond with the state 
of stress which is predicted by the St. Venant 
torsion theory, as was later pointed out by 
Hetényi.” 

Oppel'’® made a similar experiment with a 
block of Trolon. He heated the block, pressed a 
metal sphere against it, and cooled the block 
under load. He observed that slices cut from the 
cooled block exhibited a fringe pattern corre- 
sponding to an elastic state of stress. 

The recent experiments of Hetényi (see Sec- 
tion 12), and a discussion by Kuske™ furnish an 
explanation of the phenomenon. At the elevated 
temperature, the elastic, infusible skeleton alone 
resists the applied forces. When the material is 
cooled under load, the fusible part freezes around 
the deformed skeleton and essentially maintains 
the the load is removed. 


deformation when 
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The resulting material is thus in a state of 
deformation corresponding to the elastic state of 
stress at the elevated temperature. This state 
of stress is not seriously disturbed by carefy| 
slicing of the model because the equilibriym 
between the fusible and infusible parts exists 
over regions of molecular proportions. A difficulty 
in the procedure is that, while Young’s modulys 
is diminished in 640 :1 at 


the relative 


the ratio of about 
the 
optical coefficient is increased only in the ratio 
of 26:1. the must be 
applied to obtain high double refraction is so 


elevated temperature, stress 


Hence, strain which 
large as to change seriously the shape of the 
model. 

The interpretation of the double refraction in 
the slice is rather more complicated than for 
two-dimensional stress systems. In general, a 
principal plane of optical symmetry does not 
coincide with the plane of the slice, so that the 
relative retardation is not a measure of a differ- 
ence between principal indices of refraction, and 
is accordingly not proportional to a differenc 
between principal stresses. The three principal 
stress differences may, however, be determined 
the additional 
measurement, at each point of the slice, of the 


as Hiltscher has shown," by 
optic axial angle 22 and the angles 6, and 4; 
between the normal to the slice and the opti 
axes (see Sections 8 and 28). 


35. SPECIAL CASES OF THREE-DIMENSIONAL 
SYSTEMS 


There are certain special cases''* for which it 
is not necessary to measure all of the three 
angles Q, 0;, and 42. 

Most important is the case where the stress 
distribution has a plane of symmetry. This will 
occur when both the model and the loading 
system are geometrically symmetrical about 4 
plane. The plane of symmetry is a principal plane 
of stress so that a slice which contains this 
plane will have a relative retardation propor 
tional to one of the three principal stress differ- 
ences. The plane which contains the optic axes 
will be either parallel or perpendicular to the 
plane of symmetry. In the latter case, the 


bisector of the optic axial angle is either paralle! 
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endicular to the plane of the slice. Hence 
ve 6,=0c0=7 2 or 


or pr 
6,;=0.=7/2-—2 or 

Q depending on which principal plane of 
optical symmetry coincides with the plane of the 
slice. Hence, if we measure 4; and Q in addition 
to J, we may calculate the three principal stress 
differences from Eqs. (28.4). 

In a prismatic or cylindrical bar subjected to 
St. Venant torsion,* one of the three principal 
stresses at any point in the bar is zero and the 
oth I 

sien. The zero principal stress direction is 


‘wo are equal in magnitude and opposite 


perpendicular to a generator of the cylinder and 
the other two principal stress directions make 

eles of 45° with the generator. This state of 
stress corresponds to a biaxial crystal with an 


optic axial angle of 7/2, 


the optic axes being 
parallel and perpendicular to the generator. 
\ slice of the bar normal to a generator there- 
fore gives zero phase difference for light at 
normal incidence; but a slice at any other angle 
to a generator will give a phase difference from 
which we may calculate the principal stresses. 
lo effect this calculation, it is necessary to 
know the angles 6; and 6. which the wave normal 
makes with the optic axes. One of these will be 
mstant and equal to the angle between the 
iis of the bar and the normal to the slice, but 
the other will vary over the section. For each 
point in the slice we need therefore measure only 
‘latter angle (#2) and the phase difference A. 
I.qs. (28.4) give the principal stresses di- 

tly, provided the wave-length \, the slice 
kness d, and the stress optical coefficient C, 

re known. It was observed by Hetényi”® that a 
exhibits the maximum number of 
but it should be noted that the fringes 
epresent loci of constant principal stress 
e. It is necessary to know the angle 62 
point in the 45° slice. This angle may be 


SLICE 


ed in the following manner: Let O be a 
e 45° slice, at which 42 is to be found. 
nal to the slice at O and the line through 
to a generator of the bar determine a 
ch intersects the plane of the slice in 
Let OB be a polarizing axis of the 
determined from the isoclinics in the 


henko, Theory of Elasticity (McGraw-Hill 
ny, New York, 1934), p. 228. 
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usual manner) and let angle AOB=a. Then it 
may be shown!'‘ that 


sin 42= (1+ cos? 2a)!. 


Similar formulas hold for slices cut at other 


angles or viewed obliquely. 


36. DETERMINATION OF o3, 2 AND o; IN THREE- 
DIMENSIONAL SYSTEMS 


As Hiltscher®*' pointed out, measurement of 
A, 2, 6; and 62 in any number of slices will yield 
only the directions of the three principal stresses 
and their differences except on a free boundary. 


To determine the principal stresses independ- 
ently in the general case, it is necessary to 
measure ,, m and , or to resort to other than 
optical measurements. 

There are three standard methods for meas- 
uring the principal refractive indices, but none of 
them appears to be suitable for the current de- 
velopment in three-dimensional photoelasticity. 

In Stokes’ method* the slice would be placed 
over a rectangular grid of ruled lines and the 
latter observed through a microscope. On account 
of the double refraction in the plate, each family 
of parallel lines will appear in focus at two 
positions of the microscope objective. The three 
indices can be calculated from the relative focus 
positions. However, the photoelastic double re- 
fraction is so weak that this method does not 
appear feasible. 

The prism method** would require the cutting 
of a small prism from the model for each point 
at which the stress is desired. 

The total reflection methods are very accurate 
but they measure the optical properties of a 
surface layer of the specimen. In a photoelastic 
slice, the surface properties are not a reliable 
indication of the internal optical characteristics. 

Kuske*® has suggested that by heating the 
slice and measuring changes in thickness accom- 
panying the resulting relief of stress, there will 
be sufficient data, combined with that previously 
determined, to calculate o3, o2 and a. 


*J. Walker, Analytical Theory of Light (Cambridge 
University Press, 1904), p. 225. 
** J]. Walker, preceding reference, p. 230. 


+ J. Walker, preceding reference, p. 242. 





37. DyNAmiIc PROBLEMS 


There are three classes of problems involving 
motion to which the photoelastic method has 
been applied : (1) Steady state of stress, (2) cyclic 


variation of stress, (3) transient state of stress. 











Fic. 37.1. Frings photograph of a canti 
lever beam vibrating at 60 cycles per 


second (W. M. Murray). 


Frost and Whitcomb’ have investigated the 
stresses in a disk rotating with constant angular 
speed. This is a case in which the stresses do not 
vary with time. 

Investigations of the cyclical stress in the 
teeth of rotating gears were made by Himes and 
Baud.” Stresses in vibrating beams have been 
studied by Ravilly'’® and Murray.''® Fig. 37.1 is 
a photograph made by Murray of a cantilever 
beam vibrating at sixty cycles per second. 

Transient stresses due to impact have been 
examined by Tuzi,'®* Tuzi and Nisida,!?* "7! 
Frocht,” and Thouvenin.'® Murray''® has given 
a critical review of the difficulties involved in 


transient stress analysis. 


38. MISCELLANEOUS APPLICATIONS 


The photoelastic method has been applied to 
the analysis of a wide variety of problems in 
the fields of machine design, steel and concrete 
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Fic. 38.1. Stress concentration around a 


hole in a tension bar (M. M. Frocht). 


structures, mining engineering, naval architec- 
ture, aeroplane structures, materials testing, 
physical research, glass working, machine shop 
practice, soil mechanics and foundations, railway 
engineering, and many other related branches of 
engineering and applied physics. Some of the 
applications which have been made in the last 
few years are given in the following table. The 


numbers following each item refer to the bibli- 








Fic. 38.2. Stress concentration due to 
semi-circular notches in a tension bar 
(M. M. Frocht). 
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Ribbed slabs"? 

Rotating disks”! 

Rubber manufacture! 1 
Screw threads”: ®* 
Shrinkage stresses*® 

Soil mechanics™ 

Stress concentration®™: ® 157, 
Structural elements**: *; 1° 
Thermal stresses* * 153, 179 
Torsion??: 78 

Tunnel linings™® 

Welded joints*® 48 1 


1G. 38.3. Stress concentration at the fillets 
of a cantilever bar (E. E. Weibel). 


























38.4. Curved beam in pure bending (E. E. Weibel). 


ography in the next section, and the papers 
given there should be consulted for details of 
the different investigations. 


\rch bridges? 107, 127 
bridge trusses” 

r wheels*?. '4! 
crete stresses®*: '5 
da is33) 41, 160 

igue testing™! 1% 

It slabs! 

ud flow7s, 138, 139 


82 
cars - 





ss working*®. 154, 155 
eterminate structures*’ (a) 
s and keyways!? 


‘rials testing: 101, 120, 156 _ 
forced concrete”: 1% Fic. 38.6. Stresses in gear teeth (T. J. Dolan). 
ining walls* 
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F1G. 38.7. Stresses in the main plate of a riveted 


joint (W. M. Murray) 











Space limitation prohibits a detailed descrip- 
tion of so great a number of special studies but a 
few applications for which striking fringe photo- 
graphs have been obtained are noted briefly in 
this section. Figs. 38.1, 38.2 and 38.3 are typical 
of a great many investigations of stress con- 
centrations resulting from holes, notches and 
fillets. Fig. 38.4 shows a curved beam in pure 
bending. The fringe photograph for a loaded 
chain link is shown in Fig. 38.5. Photographs 
from a study of the effect of tooth contour on 
the stresses in gear teeth are shown in Figs. 
38.6(a) and (b). 


The stresses in the main plate of a riveted 


joint are illustrated by the fringe photograph in 


Fig. 38.7. Figs. 38.8(a), (b) and (c) were made in 
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Fic. 38.8. Study of stresses in a bridge rocker (J. H. A. Brahtz): (a) Rocker centrally 
loaded. (b) Enlarged view of region immediately under the concentrated load. (c) Rocker 


loaded near its edge. 


JOURNAL OF APPLIED PuysIcs 





OQ em ne -  -lSO 
EE en 


Des 





i SLO] 
(a) (b) 


Vo] 





Acknowledgments 


I wish to thank Dr. M. Salvadori for his 
kindness in reading the manuscript and for his 
many criticisms and suggestions. 

Mr. D. C. Drucker was of great assistance in 
the preparation of Part III. 

Many of the leading investigators in photo- 
elasticity in this country were asked to con- 
tribute photographs of apparatus and fringe 
patterns for use as illustrations. For their 
courtesy in responding to this invitation I am 
pleased to acknowledge my indebtedness to 
Dr. J. H. A. Brahtz, U. S. Bureau of Reclama- 
tion; Professor T. J. Dolan, University of 
Illinois; Professor M. M. Frocht, Carnegie 
Institute of Technology; Dr. M. Hetényi, West- 
inghouse Research Laboratories; Dr. O. J. 
Horger, Timken Roller Bearing Company; Dr. 
W.M. Murray, Massachusetts Institute of Tech- 
nology; Professor E. E. Sechler, California 
Institute of Technology ; Professor F. G. Switzer, 
Cornell University; Professor E. K. Timby, 
Princeton University; Dr. E. E. Weibel, Uni- 
versity of California. 

Bibliography 

An extensive bibliography, which covers the 

field of Photoelasticity up to about 1930, appears 
38.9, Stresses in a railway car wheel (T. J. in Coker and Filon’s Treatise. The following 
Dolan). list includes papers which have appeared in the 

interval, 1930-1938, and also a few earlier papers 

i study of the stresses in a bridge rocker. Fig. mentioned in the body of this review. The first 
.8(a) shows the rocker centrally loaded and nineteen references are to textbooks and reviews, 


Ky 38 X(b) 




















is an enlarged view of the region while the remainder are recent research papers 
mmediately under the concentrated load. Fig. dealing with theory, technique, apparatus, and 
38.8(c) shows the rocker loaded near its edge. applications. 

A photograph from an investigation of 


tne tr 


esses in a railway car wheel under Motor Control 


| Voltmeter 
ifferent . °,° ° ° Speed Readings 
different loading conditions is shown in : end . 


Fig. 38.9 


An 





Synchronizing 
Switch 
teresting device for the photo- 
ilysis of stresses resulting from 
nal loading has been developed 

Solakian and Baldin.” The 
‘tated in a centrifuge to simu- 
tational acceleration and the 
ern is observed with a strobo- 

iscope (see Fig. 38.10). 





mode] 


. 38.10. Centrifuge-polariscope apparatus (P. B. Bucky, A. G. 
Solakian and L. S. Baldin, reference 42). 


VOLUME ‘0, MAY, 1939 





BIBLIOGRAPHY 


R. V. Baud, “Entwicklung und heutiger Stand der 20. I. Arakawa, “On the determination of stress-optica| 
Photoelastizitat und der Photoplastizitat im Rahman coefficient of Bakelite with initial stress,” Proc 
der Gesamt-Experimentalelastizitat,”’ Schweizer Phys.-Math, Soc. Japan 5, 117-136 (1923), 
Archiv fiir angewandte Wissenschaft und Technik . I. Arakawa, ‘Some advantages obtained frop 
4, 1-15, 48-53 (1938). Bakelite in photo-elastic observation,” Proc, Phys. 
R. V. Baud, ‘Technische Methoden Photoelastischer Math, Soc. Japan 7, 160-180 (1925), 
Forschung,”’ Schweizerische Bauzeitung, 100 (1932). . L. Baes, “La Poutre Vierendeel,” L’Ossature 
». G. Coker, ‘‘Photo-elasticity for Engineers,” Gen. Métallique 5, 457 (1936), and 6, 432 (1937). 
Elec. Rev. 23, 870-877 (1920); 23, 966-973 (1920); 23. L. Baes, “‘L’are élastique parfaitement encastré A ses 
24, 82-88 (1921); 24, 222-226 (1921); 24, 455-466 deux extremités,””’ La Technique des Travaux 13 

1921). 442 (1937). 
>. G. Coker, “‘Photo-elasticity,” J. Frank. Inst. 199, . R.V. Baud, “Further development in photoelasticity.’ 
289-331 (1925). J. Opt. Soc. Am, and Rev. Sci. Inst. 18, 422 437 
G. Coker and L. N. G. Filon, A Treatise on (1929), 
Photoelasticity (Cambridge University Press, Cam 5. R. V. Baud, “Effect of elliptical polarization o; 
bridge, 1931). energy transmission,” J. Opt. Soc, Am, 21, 119-12) 
N. G. Filon, Manual of Photo-Elasticity for Engt- 1931 
neers (Cambridge University Press, Cambridge, 26. R. V. Baud, “On the determination of principal stress 
1936). from crossed Nicol observations,” J. Frank, Inst 
. Féppl, ‘“Spannungsmessung mit Hilfe der optischen 211, 457-474 (1931 
Doppelbrechung,” Zeits. f. tech. Physik 15, 430- . R. V. Baud and W, D, Wright, “Analysis of colours 
436 (1934). observed in photoelastic experiments,” J. Opt, So 
Féppl, “Die Spannungsoptik im Dienste des Am, 20, 381-395 (1930), 
Bauingenieurs,”” Der Bauingenieur 19, 341-5 (1938), 28. R. V. Baud and F, Tank, “Ermittlung von zwei- 
Féppl and H. Neuber, Festigkeitslehre mittels dimensionalen Potentialfeldern bei beliebigen Ran 
Spannungsoptik (R, Oldenbourg, Miinchen and bedingungen,”’ Schweizerische Bauzeitung III, 176 
Berlin, 1935). 177 (1938). 

P. Laurent and A. Popoff, “‘La photoélasticité,” Rev. 2 .. Bergmann, ‘‘Messung elastischer Konstanten mit 
de Metallurgie 35, 363-78, 407-24, 448-74 (1938). Ultraschall,” V. D. I. Zeits. 81, 878-82 (1937). 

G. Mesmer, “Das Spannungsoptische Verfahren,” . 30. A. H. Beyer and A. G,. Solakian, “Photoelast 
die Messtechnik, 11, 217-221, 238-241 (1935). analysis of stresses in composite materials,” Pro 

A. Mesnager, Détermination expérimentale des efforts Am. Soc. Civil Eng. 59, 1121-1132 (1933). 
intérieurs dans les solides, Fascicule X du Mémorial 31. M. A. Biot, “Contribution A la technique phot 
des Sciences Physiques (Gauthier-Villars, Paris, élastique,”” Annales Soc. scientifique Bruxelles BS3 
1929). 13-15 (1933). 

G. Oberti, ‘“‘Fotoelasticita,”” Rendiconti del Seminario 32. M. A. Biot, “Propriété générale des tensions 
di matematico e fisica d. R. Univ. Milano 6, 217-251 thermiques en régime stationnaire dans les corps 

1932). cylindriques. Application A ta mésure photo 

G. Oberti, “Indagini sperimentali sulle costruzioni con élastique de ces tensions,” Annales Soc, scientifique 
l'uso dei modelli,” Atti, Ricerche e Studi, R. Ne en, SSE Ce 
Politecnico di Milano, Hoepli, Milano (1935). 33. M. A. Biot, “Distributed gravity and temperature 

loading in two-dimensional elasticity replaced by 

di sollecitazione in strutture: concetti basiliari e boundary pressures and dislocations,” J. Ap 

: ; it wees : , a Mech. 2, 41-45 (1935). 

recent! applicazioni,” Energia Elettrica 13, 16-23 . = 
1936 34. M. A. Biot, “A general property of two-dimensiona 

thermal stress distribution,” Phil. Mag. 19, 540-54 
1935). 

M. A. Biot and H. Smits, “Etude photo-élastique des 


G. Oberti, ‘‘L’uso dei modelli per lo studio dello stato 


Photoelastic Journal, 55 West 42 Street, New York, 
New York. A. G. Solakian, Editor. 

M. Salvadori, “‘Recenti progressi della tecnica 
fotoelastica,”’ Annali dei Lavori Pubblici 75, 90-110 
1937). 


tensions de contraction dans un barrage,” Bu 
Tech. No. 4, p. 10, Bulletin technique de I’ Union ces 


Ingénieurs sortis des Ecoles spéciales de Louva't 
\V. Tesaf, “‘Les progrés récents et les applications de la 1933) 


photo-élasticimétrie,” La Technique Moderne 30, 


C. B. Biezeno and J. J. Koch, “Uber einige Beispie' 
259-265 (1938). 


zur elektrischen Spannungsbestimmung,” Ingenieu! 
E. Weibel, ‘Developments in Photoelasticity,”’ S. Archiv 4, 384-393 (1933). 
limoshenko, 60th Anniversary Volume (Macmillan 


P. H. Black, “An investigation of relative stresses " 
Company, New York, N. Y., 1938), pp. 257-267. 


solid spur gears by the photoelastic method, 
University of Illinois Engineering Experimet 
Station, Bulletin 288, 1-29 (1936). 


JOURNAL OF APPLIED PHYSICS 





\ 


anjean and F, Temmermann, “Etude de la 
on voisine du point d’inflexion d’une piéce droite 
iatique en flexion plane composée avec effort 
nchant,”” Soc. Roy. Belge des Ingénieurs et des 
ustriels, No. 9, 823-863 (1937). 
Blom, ‘‘Der Spannungs-opiische Koeffizient als 
terialkonstante,”’ Kolloid Zeits. 80, 212-215 


937). 


H 


A. Brahtz, “Photoelastic apparatus at the 
fornia Institute of Technology,” Rev. Sci. Inst. 


8-83 (1934). 


H 


\. Brahtz, ‘Stress function and photoelasticity 


plied to dams,”’ Proc. Am. Soc. Civil Eng. 61, 983- 
100 1935). 


Bucky, A. G. Solakian and L. S. Baldin, 


Centrifugal method of testing models,”’ Civil Eng. 


5, 


) 


287-290 (1935). 


k B. Carleton, “Suitability of materials for photo- 
elastic investigations,” Rev. Sci. Inst. 5, 30-32 
1934). 

E.G. Coker, ‘‘Photoelasticity,”” Engineering (London) 


90, 


ik. G 


london) 129, 465-467 (1930). 


BE. G 


1-4 (1911 
Coker, “Lateral extensometers,”’ Engineering 


Coker and R, Russell, “Stress distributions in 


fusion joints of plates connected at right angles,” 
Institution of Naval Architects 5, 1-8 (1933). 
EE. G. Coker and M, Salvadori, ‘‘Stress waves in the 


vres ol locomotives,” Proc. Inst. Mech, Eng. 131, 
493-512 (1935). 


Coker and B. P. Haigh, “An experimental 


nvestigation of cracking in mild steel plates and 


velded seams,”” Trans. Roy. Inst. Naval Architects, 


77, 193-207 (1935). 


|. W. Cookson and H, Osterberg, ‘‘A formula for the 


reft 


f 


ngence of vibrating media,’’ Physics 7, 166 


1936 


Le 


Ss 


( ’ 


VOLUM] 


t 


Den Hartog, ‘““Experimentelle Lésung des ebenen 
inungsproblems,” Zeits. f. angewandte Math. 
Mech. 11, 156 (1931). 

Edmonds and B. T. McMinn, ‘Celluloid as a 

um for photoelastic investigation,” Trans. Am. 

Mech. Eng. 54, 77-82 (1932). 

ry, “Sur une nouvelle méthode pour l'étude 
nentale des tensions élastiques,” 

is, Paris 190, 457-460 (1930). 
re, “Méthode purement optique de détermi- 


Comptes 


des tensions intérieures se produisant 
es constructions,”’ Schweizerische Bauzeitung 
i 
re, “Sur une nouvelle méthode optique de 
nation des tensions intérieures,” Rev. 
que 8, 193-213, 241-261, 289-307 (1929). 
e, “Sur une méthode optique de détermination 
nsions intérieures dans les solides A trois 
ions,’’ Comptes rendus, Paris 190, 1182-1184 


‘La detérmination optique des tensions 
ires,’’ Rev. d'Optique 11, 1-21 (1932), 


, MAY, 1939 


B. Ferretti, “Optical method for the determination of 
variable elastic stresses,’"” Nuovo Cimento 15, 77-87 
(1938). 

N. G. Filon, “On the graphical determination of 
stress from photo-elastic observations,” Engineering 
(London) 116, 511—512 (1923). 

N. G. Filon, “Exploration of stress near an 
isotropic point in an elastic plate,’’ Phil. Mag. 22, 
187-206 (1936). 

N. G. Filon and H. T. Jessop, ‘‘On the stress- 
optical effect in transparent solids strained beyond 
the elastic limit,’”’ Trans. Roy. Soc. London A223, 
89-125 (1923). 

L. Féppl, “Untersuchung der ebener Spannungs- 
zustande mit Hilfe der Doppelbrechnung,”’ Sitz. d. 
math. naturwiss, Abt. d. Bayerischen Akad. 
Miinchen, 247-265 (1928). 

Féppl, ‘“Neue Erfolge in der Spannungsoptik,” 
Zeits. V. D. I. 81, 137-141 (1937). 

P. Frenkel, ‘‘Photoélasticimétrie et ses applications 
aux constructions civiles,’”’ Travaux Architecture- 
Construction, 19, 33-38, 149-155, 182-187 (1935). 

M. M. Frocht, ‘“‘Recent advances in photoelasticity,” 
Trans. Am. Soc. Mech. Eng. 53, 135-153 (1931). 

M. M. Frocht, ‘‘Kinematography in photoelasticity,” 
Trans. Am. Soc. Mech. Eng. 54, 54-59 (1932). 

M. M. Frocht, “On the application of interference 
fringes to stress analysis,”’ J. Frank. Inst. 216, 73-89 
(1933). 

M. M. Frocht, “Factors of stress concentration 
photoelastically determined,” J. App. Meck. 2, A67, 
A68 (1935). 

M. M. Frocht, ‘“‘Photoelastic studies in stress concen- 
tration,”’ Mech. Eng. 58, 485—9 (1936). 

M. M. Frocht, ‘“‘A rapid method for the determination 
of principal stresses across sections of symmetry 
from photoelastic data,’’ J. App. Mech. 5, 24-28 
(1938). 

M. M. Frocht, “On the optical determination of 
isopachic stress patterns,’ Proceedings of the 5th 
International Congress for Applied Mechanics, 
Cambridge, Massachusetts (1938). 

T. H. Frost and K. F. Whitcomb, ‘‘The stresses in 
rotating disks,’”’ Trans. Am. Soc. Mech. Eng. 53, 
1-11 (1931). 

A. Goetz, “A modified optical arrangement for 
photoelastic measurements,” Rev. Sci. Inst. 5, 84 
(1934). 


. J. N. Goodier, ‘“‘An analogy between the slow motions 


of a viscous fluid in two dimensions and systems of 
plane stress,” Phil. Mag. 17, 554 (1934). 

S. G. Hall, ‘‘Determination of stress concentration in 
screw threads by the photo-elastic method,” 
University of Illinois Engineering Experiment Sta- 
tion Bulletin No. 245 (1932). 

F. C. Harris, “The photo-elastic constants of glass 
as affected by high temperatures and by lapse 
of time,’’ Proc. Roy. Soc. London 106, 718-23 
(1924). 


291 





F. C. Harris and B. R. Seth, ‘““The variation of double 
refraction in celluloid with the amount of permanent 
stretch at constant temperature and at different 
temperatures,’” Proc. Phys. Soc., London 48, 477- 
487 (1936). 

M. Hetényi, ‘‘Photoelastic studies of three-dimen- 
sional stress problems,” Proceedings of the 5th 
International Congress for Applied Mechanics, 
Cambridge, Massachusetts (1938). 

M. Hetényi, ‘“The fundamentals of three-dimensional 
photoelasticity,”” J. App. Mech. 5, A1l49-A155 
(1938). 

M. Hetényi, ‘“Photoelastic stress analyses made in 
three dimensions,’’ Machine Design 10, 40-1 (1938). 

W. Herzog and G. Szivessy, ‘“‘Uber die Messung von 
Phasendifferenzen mit einer empfindlichen Methode 
von uneingeschriinktem Messbereich,’’ Physik. 
Zeits. 38, 129-133 (1937). 

Rk. Hiltscher, ‘‘Polarisationsoptische Untersuchung des 
raumlichen Spannungszustandes im konvergenten 
Licht,’’ Forschung auf dem Gebiete des Ingenieur- 
wesens 9, 91-103 (1938). 

W. H. Himes and R. V. Baud, ‘‘The use of motion 
pictures to illustrate the generation of involute 
gearing and the use of motion pictures and polarized 
light as a method of studying gear stresses,”’ Iron 
and Steel Engineer 6, 372-380 (1929). 

QO. J. Horger, *‘Photoelastic analysis practically ap- 
plied to design problems,"’ J. App. Phys. 9, 457-464 
(1938 

B. K. Hough, Jr., “Stability of embankment founda- 
tions,” Proc. Am. Soc. Civil Eng. 63, 1340-1357 
(1937). 

A. W. Hull and E, E. Burger, ‘Strain analyser for 
glass seals,’” Rev. Sci. Inst. 7, 98-100 (1936). 

H. Jehle, ‘Polarisationsoptische Spannungsunter- 
suchungen an einer Schrauben-verbindung und an 
einzelnen Gewinde Zahnen,” Forschung auf dem 
Gebiete des Ingenieurwesens 7, 19-30 (1936). 

A. Jona, “La foto-elasticita e sua applicazione allo 
studio delle strutture iperstatiche,"’ Aerotecnica 15, 
180-193 (1935). 

J. Kuno, “Law of photoelastic extinction,” Phil. Mag. 
12, 503-511 (1931). 

J. Kuno, “Law of photoelastic extinction,” Phil. Mag. 
13, 810-824 (1932). 

J. Kuno, ‘Relation between strain and photoelastic 
effect in phenolite,’’ Phil. Mag. 16, 353-362 (1933). 

J. Kuno, ‘Practical methods of determining the 
coefficient of photoelastic extinction,’’ Phil. Mag. 
19, 457-466 (1935). 

J. Kuno, “Effect of age of phenolite on the coefficient 
of photoelastic extinction,” Phil. Mag. 23, 63-64 
(1937). 

A. Kuske, ‘‘Das Kunstharz Phenolformaldehyd in der 
Spannungsoptik,’’ Forschung auf dem Gebiete des 
Ingenieurwesens 9, 139-149 (1938). 

A. Kuske, ‘“‘Berechnung von Fehlern Infolge von 
Vorspannungen in der Spannungsoptik,” V. D. I. 
Zeits. 82, 1455-8 (1938). 


M. v. Laue, ‘Einsatz iiber Eigenspannungen,” Sit, 
der preussischen Akad. der Wissenschaften Ber|j, 
21, 377-382 (1931). 

G. H. Lee and C. W. Armstrong, “Effect of tempera 
ture on physical and optical properties of photo. 
elastic materials,” J. App. Mech. 5, 11-12 (1938). 

H. Liebmann, ‘Die angenaherte Ermittlung har. 
monischer Funktionen.  Sitz. der math.-phys 
Klasse der Bayer. Akad. Miinchen, p. 385 (1918 

G. Mabboux, “Applications de la photoélasticimétrie 
a l'étude des ouvrages en béton,” Rev. d'Optique 11, 
501-507 (1932). 

W. MacGregor, “Concerning the distribution 
stress in.a laterally compressed strip,”’ Physics 5 
140-145 (1934). 

L. Malavard, ‘“L’analogie electrique comme méthode 
auxiliaire de la photo-¢lasticité,”” Comptes rendus 
Paris 206, 38-39 (1938). 

H. B. Maris, ‘“Photoelastic investigations of th 
tensile test specimen, the notched bar, the shi 
propeller strut, and the roller path ring,” J. Opt 
Soc. Am. 15, 194-200 (1927). 

H. E. L. Martin, “Some experiments with reinforce 
materials examined by aid of plane polarized light 
Trans. Liverpool Eng. Soc. 36, 59-98 (1915). 

Tr. V. Matthew, “Stress distribution in irregula: 
sections using photo-elastic models,”’ J. Roy. Tect 
College Glasgow 4, 121-134 (1937). 

J. C. Maxwell, “On the equilibrium of elastic solids, 
l'rans. Roy. Soc. Edinburgh 20, 87 (1850). 

G. Mesmer, ‘Vergleichende spannungsoptische Unter- 
suchung und fliess Versuche unter konzentrierte 
Druck,” Zeits. f. Tech. Mech. Thermodynamik, 
V. D. I. 1 (1930). 

A. Mesnager, ‘‘Contribution a l'étude de la déforma- 
tion élastique,’’ Ann. des Ponts et Chaussées 4, 128 
190 (1901 

A. Mesnager, “Détermination compléte sur 
modéle réduit des tensions qui se produiront dans u 
ouvrage,”” Ann. des Ponts et Chaussées 16, 133-16 

1913). 

A. Mesnager, “Sur la détermination optique ces 
tensions intérieures des solides A trois dimensions 
Comptes rendus, Paris 190, 1249 (1930). 

H. Meyer, ‘“Spannungsoptische Untersuchung ebener 
Schwingungsvorgiinge,”” Ingenieur-Archiv 7, 273 
293 (1936). 

H. Meyer and F. Tank, “Uber ein verbessertes 
elektrisches Verfahren zur Auswertung der Gleic! 
ung A@=0 und seine Anwendung bei photoela> 
tischen Untersuchung,” Helv. Phys. Acta 8, 315-3!/ 
(1935). 

R. D. Mindlin, “A reflection polariscope for photo 
elastic analysis,” Rev. Sci. Inst. 5, 224-228 (1934 

R. D. Mindlin, “Analysis of doubly refracting 
materials with circularly and elliptically polarize 
light,” J. Opt. Soc. Am, 27, 288-291 (1937). 

R. D. Mindlin, ‘Distortion of the photoelastic fringe 
pattern in an optically unbalanced polariscope, / 
App. Mech. 4, A170—-172 (1937). 


JOURNAL OF APPLIED PuHysIcs 





Mindlin, “Optical aspects of three-dimensional 
elasticity,” Eighth Semi-Annual Meeting, 
ern Photoelasticity Conference, N. Y. (1938). 
lueller, “Theory of photoelasticity in amorphous 
s,"’ Physics 6, 179-184 (1935). 
‘lueller, ‘““Theory of the photoelastic effect of 
c crystals,”” Phys. Rev. 47, 947-957 (1935). 
‘lueller, “Light diffracted by supersonic waves in 
s.’ Phys. Rev. 52, 223-229 (1937). 
\lueller, ‘‘Determination of elasto-optical con- 
ts with supersonic waves,” Zeits. f. Krist. 99, 
141 (1938). 
\l. Murray, ‘“‘Some aspects of dynamic problems 
yhotoelasticity,’’ Eighth Semi-Annual Meeting, 
astern Photoelasticity Conference, N. Y. (1938). 
11. Mussmann, “Die elastische Nachwirkung und ihre 
isammenhang mit dem optischen Nachwirkung,” 
d. Physik 31, 121-144 (1938). 
McGivern and H. L. Supper, ‘‘Membrane 
ogy supplementing photoelasticity,”’ J. Frank. 
217, 491-504 (1934) 
McGivern and H. L. Supper, “A membrane 
ogy supplementing photoelasticity,’’ Trans, 
Soc. Mech. Eng. 56, 601-604 (1934). 
|. G. MeNally and S. E. Sheppard, “‘Rate of deforma- 
1 of cellulose nitrate-camphor films under static 
resses,”” J. Phys. Chem. 35, 2498-2507 (1931). 
Masataka Nisida, “On the physical properties of 
toelastic material ‘Phenolite,’”’ Sci. Pap. Inst. 
hys. Chem. Research 22, 269-283 (1933). 
H. Neuber, ‘“‘New method of deriving stresses 
iphically from photoelastic observations,” Proc. 
Soc., London 141, 314-24 (1933). 
Hl. Neuber, “Exact construction of (o,;+02) network 
photoelastic observations,’ Trans. Am. Soc, 
Mech. Eng. 56, 733-737 (1934). 

(Gs. Oberti, “Studi sul comportamento statico di archi 
lari considerati come elementi di dighe a volta,” 
gia Ellettrica 13, 578-584 (1936). 
pel, ‘‘Polarisationsoptische Untersuchung raum- 

spannungs- und dehnungszustinde,” Fors- 
g auf dem Gebiete des Ingenieurwesens 7, 
1936). 
ppel, “‘Photoelastic investigation of three- 
nsional stress and strain conditions,” Nat, 
ry Com. for Aeronautics—Technical Memo. 
$24 (1937). 
, ‘Das polarisationsoptische Schichtverfahren 
\lessung der Oberflachenspannung am_ bean- 
ten Bauteil ohne Modell,” V. D. I. Zeits. 81, 
1937). 
eterson and A, M. Wahl, ‘‘Fatigue of Shafts at 
Members with a Related Photoelastic 
s,"’ J. App. Mech., Trans. Am. Soc. Mech, 
7, Al-All1 (1935) 
eterson and A, M. Wahl, ‘‘Two- and three- 
onal cases of stress concentration and 
son with fatigue tests,” J. App. Mech. 3 
22 (1936). 


, MAY, 1939 


H. Poritsky, ‘‘Analysis of thermal stresses in sealed 
cylinders and the effect of viscous flow during 
anneal,” Physics 5, 406-411 (1934). 


. G. M. Pugno, “Studio sperimentale sugli stati piani di 


tensione,” Atti della Accad. della Science di 
Torino, 71, 55-71 (1935). 


. G. M. Pugno, “Sulla rappresentazione degli stati piani 


di tensione dedotta mediante l’esperimento,” Atti 
della Acad. della Science di Torino, 71, 110-115 
(1935). 
Ravilly, Publications Scientifiques du Ministere 
d l'Air, Paris, No. 120 (1938). 


. J. J. Ryan and L. J. Fischer, ‘‘Photoelastic analysis of 


stress concentration for beams in pure bending with 
central hole,”’ J. Frank. Inst. 225, 513-526 (1938). 


. Ch. Sadron, “Sur une nouvelle méthode optique 


d’exploration d’un champ de vitesses bidimensional,”’ 
Comptes rendus, Paris 197, 1293-1296 (1933). 

Ch. Sadron and E. D. Alcock, ‘“‘An optical method for 
measuring the distribution of velocity,’”’ Guggenheim 
Aeronautical Laboratory, California Institute of 
Technology, Publication No. 44 (1934). 

M. Salvadori, ‘“Alcune richerche fotoelastiche nel 
Laboratorio dello University College di Londra,” 
Annali dei Lavori Pubblici 7, 865-879, 1050-1091 
(1934). 

M. Salvadori, ‘“‘Sulle tensioni interne nelle ruote 
motrici delle locomotive,”’ Ricerchi d’Igegneria, 3, 
113-118 (1935). 


. G. H. Shortley and R. Weller, ‘‘The numerical solu- 


tion of Laplace’s equation,” J. App. Phys. 9, 334- 
348 (1938). 


. J. E. Soehrens, R. T. Cass and J. E. Sower, ‘Elimi- 


nating local stress concentrations in photoelastic 
model,”’ Civil Eng. 6, 594-5 (1936). 

A. G. Solakian, “Stresses in transverse fillet welds by 
photoelastic methods,”’ J. Am. Welding Soc. 13, 22- 
29 (1934). 

A. G. Solakian, ‘“‘A new photoelastic material,”” Mech. 
Eng. 57, 767-771 (1935). 

A. G. Solakian, ‘“‘Photoelastic models with cemented 
elements,” Photoelastic J. 1, 14-17 (1938). 

A. G. Solakian and G. B. Karelitz, ‘‘Photoelastic 
study of shearing stresses in keys and keyways,” 
Trans. Am. Soc. Mech. Eng. 54, 97-123 (1931). 

H. Souillot, “La photo-élasticimétrie méthode G. 
Mabboux dans l'étude des déformations des 
ouvrages en béton armé,” Tech. des Travaux 19, 
157-9 (1935). 

M. Suquet, ‘‘Recherches photoélasticimétriques sur 
les dalles nervurées poursuivies au laboratoire de 
l'Ecole National des Ponts et Chaussées,’ Ann. des 
Ponts et Chaussées,” 104, 34-39 (1934). 

Tank, “Die Tatigkeit des photo-elastischen 
Laboratoriums der Eidg. Technischen Hochschule 
1927 bis 1933,’” Schweizerische Bauzeitung 104, 45- 
48 (1935). 

F. Tank, ‘‘Eine neue Methode der interferometrischen 
Spannungsmessung,’’ Helv. Phys. Acta 9, 611-616 
(1936). 


293 





lank, R. V. Baud and E. Schiltknecht, ‘Die neuen 
Einrichtungen des photoelastischen Laboratoriums 
an der EKidg. Techn. Hochschule und an der Eidg. 
Materialprufungsanstalt,’’ Schweizerische Bauzeit- 
ung 109, 249-52 (1937). 

lank, R. V. Baud Schiltknecht, ‘‘Les 
nouvelles installations du laboratoire de photo- 


and E. 


élasticimétrie de I'Ecole Polytechnique de Zurich,” 
Genie Civil 111, 482-3 (1937). 
P. Tarbés, ‘Méthode pour |’étude de la soudure et de 
e biré- 


224 


la dilatation des verres. Compensature 


Paris 


d 
195, 222 


fringence,”’ Comptes rendus, 

1932 
Pas: Bey Bs 
meésure de la biréfringence des verres d’optique,” 
d’Optique 8, 59-69 (1929). 


lemmermann and L. Blanjean, “ 


Tardy, ‘‘“Méthode pratique d’examen de 


Rev 
Un probléme de 
de photoélasticité,” 
1938). 


matériaux et 
Ossature Metallique 7, 132-7 


\V. Tesa¥, ‘Photoélasticité méthode purement optique 


résistance des 


pour déterminer les déformations d’épaisseurs des 
d’Optique 11, 97-104 (1932). 


lesafr, “Détermination expérimentale des tensions 


modeles, etc.,’” Rev 


dans les extrémités des piéces prismatiques munies 


Assoc. Int. des Ponts et 


1932). 


experimentale 


d'une semiarticulation,” 
Charpentes 1, 497-505 
Vv. Tesaf, 


produites dans une poutre par des charges con- 


“Etude des contraintes 
centrées,”” Int. Verein. der Brucken-und-Hochbau 4, 
543-570 (1936). 

Tesaf, “Les procédés photoélasticimétriques et 


leur application aux barrages-voutes,”’ Ann. des 
Ponts et Chaussées 107, 627-662 (1937). 

W. E. and A. T. McPherson, ‘‘Photo- 
elastic properties of soft vulcanized rubber,” Nat. 
Bur. Stand. J. Research 13, 887-896 (1934). 

W. E A. Wood, ‘‘Photoelastic 


determination of stresses around a circular inclusion 
Nat. Bur. Stand. J. Research 20, 393- 


Phibe xleau 


Thibodeau and L. 


in rubber,” 
409 (1938), 
Thouvenin, “Application de la photoélasticité a 
l'étude des percussions,’ 
769-771 (1935). 
Kk. Timby and I. G. Hedrick, Jr., ‘‘Photoelastic 
analysis broadened,”’ Eng. News-Record 121, 179-81 
1938). 
Z. Tuzi, 
elasticity,” 
79-96 


"Comptes rendus, Paris 201, 


“A new material for the study of photo 
Inst 


1927). 


Phys. Chem. Research Japan 7 


168. 


Z. Tuzi, ‘“Photoelastic study on a specimen of three. 
dimensional form,” Sci. Pap. Inst. Phys. Chey 
Research Japan 7, 97-103 (1927). 

Z. Tuzi, “On the experimenta 
methods in photoelasticity,” Sci. Pap. Inst. Phys 
Chem. Research Japan 12, 21-69 (1929), 

Z. Tuzi, “Stress distribution in an angle plate,” & 
Pap. Inst. Phys. Chem. Research Japan 16, 140-145 
(1931). 

. Z. Tuzi, ‘Photographic and kinematographic study of 
photoelasticity,”’ Sci. Pap. Inst. Phys. Chen 
Research Japan 8, 247-267 (1928). 

Z. Tuziand M. Nisida, ‘*Photoelastic study of stress 
due to Sci. Pap. Inst. Phys. Chem, 
Research Japan 26, 277-309 (1935). 

Z. Tuziand M. Nisida, ‘‘Photoelastic study of stresses 
due to impact,”” Phil. Mag. 21, 448-473 (1936). 

Z. Tuzi and M. Nisida, “On a new apparatus for 


development of 


impact,” 


photo-elasticity-stress-type quarter-wave-plate of a 
large field,’ Sci. Pap. Inst. Phys. Chem. Researc! 
Japan 31, 99-107 (1937). 

RK. von Mises, ‘Das Verhalten der Hauptspannunge: 
in der Umgebung einer Verzweigungstelle,” Zeits. { 
angewandte Math. Mech. 18, 74-76 (1938). 

R. von Mises, ‘‘Uber den singuliren Punkt zweiter 
Ordnung im ebenen Spannungsfeld,”’ S. Timoshenk 
60th Anniversary Volume (Macmillan Company 
New York), pp. 147-154. 

R. W. Vose, “An application of the interferometer 
strain gage in photoelasticity,” J. App. Mech. 2 
A99-A102 (1935 

A. M. Wahl and R. Beeuwkes, ‘Stress concentratior 
produced by holes and notches,”’ Trans, Am. So 
Mech. Eng. 56, 617-625 (1934). 

E. E. Weibel, ‘“‘Studies in photoelastic stress detern 
nation,”’ Trans. Am. Soc. Mech. Eng. 56, 637-658 
1934). 

E. Weibel, ‘Application of soap-film-studies t 
photoelastic stress determination,”’ Assoc. Int. des 

Mémoires 3, 421-38 (1935 


E. E. Weibel, ‘‘Thermal stresses in cylinders by th 


Ponts et Charpentes 


photoelastic method,”’ Proc. 5th Int. Cong. for Apt 
Mech., Cambridge, Massachusetts, 1938. 
G. Wertheim, ‘‘Mémoire 


temporairement 


sur la double refract 

produite dans les corps 1% 

tropes. Ann. de Chimie et de Phys. 60, 15 
1854 

H. E. Wessman, ‘‘New universal straining frame a 
photoelastic research,” Civil Engineering 8, 61+ 

1938). 


JOURNAL OF APPLIED PHYSICS 





The Application of Hardening Resins in Three- 
Dimensional Photoelastic Studies 


By M. HETENYI 


Westinghouse Research Laboratories, East Pittsburgh, Pennsylvania 


ning resins at elevated 


temperatures 


(ENT photoelastic tests made with hard- 


| 
R’ 
furnish a striking demonstration of the current 
sheory for the molecular structure of these resins 
vhile the results of such tests lay the foundation 

a new line of photoelasticity dealing with 
three-dimensional states of stress distribution. 

lhe fundamentals of the new method can be 
easily demonstrated by subjecting a bar of such 
hardening resin, phenol-formaldehyde resin for 
stance, to some simple loading at a temperature 
f about 110°C. We can observe that while the 
material becomes very flexible at this tempera- 
that is it 
issumes its final deformation in a few seconds 


it remains perfectly elastic, 


upon application of the load, while upon unload- 
ng the piece the original shape is also recovered 
most instantaneously. If we leave the loading 

the test piece at the elevated temperature 

| by slow cooling return to room temperature, 
ter the removal of the loading we obtain a 
permanently deformed sample which, however, 
xhibits in polarized light an  isochromatic 
ige pattern corresponding to an elastic state 
i stress. This observation, that a permanently 


deformed body can exhibit an 


elastic stress 
tribution, is quite unusual in itself, though 
that the obtained test piece has also 
curious properties. We will observe 

eful sawing or machining of the sample 

ot disturb its isochromatic pattern and the 

| fringes will be preserved in every small 

it out of the test piece. Thus we can 

that the annealing has 
us a_ solidified, form of 
iormation, which is completely pre- 
ugh the whole bulk of the test piece. 
that 


procedure 
“frozen” 


the isochromatic lines are due 
| the material, their preservation will 
ibly to the same particular properties 


ucture of hardening resins which 
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produced the obtained frozen form of the elastic 
strains. Accepting the validity of the above 
observations they will evidently provide means of 
investigating photoelastically three-dimensional 
states of stress distribution. Subjecting any 
three-dimensional model to the above-mentioned 
annealing procedure we obtain a sample where 
the elastic strains produced by the loading at the 
elevated temperature are preserved in a frozen 
form. Since cutting of the test piece does not 
disturb this solidified state, we will be able to 
cut out thin slices of the sample and each of the 
slices will show in polarized light the stress 
pattern existing in the original three-dimensional 
piece in the corresponding plane. 

All the above-mentioned experimental ob- 
servations can be well explained on the basis of 
the current molecular theory of heat-hardening 
resins. We will discuss therefore briefly this 
molecular theory! taking the phenol-formalde- 
hyde resin as a prototype for our material. 

The phenol-formaldehyde resin is a result of 
interaction between phenol and formaldehyde. 
Reacting these two chemicals at suitable temper- 
ature, we obtain a clear resin-like material which 
is in a solid state at room temperature, can be 
fused, however, and will dissolve in solvents such 
as alcohol and acetone. The assumed chemical 
reaction scheme of phenol-formaldehyde resins 
is shown in Fig. 1. In the first (A) stage the 
material is fusible and soluble and its molecules 
adhere only by the secondary bonds of the van 
der Waals forces. The benzene nuclei, however, 
are assumed to possess three active spots (shown 
by crosses in Fig. 1) which gives a possibility for 
a three-dimensional molecular growth. 
heating this A the assumed 
react with one 
another, creating thus a strong primary bond 


Upon 
stage material, 
chemically active spots will 


among them. As this reaction proceeds a strong 
molecular network will be built up throughout 
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the whole bulk of the resin. The molecules of this 
network will be then in primary bond with one 
another, being thus in an infusible and insoluble 
state. Because of some particular features of the 


OH OH 
x x HOH.,C CH,OH 
O + 3CH,(OH) > -_ © 2 


x CH30H 
PHENOL FORMALD ¢ HO PHENOL ALCOHOL 


OH 


ie 


H 
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eo 


He) 
1 


5 “Ge @ OH 


HOH, CHa 
OH 
® Bond -CHa— CH>— WC 
OH 


Fic. 1. Areaction scheme of the phenol-formaldehyde resins 
(hardening type). 


CH,0H 


structure, however, not all the possible chemi- 
cally active spots will be able to react and within 
the globular micelles of the 


network in the 


three-dimensional 


resin there will be material 


which is still in the fusible and soluble state. 
This 
through the B (partly fusible and soluble) stage 
to the 


as polymerization. The 


process, as the A stage material arrives 
C (infusible and insoluble) stage is termed 
synthetic resins of the 


hardening type used in photoelastic tests are 


largely in this latter stage, where their molecular 
structure can be depicted as shown in Fig. 2. 
This state is termed isogel, because there is no 
difference in the chemical composition between 
the fusible and infusible part of the material. 
When a piece of such resin is heated, the A 
stage material within the cells will be gradually 
fused, becoming liquid at a certain temperature. 
At that point the only load-carrying element of 
the material will be the 
network throughout the resin. 
The strength and elasticity of this network will, 
the 
the A 


we load such 


its skeleton given by 
strong primary 
however, not be affected by 
melt 


temperature 
able to 
When 


elevated temperature, 


which was stage filling 
a piece at the 
the load will be carried 


completely by the primary network and conse- 


material. 


quently the deformation in all its characteristics 
will be purely elastic. If we have now the load 
on the sample and lower the temperature, the 


fusible A stage material will return gradually to 
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a solid state again, and surrounding with jt 
large mass the elastically deformed skeleton oj 
the material will keep the piece in its deformed 
position even when the load is removed. In this 
between the 
elastically deformed network and the restraining 


mass of the filling material. Since the equilibrium 


state there is an equilibrium 


is of molecular dimensions, we have the explana. 
tion of why the deformation together with the 
accompanying birefringence can be preserved in 
any arbitrarily small piece cut out of the original 
sample. 

In order to obtain quantitative results of the 
properties of some hardening resins at elevated 
temperature the writer has conducted simple 
bending and Bakelite and 
Marblette, both of the heat-hardening type 
Fig. 3 shows the results of some bending tests 
with BT-61-893 Bakelite 
anhydride, 


tension tests with 


This material js 
but 


bars.” 


glycerine phthalic being also 


MICELLE 
(COMPARE TO IZ) 
IN FIGI 


“SOLID” MICELLES. 
“LIQUID” MICELLES. 


A hardening resin in the isogel state (according | 
R. Houwink in the British Plastics, 1935). 


three-dimensional polymer its molecular struc- 
ture has the same characteristics as the ones 
described above for phenol-formaldehyde resins. 
The dimensions of the test bars and the loading 
were the same in each case, only the applied 
temperature was different as is shown in Fig. 3 
It is seen that while the speed of deformation 
depended on the applied temperature the final 
amount of deformation was the same in each 
Above 110°C the bar assumed its final 
deformation in a few seconds after the application 
of the load; the deformation then remained the 
same and no creep phenomena could be observed. 
At 90°C it took about 8 hours for the sample 
arrive to the same amount of deformation whic! 
it assumed almost instantaneously at 120°C. 
On the basis of the above-described theory the 


case. 


. . 2 ole 
structure of hardening resins can be symbolicall) 
represented by a parallel coupling of an elasti¢ 
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ruc- 
ones 
ins. 
ding 


plied 


ation 
final 
eacl 
fina! 
ation 
| the 
rved 
le to 
which 


he primary network) and a dashpot 
le filling material) as shown in Fig. 4. 
ning various viscosity coefficients de- 
on temperature to the viscous material 
‘1 the dashpot, the curves of Fig. 3 could be 
reproduced by such a model. At the moment 
when an external force P is applied to this 
model, the load will be carried first mainly by 
the viscous component, on whose gradual yield- 
hen the elastic spring will come in action, 
isymptotically approaching the limit where all 
the load will be carried by the elastic part. The 
transition will depend on the viscosity coefficient 
of the material in the dashpot. For phenol- 
formaldehyde resin in the fusible stage the 
viscosity coefhicient decreases approximately to 
1 10 of its value by each 10°C rise in tempera- 
iure.! This means that even if the final defor- 
takes place at 120°C in 1 to 2 seconds, 
room temperature the time required for the 
sample to arrive to the same amount of defor- 
mation under the same loading would be of the 
ler of ten thousand years. 

\t the annealing temperature, where the 
resistance of the large fusible part is practically 
eliminated, elastic deformation of the 
samples can be observed. 


great 
Bakelite material 


ing at room temperature a modulus of 
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e deflection curves for Bakelite at elevated 
temperature. 


710,000 Ib./sq. in. and an optical 
such that 1 isochromatic fringe =85 
per inch will give at the annealing 

E=1080 |b./sq. in. and 1 fringe 
sq. in. per inch which means a 
f the modulus of elasticity by 640 

in the optical sensitivity of the 
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Fic. 4. Mechanical analogy to the structure 
of hardening resins. 


material a 26-fold increase was observable. This 
means that after the sample has been annealed 
under the load, upon removal of the loading 
1/640 of the deformation and 1/26 of the number 
of isochromatic bands will be recovered from 
their respective values produced at the elevated 
temperature. 

The experimental fact that a “frozen’’ sample 
preserves its fringe pattern even if it is cut in 
places is illustrated in Fig. 5. The upper picture 
shows the fringe pattern in a Bakelite bar after 
it has been annealed under a concentrated load 
applied at the center of the span. The lower 
picture shows the same bar after it has been 
cut in at two places. During the cutting process 
some precautions should be taken not to warm 
up the sample by frictional heat produced by 
such action. An ordinary, not too fine, mechani- 
cally driven jig saw with use of abundant 
lubricating oil assures good results. 

The application of the method to the solution 
of three-dimensional stress problems is illustrated 
in Fig. 6. The upper picture shows the Bakelite 
model, a shaft with a keyway which has been 
subjected at 115°C to simple bending. After 
the annealing process thin slices were cut out of 
the test piece which were then examined in 
polarized light as is done in the usual photoelastic 
method. The fringe pictures show the amount 
of stress concentration (that is the ratio of the 
peak stress to a properly chosen nominal stress) 
at two ends of the keyway. At the sled-runner 
end the maximum stress was observed in the 
vertical central slice taken through the keyway 
which gave a stress concentration factor k= 1.38 
while the peak stress at the profiled end was 
observed in the horizontal surface slice giving 
k=1.79. Both of these concentration factors are 


297 











ic. 5. Illustration showing that careful sawing does not 


disturb appreciably the frozen fringe pattern. 


in good agreement with the results derived from 
fatigue tests by R. E. Peterson.® 

This new photoelastic method can be applied 
in all cases where the test piece is subjected to 
constant loads. This is the case with rotating 
disks when the centrifugal loading is constant 
as long as the speed of rotation of the disk is 
the the this 
consideration the method has been used by the 


maintained basis. of 


same. On 
writer to determine the distribution of centrifugal 
stresses in rotating disks. The procedure con- 
sisted in rotating the test piece in an oven at 
elevated temperature (115°C) for about one hour 


and then slowly cooling it down while main- 











b) 


Fic. 6. Problem of shaft with keyway. (a) View of Bakelite test piece. (b) Ver- 
tical slice through sled-runner keyway stress concentration factor k = 10.00/7.25 
Horizontal slice through profiled keyway k =5.9/3.3 =1.79. 


1.38. (c) 
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taining the same speed of rotation. This way , 


sample was obtained which exhibited a solidified 
form of the deformations and _ birefringence 
produced by 


the centrifugal stresses at the 


annealing temperature. 


iB 


lic. 7. Dimensions of test piece. 


The dimensions of one of the test pieces, a 
disk with two holes, are shown in Fig. 7. The 


disk 


during the annealing process, and the resulting 


was rotated at 


a speed of 3100 rpm 
frozen sample gave then in polarized light an 
isochromatic fringe pattern as shown in Fig. 8 
the 0.40-inch thick material 
(Bakelite 61-893) of the sample it was found that 
one isochromatic line corresponds to 8.27 Ib. 


By calibrating 


sq. in. difference in principal stresses (¢;—#:). 
Knowing this from Fig. 8, the difference o! 
principal stresses could be determined at ever 
point of the test piece. In order to derive values 
for each of the principal stresses, lateral extensio 
measurements were made 0! 
the sample determining the 
change of thickness produced 
by the centrifugal forces. I 
is known that the change 0 
thickness At=—(ut/E)(oit% 
is proportional to the sum 0! 
principal stresses where # © 
Poisson’s ratio and E is the 
modulus of elasticity. Measu 
ing the thickness of the te* 
piece before test and after the 


(c) 


annealing process, these + 
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ere determined along section A—A and 

several points. Since at the annealing 
temperature E=1080 Ib./sq. in. and w=1/2 
due to one fringe value (8.27 lb./sq. in.) a 
t= (0.50X0.400"/(1080 Ib./sq. in.))8.27  Ib./ 
sq. in. 0.00165” lateral extension of the 0.400- 
‘nch thick sample will be produced. With a ball 
pointed micrometer’ reading 0.0001 inch, this Af 
value can be measured with about 6 percent 
curacy, but for higher values the 
iecuracy will be proportionally increased. This 
wav, for both of the symmetry sections A—A 
nd B-B, the values for each of the principal 
tresses could be determined. The final results 


stress 


re shown in Figs. 9 and 10, respectively. 


‘. lsochromatic lines due to centrifugal stresses in a 
disk with two holes. 


de 0 For 


: he solid A—A section the stresses were 
g in 


omputed also analytically.’ The good agreement 
found between the theoretical and experimental 
for section A-—A is shown in Fig. 9. 
on B-B no theoretical solution is 
present. Considering as concentra- 
t k for the present case the ratio of the 
stress at the hole in section B—B, to 
lum stress which would be produced 
section (A—A) without the hole, a 
95.1/49.6=1.92 can be derived. 
rged picture of the isochromatic 
und the hole is shown in Fig. 11. 


duced 
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mMaximut 
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Fic. 9. Principal stresses in section A-A (Fig. 7). 


Using plane polarized light, the isoclinic lines 
can be determined at any part of the frozen 
sample, from which data then the directions of 
the principal stresses can be derived. This has 
been carried out on a portion of the sample 
around the hole and the results are shown in 
Fig. 12. The results of this test show that from 
the frozen sample a complete knowledge of the 
95.1 


EXPERIMENTAL DATA 


FROM ISOCHROMATIC 
PATTERN (Op -6;) 





FROM LATERAL EXTENS. 
MEASUREMENTS (6, + 6; ) 
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Fic. 10. Principal stresses in section B-B (Fig. 7). 
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Fic. 11. Isochromatic lines. Detail of Fig. 8. 


SOCLINIC LINES STRESS TRAJECTORIES 


Fic. 12. Direction of principal stresses 


stress system, including the values of both of 
the principal stresses and their direction at any 
required point, can be derived. 

The annealing process was also applied to 
rotating disks of various other designs. Fig. 13 
shows the frozen fringe pattern obtained in a 
disk with 8 holes. The samples used in these 
tests were of constant thicknesses but the same 
method can be applied to disks of variable 
thickness. Since with some approximations 
gravity forces can be replaced by centrifugal 
loading, the above-described method can be also 
applied to determining stress distribution due to 
gravitational forces. 
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The Operational Calculus. IIT 


By Louis A. PIPES** 


Electrical Engineering Department, The University of Wisconsin, Madison, Wisconsin 


’ Solution of Linear Partial Differential Equa- 
tions with Constant Coefficients 


A. Tue VIBRATING STRING 


\s a simple example of the use of the opera- 


onal method of solving partial differential 
quations with constant coefficients subject to 
e-point boundary conditions in one of the 
the 


homogeneous Wave equation in one dimension 


dependent variables, let us consider 


(132) 


07 / Ot? = a*d?u/ dx. 


o take a concrete example, let us consider the 

i vibrating string. In such a case u is the 

teral displacement of the string and a?=T7)/m, 

ere 1] is the tension of the string and m is the 

mass per unit length. Let us consider the string 

ed at the extremities x =0 and x=L, and let 

e disturbance be initiated with a displacement 

| at t=0. We have then, the following 
undary conditions 


u=Qatx=0 and x=L. (133) 


Let us place 


U(x, p)=u(x, t) (134) 


‘use of the transforms of the second 
by Eq. (6). 


obtain the transform of Eq. (132) in 


Ty 1 9 9 Py 


dx*— p*U/a®?=—p*F(x)/a?, (135) 


have made use of the boundary con- 


ditions (133) which depend upon the time. We 
may consider the transformed equation as an 
ordinary differential equation in x since ~ enters 
only as a parameter. The method of variation of 
parameters may be used to solve this ordinary 
differential equation with constant coefficients 
by assuming a solution of the form 


U=A(x) cosh (px/a)+ B(x) sinh (px/a), (136) 


where we may impose the restriction 


A'(x) cosh (px/a)+B’(x) sinh (px/a)=0 = (137) 


on the derivatives of the functions A(x) and 
B(x). Substituting (136) into (135) we obtain 


A’(x) sinh (px/a)+B’(x) cosh (px/a) = 


—pF(x)/a. (138) 


Solving (138) and (136) for A’(x) and B’(x) and 
making use of the boundary conditions (133) on 
x, we obtain 


A(x) = / (p/a)F(v) sinh (pv/a)dv (139) 
and 
- 
B(x) = —coth (pL/a) | (p/a) F(v) sinh (pv/a)dv 
aL 
+] (p/a)F(v) cosh (pu/a)dv. (140) 


-< 


These values of A and B are substituted into 
Eq. (125), and we obtain, after some reductions, 


( (pa) F(x) sinh (pv/a) sinh p(L—x)/adv 


sinh (pL a) 





*L (p/a) F(v) sinh (px/a) sinh p(L —v)/adv 
+| 
J, sinh (pL /a) 


(141) 


last of a series of three papers. The first and second parts appeared in the March and April issues, 


irvard University, Graduate School of Engineering, Pierce Hall, Cambridge, Massachusetts. 
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Before the integration with 


respect to v, we may first obtain the function 


carrying out 


u(x, t) by means of the basic integral 


* 1 2 — 
u(x, t)= 
2rie 


e”'U(x, p)dp/p. (142) 


e— ta 


kor this purpose, we must compute the 


residues of the expression 
F(v) e? sinh (pv/a) sinh p(L—x) a 


- —, (143) 
a sinh (pL /a) 


The zeros of the denominator, and hence the 


(p,a) sinh (pv/a) sinh p(L—x) a » 


— 


sinh (pL a) 


By symmetry, the second term of (141) has 
the same transform and we may, accordingly, 
combine the two integrals and write 


- ol 
u(x, t)=>> (2/L) F(v) sin (rrv/L)dv 


Xsin (rrx/L) cos (rmat/L). (147) 


This is the solution obtained by the classical 
method of normal coordinates. We note that if 
we place {=0 and hence u(x, t)= F(x) we have 


I 


F(x) = > (2/L) [ F(v) sin (rmv/L) 


1 J 


Xsin (rrx/L)dv, (148) 
the usual Fourier sine series expansion for an 
arbitrary function F(x) valid in the interval 
0O-L). 

The procedure followed in this example is 
typical of the way in which partial differential 
equations are solved by using the Laplacian 
transform theory. We note that the boundary 
conditions on the independent variable, ¢t, are 
given at the single point t=0. This makes the 
the formally 


transformation § of derivatives 


possible. 
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(rix/L) sinh (rixv,L) sinh rix((L—x)/L)e"itet/£ 


=(2/L)> sin (rxv/L) sin (rxx/L) cos (rrat/L). (146 
r=1 








possible poles are at the values of p for which 





sinh (pL /a) =0. (144 





Now a product series expansion for sinh (z) ; 


#/) is 





Fo 


sinh (zs) =2 Il (1+ ; :). (145 
S=1 S*x? 


Hence, sinh (s)=0 at s=rim where r=0, +1 
+2, etc. 

Therefore the integrand has simple poles a 
pL/a=rix where r has the above values. We may, 






therefore, write 
















(arixr LL) cosh (rix)L/a 







B. THe GENERAL TRANSMISSION LINE 





As a more complicated example of the method 





let us consider the problem of determining th 





voltage and current distribution of a finite trans- 


mission line with terminal impedances and 


constant distributed parameters. If i(x, t) denotes 





the current and v(x, ¢) the voltage distributi 





at a point x on the line, the differential equations 





giving the voltage and current along the line ar 






01 ov Ov O1 
L—+Ri=-—, C—+Gr=-——, (4 
at Ox al Ox 







where R and L are the resistance and inductance 





coefficients per unit length (including both sides 
of the line) and G and C, are the leakage con- 


ductance and capacitance coefficients betwee! 







the conductors per unit length. Let us now platt 






v(x, t)=Vi(x, p) and 1(x, t)=1(x, p). 130 





Then on the assumption that the initial curren! 





and voltage distribution of the line at f=" * 
zero, we have the transformed equations 







—dV/dx=(R+Lp)I, 





(151 





—dI/dx=(G+Cp) V. 







Eliminating J between these two equatiol: 
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@?V/dx*=a?V, (152) 


is the operational propagation constant 
lefined by 
a’=(Lp+R)(G+Cp). (153) 


ting the two differential equations, we 


V(x, p) =A e~-*7+A ce", 
(154) 
I(x, P) =A se art 4 a. 


If we use the Eq. (151), we may eliminate the 
nts As and A, and obtain 
V(x, p) =A,e-** +A 2e** (155) 
I(x, p) = A 1€ “- LZytA 2€°X, Zo; 


} eTe Z 
pedance of the line and is defined by 


is the characteristic operational im- 


Zo=((R+Lp)/(G+Cp))'. (156) 


Ao, we 


the and 


ssume an arrangement as shown in Fig. 5. 


lo determine constants A, 

\Ve assume an e.m.f. e(¢) impressed on the 

0 through a terminal impedance Z, 

nd the line closed at x=S through a second 

terminal impedance Z. We have, therefore, the 
oundary condition 


v(0, t) =e(t) —Zy(t) 
v(S, t)=Zat(S, t). 


(157) 


ii we assume the impedances Z, and Zz to be 

nitially without charges or currents at t=0, and 
( k\p)=e(t), we have the boundary conditions 
be satisfied by the transforms 


Zol\S, p)=V(S, p) atx=S (158) 


at x=Q. 


0, p)=E(p) —Z,1(0, p) 

uting the expressions (154) for V and 

e boundary conditions (158), we obtain 

}equations in the unknown constants A, and 
hese may be solved and we obtain the 
xpressions for A; and Ag: 


(1+w)E(p) 


1+w)—(1—q)(1—w) exp (—2aS) 
(159) 


1 
tt 
) | 


-(1—w)E(p) exp (—2aS) 


ations I+. (1+w)—-(1 —q)(1—w) exp (—2aS) 
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where 


g=Zi:/Z, and w=Z2/Z,. (160) 


If these constants are substituted into the Eq. 
(154) for I(x, p) we find, after some reductions 


_ Ep) 
(Zo+Z1) 


{Lexp (—ax)+b; exp (—a[2S—<x]) ]} 
~—- —__—__—__ ——, (161) 


I(x, p) 


(1—0,b2 exp (—2aS) 
where 
by=(Zo—Z)/(Z0+Z3) (162) 


and 


be=(Zy—Z2), (Zo+Z2). 


We then have directly 


1 
(x, 1) =— | 
2rie 


c—tco 


acti 


e”'I(x, p)dp/p. (163) 


Many interesting special cases may be worked 
out directly. We note that if e(t)=A sin (of), 
then E(p)=Awp/(p?+?). Evaluation of the 
integral at the poles p= +iw, gives the steady- 
state solution of the problem due to an impressed 
sinusoidal electromotive force. In general, the 
evaluation of the residues of the integrand to 
obtain the transient terms is a formidable task. 
Considerable simplification is effected in case 
Z,=Z,=Z2 or in case the terminal apparatus 





x0 X=5 











Z, 





























Fic. 5. 


has an impedance that is equal to the character- 
istic impedance of the line. Then 6;=).=0, and 
we have 


I(x, p) =E(p)e-**/2Zp. (164) 
In case we have the further condition that the 
line parameters satisfy the equation 
RC=GL, (165) 

then: 


E(p) exp{| —Rx(C/L)'} exp{ —xp(LC)!} 
I(x, p)= $$, 
Zo (166) 
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Now since e(t)=E(p), we may. use our 


Theorem V, Eq. (10), which states that 
if g(p)=h(t), (10) 


kPg(p)=h(t—k) 
=() 


then e t>k, 
t<k. 
Therefore, 


i(x,t)=0 for t<x(LC)}, 
e(t—x(LC)') exp| —Rx(C/L)*} 


2Z 


a(x, t)= (167) 


for t>x(LC)'. 


In this case there is no current at the point x 
until ¢=x(LC)} 
with a velocity of 1/(LC)! and the current has 


so the disturbance is propagated 


the same form as the voltage impressed at the 
the 
exponential term. Such a line is called a dis- 


sending terminal but is attenuated by 
tortionless line. Many other interesting special 
cases may be worked out such as open-ended 
line, lines short-circuited at the terminals, etc. 

We shall now turn to a procedure which will 
give the solution in the general case. Let us 
return to our general equation for the transform 


of the current distribution along the line. 


exp (—ax)+2 exp [ —a(2S—x) ] 
——- —-——- —. (161) 
[1—b,b2 exp (—2aS) | 
Let us place 
M=Z,)/(Z0+-Z)). 
We may then write 
I(x, p)=E(p)A(p), (169) 
where 
M exp (—ax)+b. exp [—a(2S—x) | 
A(p)=—- ——_—_—__—_——. (170) 


Zo 1 —b,b2 exp (—2aS) 
Expanding A(p) we obtain 
M 
A(p)=— le art Doe 
Zo 
4b, be% 


a(2S—z +b boe a(28+2) 


(AS—2) 4 fh 2hore-@CASta) 4 .-.- +}, (177) 


In obtaining the inverse transform of this 
expression by our basic formula we must have 
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the real part of £ positive along the path of 
integration. Hence the series is absolutely cop. 
vergent and term by term integration is justified 
Therefore the expansion is legitimate. 

The quantity A(p) is the operational equiva. 
lent of the current distribution along the line dye 
to the application of a constant unit electro. 
motive force at the terminal. Such a quantity is 
known as the operational indicial admittance of 
the circuit. The various terms in the series have 
interesting physical interpretation. 

The term exp (—ax)/Z» is the operational 
formula for the current at a distance x of an 
infinitely long line with unit e.m.f. impressed 


directly on the line at x=0. Let 


exp [ —a(2S—x) ]/Zo+des_2(t). (172 


Then exp [ —a(2S—x) |/Z is the operational 
formula for the current flowing at a point (2S—. 
of an infiniteiy long line with a unit  electro- 
motive force impressed directly on the line at 
the terminal. Let 

exp (—ax)/Zo+a,(t). (173 


We now introduce a series of functions, : 


V1, Ue, etc. that satisfy the operational equations 


(174 


Vo= M P), Vo= Mbybe, 
= i 


Vi= Mhbe, U3 [b,b-?, etc. 


Then we may use our basic Theorem III (8b 
that if 
go(p) = ha(t), 


gi(p)=h(t) and 


then 
tf 


1 
} hy(t—u)ho(u)du. 
1 te 


dl 0 


£:\P)g2(p) = 


Hence if we let 


B(x, t)=A\(x, p), 


we have 


d ;** 
t)= [vo(t—u)a,(u) +0, (t—u)des-z 


B(x, t)h= 
dte 0 


+2(t—u)desr2+ +++ ]du. (176 


Now the first term in the above expressic! 
represents the current at a point x of an infinitely 
long line in response to a unit e.m.f. impresse’ 
at x=0 through an impedance Z;. The secon 
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(176 


essi( if 
nitel’ 
ressed 


second 


rysics 


‘erm is a reflected wave from the other terminal 
produced by the irregularity which exists there. 
The ‘ird term is a reflected wave from the 
«ending end terminal, etc. 
© current distribution 7z(x,/) due to the 
ry e.m.f. e(x, f) is obtained by our basic 
m in the form 


d * 
e(t—u)B(u)du. 


dte 0 


i(x, t)= (177) 

[his so-called wave method of solution is most 
laborious but entirely possible and may be 
carried out if the importance of the problem 
iustifies the labor involved. 


C. APPLICATION TO HEAT CONDUCTION 


lhe well-known equation which determines 
the distribution of temperature in a homogeneous 
substance of density D(g/cm*) specific con- 
ductivity K(cal./em sec. °C) and specific heat, 
cal 4 ) is 


dU /dt=h?V?U, (178) 

-['(x, y, 2, t) is the temperature distribution 
in degrees centigrade and h? = K/CD. h?(cm?/sec.) 
is the thermometric conductivity and we have 
issumed that there are no sources of heat in the 
substance. As an application of the operational 
method of solving this equation, let us consider 
the problem of the cooling of a uniform rod 
whose sides are thermally insulated. To fix the 
problem, let us consider the boundary conditions 


is follows: 
r=0, V=Satx=L, V=Sati=0. (179) 


Chat is, the initial temperature of the rod is 
+. Then at time t=0, the end x=0 is cooled 
‘0 temperature zero and maintained at that 
temperature, The end x=L is kept at tempera- 
lhe problem is to find the variation of 

re at other points of the rod. 
ne-dimensional case, the equation for 

rature distribution reduces to 


dU/dt—h?02U/dx? =0. (180) 


l’ and we obtain the transformed 


pV —h?d*V /dx? = pS, (181) 
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where we have made use of the boundary con- 
dition U=S at t=0. If we let 


g = p?/h?, 


we must solve the inhomogeneous equation 


(182) 


@V/det—g@V = —g?S (183) 


subject to the boundary conditions (179) in x. 

The solution of this equation subject to these 
boundary conditions is easily obtained by the 
method of variation of parameters and may be 
written in the form 


V=S[1—sinh g(L—x)/sinh (qZ)] (184) 


and we have from our basic integral 


U(x,t)= > (Se?*) 


residues 


/pl1—sinh g(L—x)/sinh (qZ)]. (185) 


The residue at p=0 is 


S[1—(L—x)/L]=Sx/L. (186) 


If we now recall the product series expansion 
of the function sinh (z) as in Eq. (145), we see 
that the other zeros of the denominator are at 


2 


g=—n7r/L?, r—1,2,3,--- (187) 


or p= —rrh?/L?. 


Summing the residues, the solution may be 
written in the form 


U(x, t)=S[x/L+ > (2/nzx) sin (nxx/L) 
n=1 


Xexp (—n*r*h*t/L?) ]. (188) 


This is a rapidly converging series. We see that 
as [+o the temperature tends to the steady 
distribution, U=Sx/L. 


D. GENERAL CONSIDERATIONS 


From the above examples illustrating the use 
of the Laplacian transformation in the solution 
of some of the partial differential equations of 
mathematical physics, we see that the method 
is applicable where we deal with linear equations 
with constant coefficients having one-point 
boundary conditions in one independent variable, 
which in physical problems, is usually the time. 


305 





The application of the Laplacian transformation 
reduces the partial differential equation in two 
independent variables to an ordinary differential 
equation in one independent variable which may 
then be solved and then the calculation of the 
the 
problem. In the case of more than two inde- 


inverse transform gives the solution of 
pendent variables, the method of procedure is 
similar except that in this case the transform 
satisfies a partial differential equation rather 
than an ordinary one and hence the solution is 
somewhat more complicated. 


VI. Evaluation of Integrals 

Operational processes may be employed to 
advantage in the evaluation of certain definite 
integrals. Many integrals may be evaluated by 
means of the following theorem: 


If g(p)=h(t), (189) 


ge(p)dp, p= [aca t. 


then 


The proof of this theorem is as follows: 
By hypothesis 


g(p)/ p= e~?*h(t)dt. 
nasil 


Hence, 


g(p)dp p= | lJ e vn bat fd. 


e-"'dp=1/t 


“é 


and, therefore, we obtain 


g(p)dp p= | 


( 


h(t)dt/t. 


. 0 


As an example, let it be required to evaluate 


| sin (t)dt/t. 


0 


the integral 


(190) 


Now by (38), p/(p?+1)=sin (¢) therefore, 
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sin (t)dt/t= dp/ (p?+1) 


all 


-| tan-' p=7/2. (194 


As a second example, consider the integral 


(192 


[ [e-**—e bt ldt 
=~¢ t ) 


Now by (23), p/(p+a)+e~*, therefore, 


ie le at_e ot ldt i 
ze - =| (1 (p+a)—-1 (p—b) lt 


= log (b/a). (193 


Another method depends on the use of 
that 
sented. As an example, consider the integral 


parameter may be operationally repre. 


(194 


This is usually evaluated by a trick. Let (=1 
and consider 


e~ '=*dx. 


val 


Now by (23) p/(p+x"*)+e~*”’, hence 


ox 


| e dee [pds (p?-+x?) 


0 ( 


x 


us r 
=1/ptan“'x/, P| a, p. 


0 


(196 


But by (19), 4/p=1(zt)! and, therefore 


x 


| e~ "dx = (1/2)(x/t)!. 


0 


If we place ‘=1, we obtain 


x 


| e~dx=1/1/2. 


0 
As a further example, consider 


i 


0 


x 


cos (tx)dx/(1+x?). (199 
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re, let t=1/p. Now by 


) and, therefore, 


(39) cos (tx) 


tx )dx/(1+x*) 


(p?+x°)(1+x*) = (2/2) p/(p+1). (200) 
23) p/(p+1)+e~', hence, 


cos (tx)dx/(1+x?) =(1/2)e*, (201) 


1, 


cos (x)dx (1+x°*)=(2 2)1 é. (202) 


The above illustrative examples show the 
venel il procedure. 


VII. Solution of Integral Equations 


A. APPLICATION TO ELECTRIC 
CirRCUIT THEORY 


Consider a series circuit containing linear, 
lumped, constant parameters of resistance R, 
nductance, L and capacitance C, in series with 
irbitrary electromotive force e(t). The integro- 


litferential equation of such a circuit is 


at 


Ldi/dt+Ri+(1 c) | idt=e(t). (65a) 


let: I(p)=i(t) and e(t)=E(p), and 
he circuit to be initially without charge 
t, then (65) transforms into 


Z(p)I(p)=E(p), (203) 


Z(p)=Lp+R+(1/Cp). (204) 


our basic theorem (8), we see that if 


A(t)=1/Z(p), (205) 


admittance of the 
(199 ‘cuit, or the current response to a unit electro- 


is the indicial 
impressed at t=0, then we have 
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d et 
1(t) =- | A(t—u)e(u)du 


dt. 0 


t 


adr 
= -| A(u)e(t—u)du. 


dt. 0 


(206) 


illustrate the use of 
integral equations, let us consider the case of a 
variable circuit, that is, a circuit that has one 
or more elements varying with the time. To be 
specific, let us consider a series circuit containing 
the invariable elements L, R, C, in series with 
the variable elements, L,(t), R,(t), and C,(t) 
which are functions of the time, and connected 
to a source of e.m.f. e(t). 


Variable circuits.—To 


This circuit may be treated as an invariable 
circuit by eliminating the variable elements R,, 
L,, Cy, and inserting an electromotive force 
e,(t)=2(t)R, across the 
e,=(d/dt)[ L,(t)i(t) ] across 
ductance and 


variable resistance, 


the 


variable in- 


. 


e.=[1 c.)1{ i(t)dt 
“0 


across the variable capacitance. If we _ let 
€,=e_.t+er+e., then —e, represents the potential 
drop across these variable elements. We may, 
therefore, regard the total electromotive force 


impressed on the circuit to be: 


[e(t) —e,(t) ], (207) 


and obtain 
. d of ‘ * 
i(t)= -| [e(u)—e,(u) |A(t—u)du. (208) 
die 0 
In the case of a variable resistance element 
alone, we have 


d * 
i(t)= f ew @—w)du 


dt 


d ¢* 
_ | r,(uji(u)A(t—u)du. (209) 
dte 0 


In the absence of the element r,, the current 
would be given by the first term on the right of 
this equation. This, then, represents the current 
that would flow in an invariable circuit of indicial 
admittance A(t). We may, therefore, write 





d of 
r,(u)i(u)A(t—u)du. (210) 
dle 0 


i(t) =19(t) — 


Since the unknown occurs under the sign of 


integration, this is an integral equation. An 
effective method of solution is the following one. 


Let 


1(t) =t9(t) —i,(t) + 22(t) —23(t) 4 (211) 


Now if we define the terms in the series so that 


e(u)A(t—u)du, 


1,(t)= r,(u)ig(u)A(t—u)du, (212) 
dte 0 

d ;* 

r,(u)t,(u)A(t—u)du. 


dts 0 


tind 1) (4) = 


We find by direct substitution into the integral 
equation that this series formally satisfies the 
integral equation. An operational method of 
computing the terms in the series may be carried 
out as follows: Since io(f) is the current flowing 
in the invariable network, its computation 
presents no difficulty. We have A(t)=1/Z, r»(t) 
is a known time function and hence the transform 


bi(p) = rv(t)to(t) (213) 


may be computed. Then we have directly 


I\(p)= bi(p) Z(p)$11(t) (214) 


and hence 7;(f) may be computed. Now if we let 


bo(p)=r,(t)ii(t), (215) 
therefore 


I2( p) = be( p)/Z(p) = 12(t) (216) 


and i2(t) can be found. In general if 


brs ilp)=rr(t)2,(b), (217) 


then Dna i(P) = Onailp)/Z(p)Ftngi(t). (218) 


In general the operational manner of handling 
the set of integrals is simpler than direct com- 
putation. In case the magnitude of the variable 
element is small in comparison to the invariable 
elements, the series converges rapidly. 
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B. SOLUTION OF VOLTERRA’S INTEGRAL 
EQUATIONS OF THE SECOND KIND 


Let us consider the equation 


FW) =CG(t) + | G(u)K(t—u)du, (219 


0 


where F(t) is a known given function, C, js 4 
constant parameter, and K(f) isa known function 
called the nucleus. G(¢) is the unknown function 
to be determined. Eq. (219) is a Volterra’s 
integral equation of the second kind with the 
special nucleus K(t—u). If C=0, we have Vol. 
terra’s integral equation of the first kind. Let us 
introduce the transforms 


G(t)h=g(p), FU)=f(p), K(H)=kR(p). 


(220 


Then, in view of the fact that by one of our 
fundamental theorems 


at 
| G(u)K(t—u)du= g(p)k(p)/p (8a 


0 


the integral equation (219) may be transformed 
into 


f(p) =Cge(p)+2(p)k(p)/p 


(222 


g(p) =f(p)/LC+k(p)/p] 
and we have explicitly by our fundamental in- 
tegral 


1 pctim ert f(p)dp 


G(t) = : 
ico Cp+k(p) 


(223 


"= | 


As an example, let us consider Abel's integral 
equation: 
** G(u)du 


F(t)/=] — 
“9 (t—u)” 


for O0<n<1. (224 


This is a special case in which C=0 and 
K(t)=t-"*. Now by (19) we have ("°=T(1—”)/p" 


in the range 0<n<1. Accordingly 


f(p) 
p” "TU—n) 


(225 


g(p) = 


is the transformed equation. Now from our basi 
(19) we 1/pt=t*/T(st!), 
hence, we obtain 


transform have: 
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and 


1 1 
— , . (226) 
p-Trai-n) TU—n)r(n)t 
. a theorem on IT functions, we have 


1—n)n=-7/sin (mr), O<n<1 (227) 


1/py"P(1—n)=sin (nr)/rt'-" = (228) 


nd since f(p)=F(t), we have the explicit 
olution of Abel's integral equation in the form 


Ti-n) 


sin (nr) d f F(u)du 


=G(t). (229) 
T dt/Jy (t—u)'-” 


\s a second example, let us consider 


az 


cosh (x)= G(u) cos (x—u)du. 


“o 
ow from our list of transforms we have 
cosh (x)= p?/(p?—1), (41a) 
cos (x)= p*/ (p? +1). (39a) 
Hence the transformed equation is 
g( p) = p( p?+1)/(p?—1) (231) 


nd calculating the inverse transform we obtain 
the explicit solution 


G(x) =2 sinh (x). (232) 


VIII. Differential Equations with Variable Coef- 
ficients 

considering the direct operational solution 
| differential equations with variable 
is, the operational equivalents of the 
vhich satisfy the differential equations 
btained and the properties of these 

studied in an operational manner. 
imple, let the equation satisfied by the 
ctions be considered. This equation 

itten in the form 


2 1dz 
-+— —+[1—n?/u ]z=0. (233) 


udu 
ge in variables z=u-"y and u?=4x 


VOLUM MAY, 1939 


is made, the equation transforms into 
xd*y /dx*?+(1—n)dy/dx+y=0. (234) 


As in the case of equations with constant coef- 
ficients, we let 


wo 


X(p)= { e~?*y(x)dx (235) 


0 
or pX (p) = y(x). 

By integration by parts, we may construct the 

following list of formulas: 

yX, (236) 

xy— —dX /dp, (237) 

x*y—d?X /dp?, (238) 

dy/dx——y(0)+pX, (239) 


—d 
xdy /dx———(pX), (240) 
dp 


d? 
x*dy/dx——(pX), (241) 
dp* 
d*y 
——y'(0)—py(0)+p2X, (242) 
dx? 
xd*y d 
———y(0) -—-(p2X), (243) 


dx* dp 
d? . 
x*d*y/dx*’—+—(p*X), (244) 
dp? 
where we have used the notation that if 


foe) 


X(p) = [ e~P*y(x)dx 


v0 
then y—>X (p). 


If we apply these transformations to (234) and 
take y(0)=0, we obtain 


pd(pX) va 
—_—+(n—1/p)pX =0. (245) 
dp 


Integrating this differential equation in (pX) we 
obtain 


px = Cp-e-"/?, (246) 
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If we 
normalization of the Bessel functions, we obtain 


adjust the constant C to the ordinary 


the transform 


pre VPSxr?J,(2\/x) validforn>—1. (247) 


\Ve may now use our fundamental theorems 
to obtain some fundamental properties of the 
function J,. By Theorem I, (6) we have that if: 


g(p)=h(x) 
d 
ped) =- 


ax 


then 


h(x) 
if h(O)=0. Hence applying this to (247) we 
obtain 


d ‘ 
Ip [x 2J (24 x) ] 


dx 


=xle-di27 
(n/x)J n(x) +S n(x) = JT n_1(x). 


(248) 
(249) 


1(2+/ x) 


If we multiply by p" and make use of repeated 
applications of Theorem I, (6), we obtain 


=e~ P= Jy(2./x) 


(d/xdx)"[x"J n(x) ]=Jo(x). 


(250) 
(251) 
From the above result we have 
e~ VP Jy(24/x) (252) 

and if we now make use of a theorem that if 
g(p)=h(x) 


g(p/s)=h(sx) 


(253) 
then 


where s is a constant whose absolute value is 
greater than zero. If we apply this to (252), we 
obtain: 


e'/ P= Jy(2(—ix)!) = Ber(2./x)+iBei(2./x) (254) 

and 

e~/P= Jy(2(ix)!) = Ber(2./x) —iBei(2\/x). (255) 

Hence, separating reals and imaginaries, we have 
cos (1/p)= Ber(2v/x), (256) 


sin (1/p)+ Bei(2/x). (257) 


These two formulas may be used to compute 
many interesting properties of these functions, 
The plan followed in this example on Besgsel’s 
equation may be carried out in the case of 
several other linear differential equations with 
variable coefficients and the properties of their 
solutions thus studied. 


IX. Discontinuous Functions 


In many cases of physical interest, the electro- 
motive forces acting on an electrical system or 
the forces acting on a mechanical system are dis- 
continuous in nature. Usually these functions are 
periodic in nature and their effect may be com- 
puted by expanding them in a Fourier series and 
computing the effect of the harmonic terms 
separately. In case the problem is attacked 
operationally, however, the transforms of th 
discontinuous functions may be computed and 
then the computation of the inverse transform 
of the equation gives the desired solution. 














t 


Fic. 6. 





In order to give a simple illustration of the 
general procedure that may be followed in 4 
typical physical problem of this type let us 
consider the case of a pulse of voltage applied to 
a simple circuit consisting of a resistance and 
inductance in series. 

The equation satisfied by the current in such 
a circuit is 


di/dt+(R/L)i= E(t)/L. (258 
We will consider E(t) to be a pulse of voltage 
as shown in Fig. 6. It has a constant magnitude 
E and is formed at t= 7; and ceases at t= 73. Let 
us further assume that at {=0 there is an initia! 
current flowing in the circuit 7=79. Proceeding 
in the usual manner, we let 


(59 


I(p)=1(t). 4s 
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By Theorem V, (10), the transform of the 


volt pulse is 


Ele-Tip— ¢-T27) (260) 


transform of Eq. (258) is 


1 E 
( pio [e-Tir—e 
p +R L) 


4 


r]). (261) 


use of transforms (23) and (10), we have 
Rt/L. O0<t<7), 
ULt F/R—(E/R)e~RO-Tv/£, 
1, <t<To, 
WL (F/R)ye~RO-THIL 


+(E/R)e "2 Fo/£, 


(262) 


t>T>. 


\lore complicated cases may be handled in 
the same manner. The method is particularly 
valuable in case the function is not periodic and 

nce does not admit a Fourier series expansion. 
admit a Fourier series 
expansion, the application of this method yields 

solution in closed form for any desired value 
the time. 


the function does 


X. Conclusion 


The many fields in which the operational 
method is a valuable tool should be apparent 
from the above account. It has been the purpose 
of this discussion to present typical problems in 
a varied group of problems and to show how 
such problems may be attacked operationally. 
The reader who desires a more detailed account 
of the various fields of application here discussed 
so briefly will find the following bibliography to 
standard treatises and fundamental papers of 
value. It may be pointed out that some of the 
standard treatises do not base their discussion on 
the more general Laplacian transform point of 
view but follow the heuristic treatment of 
Heaviside. However, all results and theorems of 
these accounts may be translated without 
difficulty into Laplacian transform language, and 
proofs that are rather vague are seen to be clear- 
cut applications of the notion of contour inte- 
gration. It is hoped that this account will help 
the reader bridge the gap between the two 
methods of presentation, and give him a sense 
of the utility of these methods. 
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of Recent Research 





Recent advances in experimental and theoretical physics 
are described in these columns in nontechnical language. 
It is intended that sufficient background be given to 
introduce the reader to the subject. 





The Further Con- 
centration of N'° 


The early methods of 
separating isotopes de- 
pended on differences in 
physical properties related directly to mass. 
Several of these methods have been developed 
n recent years to produce substantial changes 

isotopic ratios but yield only small quantities 
of product. Because relatively large amounts of 

} isotopes are 
ork, as 


hysical properties of the isotopes and _ their 


needed for biological tracer 
well as for many investigations of 
ompounds, other methods of separation had to 
be developed. The electrolytic method proved to 
most successful for concentrating large 
mounts of deuterium but it is not practical for 
ncentrating isotopes of the other light elements. 
Small differences in the chemical properties of 
the isotopes of the light elements which were 
rst predicted by theoretical considerations of 
the isotopic molecules and later confirmed by 
\periment, made it possible for Professor Urey 
nd his co-workers to develop chemical methods 
ition.’ Recently H. G. Thode and H. C. 
reported the successful concentration of 
ogen (15) in relatively large quantities. Their 
thod depends on the isotopic exchange between 
mmonia gas and an ammonium salt solution 
ling to the equation 


+ N4“H,* solution 
NH; gas+ N"H,* solution, 
the equilibrium constant is about 1.023. 


s that if we shake ammonia gas with 


Greiff, J. Am. Chem. Soc. 57, 321 (1935); 
nd A. H. W. Aten, Jr., Phys. Rev. 50, 575 
Huffman, Thode and Fox, J. Chem. Phys. 
Thode, Gorham and Urey, J. Chem. Phys. 


le and H. C. Urey, J. Chem. Phys. 7, 34 
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in water there will be 2.3 
percent more N® in the salt solution than in the 
ammonia gas. This effect has been multiplied 
many times by using a counter current flow of 
ammonium salt in water and ammonia gas in 
apparatus of the distillation column type. The 
process consists of a flow of the liquid phase 
down through a fractionating column, liberation 
of the ammonia by the addition of sodium 
hydroxide and subsequent boiling and the return 
of the ammonia up through the column. Nitrogen 
(15) is transported down the column since it 
concentrates in the liquid phase. 


an ammonium salt 


A six-and-a-half-fold change in the concen- 
tration of N® was obtained by this method, using 





a single column. In order to produce higher con- 
centrations of N“ Thode and Urey designed a 
cascade system of columns such that a smaller 
second column operates on the product from the 
first and a third column operates on product 
from the second. This system will come to a 
steady state producing high concentrations of 
N® in much less time than a single column of 
equal length but the same size throughout. 
The three columns of the cascade reached a 
maximum concentration of 72.8 percent N® in 
four weeks of continuous operation. This repre- 
the the 
nitrogen isotopes. The columns transported one 


sents a 702-fold change in ratio of 
gram of N® per 24 hours, an amount 1000 times 
that transported by the Hertz cascade diffusion 
method. In one run Thode and Urey produced 
8.8 grams of 72.6 percent N”, an amount suf- 
ficient for many investigations of physical 
properties of N® and its compounds. In addition 
large material were 


quantities of produced 


having ten times the natural concentration of 
N® or better. This material is very valuable in 


the biological tracer work now in progress. 


Photographic Theory The grains of a photo- 
flat 


microscopic crystals of silver halide, are dis- 


graphic emulsion, 
tributed haphazardly -throughout the thickness 
of the emulsion layer. Owing to the gradual 
absorption of radiation in traversing the emul- 
sion, the grains placed at a greater depth receive 
a weaker than 


exposure those spread at or 


nearer to the upper surface of the layer. This 


changes perceptibly the photographic density 


blackening) produced by a given total exposure. 
The contribution to density due to a single grain 
depends essentially on its size (area), apart from 
its intrinsic sensitivity and since in practice every 
emulsion contains a wide range of grain sizes, the 
elementary contributions have to be summed or 
integrated over all depths as well as over all 











4 
150 1LOG(EE) 


sizes. Silber- 


stein,' based on his quantum theory of photo- 


A recent investigation by Dr. L. 


graphic exposure, shows that the integral can be 
evaluated in finite form for a_ size-frequenc) 
distribution of the exponential type, leading t 
a compact formula for the photographic density 
as function of the exposure, with the transparency 
(T) of the coated plate as parameter. This 
formula represents closely enough a number of 
observed sensitometric curves (density plotted 
against the logarithm of exposure). The effect 
of the gradual absorption may be seen from the 
diagram (Fig. 1). The curve on the right cor- 
0.1 and that o 
the left, caeteris paribus, to the ideal case of T=! 
Besides 


general depression of density there is a per 


responds to a transparency 7 


or perfect transparency. the obvious 
ceptible change of shape, and especially slope 
of the sensitometric curve. A comparison of thi 
theoretical formula 
that 


quanta, instead of one, is the minimum require 


particularly under the a> 


} 


sumption the absorption of two light: 


for the developability of a grain—with observet 


sensitometric curves leads to consistent deter 
minations of intrinsic sensitivity. 


1 L. Silberstein, J. Opt. Soc. Am. 29, 67 (1939). 


JOURNAL OF APPLIED PHYSICS 





lew Books 





Tensor 


Analysis of Networks. GABRIEL Kron. Pp. 635 
524 cm. John Wiley & Sons, New York, 1939. 


st} 


ok deals primarily with the subject of matrix 
ilgebra in a language which is simple and easily 
xl by the reader who does not possess the ad- 
of a special mathematical training. Since a 
tical treatment of the subject is thus precluded, 
ntation must content itself with a parade of 
rules, operations and concepts. The reader 
1 good exposure to all the terminology and 
tical jargon peculiar to this subject. 
k theory is used as a vehicle for illustrating the 
m of this branch of mathematics to engineering 
lhe systematic nature of the mathematical 
s claimed to be particularly effective in dealing 


such problems as are encountered in multi-winding 


VoLUM 


ers, armature windings, and polyphase networks 
litions of unbalanced operation. 
er who is unfamiliar with the state of develop- 
he art of network theory prior to the advent of 
ilysis in this field, gains the impression that this 
procedure is new and a unique characteristic of 
viewpoint. In reality the novel features intro- 
the tensor algebra are confined for the most part 
iliarities of the notation and terminology which 
[his erroneous impression is caused in part by 
ce of specific references in the text to the literature 
k theory. A bibliography is provided at the end 
e, but even this lacks a number of pertinent 
ns to the earlier literature of matrix algebra as 
network theory, for example, the work of 
ind Strecker on the analysis of vacuum tube 
reuits (1930), and that of Baerwald on 2n 
urs (1931 
h the chief aim of the volume is to demonstrate 
eness of the tensor method of approach in 
ry, the latter is almost obscured by the sheer 
he tensor mechanism itself. For the reader who 
a treatment of network theory the book is 
ble. In fact, to appreciate the value of the 
<1 it is necessary that the reader have a 
ugh understanding of fundamental network 
elore consulting this book. 
cussion of network analysis in terms of the 
esh-currents” or “junction-pair voltages,” 
fails to mention the important artifice of 
voltage sources into equivalent current 
ce versa. Instead, when dealing with net- 
h both voltages and currents are impressed, 
proach based upon what he calls the “‘or- 
ork.” His example of such a case on pp. 
inspection be immediately reduced to a 
esh or two-node problem which is solved in 
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a few lines as compared to the two pages of tensor manipu- 
lations which he gives. 

In the final chapter which is entitled ““‘The Synthesis of 
Networks,” the author works out the conditions which 
the coefficients in the equilibrium equations of a network 
must satisfy in order that certain currents or voltages 
may become independent of the applied voltages and the 
load impedances. He then gives in tensor language a 
discussion of the method of obtaining equivalent networks 
by means of linear transformations and points toward 
the greater generality of his formulation by remarking 
among other things that this transformation matrix 
(“synthesis tensor’) is not restricted (as are those of 
Cauer and Howitt) by the condition that it be real. How- 
ever, unless the transformation matrix is real the quadratic 
forms for the energy functions do not remain real and 
positive, and the transformed network is not physical and 
hence is of no practical interest! The whole question of 
physical realizability, which is the crux of the synthesis 
problem is dismissed with hardly a gesture. One glance at 
the “synthesis tensor,” Eq. (23.67) on p. 618, for the 
simple problem of Fig. 23.9 on p. 616 reveals the inade- 
quacy of the author’s grasp of what network synthesis 
implies. In the light of what the problem of network 
synthesis is generally understood to mean, one must 
conclude that there is no material of practical value in 
this chapter. 

The one significant contribution which the author makes 
is the introduction of the “‘connection tensor” into the 
mechanism of setting up the equilibrium equations for a 
network. In this way he is able to write a single tensor 
equation which is applicable no matter what the number 
or interconnection of the network branches may be. The 
connection tensor characterizes the geometrical form of 
the network or device. This unifies the theory of apparentiy 
isolated problems. On this ground the author claims, for 
example, that the salient-pole alternator and the a.c, 
commutator motor have the same theory because they 
have the same form of equilibrium equations when written 
in tensor notation, whereas the usual textbooks treat 
these devices by means of individual theories. While the 
first part of this argument may be granted, it is never- 
theless true that the theory of a device lies not in the 
general underlying equilibrium equations but in the 
detailed interpretation of them. 

In spite of some of the above remarks, I believe that the 
book has real merit if read and interpreted in the proper 
spirit. Although the tensor viewpoint in network theory 
has as yet not been instrumental in the further funda- 
mental development of the network art, it gives promise 
of becoming useful in the future and therefore certainly 
deserves recognition as an additional tool. As a first 
presentation of this subject to an engineering audience 
the book merits recognition (particularly since our col- 
leagues the mathematicians seem so far to be either 
unwilling or incapable of producing something readable 
in this field) and we should feel grateful to the author for 
all the labor which he must have put into the preparation 


of this work. E. A. 


GUILLEMIN 
Massachusetts Institute of Technology 
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The Alloys of Iron and Nickel. Vol. I. Special-Purpose 
Alloys. J. S. Marsu. Pp. 593+-xii, Figs. 290, 
McGraw-Hill, New York, 1938. Price $6.00. 


1523 cm. 


rhis is the tenth of a series of monographs on alloys of 
iron sponsored by the Engineering Foundation. The aim 
of the book, as stated by the Committee on page vi is to 
“provide a reliable foundation for further research and 
supply to the practical metallurgist, steel worker, foundry- 
man and engineer the essential information now scattered 
through more than 2000 journals and textbooks in many 
languages,” 

In the pursuit of this aim much has been included which 
is of secondary interest to the physicist, but where his 
interest is keenest he will find abundant information to 
justify purchase. The discussion of magnetism is especially 
detailed. The reviewer has been interested in this alloy 


system for eighteen years, principally because of its 
magnetic peculiarities, and has attempted to keep up 
with new information as it appeared. Upon comparing 
his own list of references with the 619 presented in the 
present monograph the completeness of coverage in the 
new review is at once apparent. To make a quantitative 
statement on this point, there are about 550 names in the 
Name Index, and only about 80 names would have to be 
added to index the combined lists. (This omits new names 
appearing in 1938 and 1939. 


papers not mentioned here shows that many of them are 


A spot inspection of the 


secondary in nature, such as reviews or discussions, or 
papers in which properties of iron-nickel alloys are reported 
along with those of quite different materials, so that their 
inclusion is incidental. 

In one respect only do omissions seem to have been 
deliberate and, in the opinion of the reviewer, ill-advised. 
Nothing is said about the electrodeposition of iron-nickel 
mixtures, the properties of such deposits, or the relatively 
simple heat treatments which may be necessary for 
rendering them homogeneous. 

lhe author’s own particular interest in phase diagrams 
is especially apparent in his discussions of several ternary 
systems. In the physical, as distinguished from the physico- 
chemical, sections the text is not nearly so critical, con- 
flicting viewpoints being presented with less guidance for 
the reader’s choice. 

It is hard to Say too much in praise of the care and 
precision with which the maierial has been assembled and 
arranged. It is easy to find what has been done, and further 
reference to original monographs will very rarely be 
and 
than in the 


The book 1S 


many of its 


necessary. This is especially true because tabular 


graphical extracts (the latter often better 


original publications) are constantly given. 
not meant to be read straight along, but 
that 


jerkiness of linked abstracts is hardly noticeable. 


sections are written so smoothly the fundamental! 


It is interesting that this book, and the only other with 


which comparison is inevitable* are both products of 


». V. Auwers, on magnetic and electrical properties of iron and its 
s. See review, J. App. Phys. 10, 186 (1939). 


industrial rather than academic enterprise. The corporate 


engineer, in fact, is the only individual who can affor, 


nowadays, to know all about a limited subject. 
L. W. McKEEHAN 
Sloane Physics Labor ory. 
Yale University 


Gas Analysis. A. McCuLtocn. Pp. 166, Figs, 3 
143 X22 cm. H. F. and G. Witherby, Ltd., London, 1938 
Price 7s. 6d. 

The author states that his Gas Analysis was writte 
primarily to meet the need for a small book on gas analysis 
suitable for students in fuel technology. It consists of siy 
Introduction, II. The Absorbents 
used in Gas Analysis. III. The Determination of Hydroge: 


chapters as follows: I. 


Carbon Monoxide and Saturated Hydrocarbons by Co 
bustion. IV. Gas Analysis Apparatus. V. Miscellaneous 
Determinations. VI. The Calorific Value of Gases. Methods 
for the determination of the various constituents of fu 
gases are described and the manipulation of five comy 
types of gas analysis apparatus is explained. The book 
contains a twelve-page appendix of useful tables, na 
and subject indices. No attempt is made to cover t! 
entire field of gas analysis. 

Ihe selection of methods is well made, especially for 
student in Great Britain. The descriptions are clear a 
sufficiently detailed in most cases so that an accurat 
trial of the method could be carried out. The inclusior 
several typical analyses and laboratory reports adds muc! 
to the value of the book as a text. The addition of materia 
on the use of a simple gas buret, the use of the gas densi! 
balance, the determination of the humidity of gases, a 
the construction of a wet test gas meter would be desirabl 
Likewise, a brief statement of the thermal conductiy 
method of gas analysis and the use of the gas interferomet 
could be added advantageously. From the point of view 
of students in the United States, a brief description 
the Shepherd,' that of 
others,? methods of the chemists of the United States Stet 


apparatus of M. Burrell 


and analysis of gases 
that of Podbielniak* and also of the National Bureau 
total 


Corporation for the sampling 


Standards method for sulfur® would increase th 
usefulness of the book. 

In any case it should be stated that the first editior 
Gas Analysis by A. McCulloch is a worth while contribu 
to the textbook literature in the field and a practica 
reference book, especially for those who do not have acces 
to the more extensive treatises on the subject. 

C. C. MELOocHI 
University of Michigar 

1M. Shepherd, Nat. Bur. Stand. J 

2 Bur. Mines Bull. 197 (1926). 

+ Pittsburgh, Pa., (1927), p. 8 

* Podbielniak, Ind. Eng. Chem., 

6 Circular Nat. Bur. Stand. No 


Research 6, 121 (1931). 


Anal 
18, p 


Ed. 3, 177 (1931 
133. 
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Here and There 





Acoustical Society of America 


ustical Society of America will hold its tenth 

meeting at the Hotel Pennsylvania in New 
m Monday and Tuesday, May 15 and 16, 1939, 
ng interesting program has been arranged: 

\ Symposium on the Measurement and Application 
\ sorption Coefficients. 
Statement of the Problem—F. V. Hunt, Harvard 

University 

\leasurements in the Reverberation Chamber— 
P. E. Sabine, Riverbank Laboraturies 

\leasurements in the Field—G. T. Stanton, Elec- 
trical Research Products, Inc. 

Open Discussion led by V. L. Chrisler, National 
Bureau of Standards, and J. Parkinson, Johns- 
Manville Sales Corporation. 
hitectural Acoustics, Its Past and Its Posstbilities. 

P. E. Sabine, Riverbank Laboratories. 

\ Modern Concept of Acoustical Design. 
Cc. <a and J. P. Maxfield, 


Research Products, Inc. 


Potwin Electrical 
e Role of the Journal of the Acoustical Soctety. 
F. Kk. Watson, University of Illinois. 

\n Ear to the Future. 
\. O. Knudsen, University of California at Los 

Angeles. 

\ Den 

Sir William Bragg, The Royal Institution, London. 


onstration of Combination Tones. 


clusions from Further Measurements of Old and 
New Violins. 
F. A. Saunders, Harvard University. 
iking Speech. 
onstrated by its Instantaneous Analysis and 
Synthesis—Homer W. Dudley, Bell Teiephone 
iboratories. 
istration Lecture on Auditory Patterns. 


Harvey Fletcher, Bell Telephone Laboratories. 


* 


Mathematical Biophysics 


number of a new booklet on Mathematical 
has just been issued. It is published as a 
to Psychometrika. It is edited by Dr. N. 
the University of Chicago. The first issue 
pers on “Mathematical Biophysics of Cell 

‘Biological Periodicities,”” ‘‘Mechanism of 
,’ “Mechanics of Viscous Bodies and Elonga- 
soidal Cells,” “‘Mathematical Biophysics of 

| “Rashevsky’s Theory of the Gestalt.” 
on seems to fulfill a very definite need and 


t 


valuable to those interested in mathematical 
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Metropolitan Section of the APS 


The program of the meeting of the Metropolitan Section 
of the American Physical Society held at 2:30 P.M. on 
Friday, March 24, in the Pupin Physics Laboratories of 
Columbia University was as follows: 

Present Status of the Analysis of Atomic Spectra— 
A. G. Shenstone. 

Cosmic Rays: The Primary Rays and the Phenomena of 
showers—K. K. Darrow. 


Cosmic Rays: The Mesotron—Enrico Fermi. 


* 


The Electrochemical Society will hold its Columbus 
meeting at the Deshler-Wallick Hotel from April 26 to 
April 29, 1939. There will be two important and timely 
symposia on “Organic Electrochemistry” and “Refractories 
for the Electrochemist.”’ 


* 


Dr. R. N. Varney, Assistant Professor of Physics at the 
Washington University, spoke to the Physics Colloquium 
of the University of Missouri on ‘The Ionization of Gases 
by Ions and Atoms.”” Dr. Varney joined the staff of the 
University of Washington in September, 1938. 


* 


Dr. E. J. Schremp, also of the Physics Department of 
Washington University, spoke to the Physics Colloquium 
of the University of Missouri on “Recent Developments 
in Electrodynamics.” 
Department of the 


Dr. Schremp came to the Physics 
Washington University from the 
Massachusetts Institute of Technology in September, 1937. 


* 


Science reports that among the men elected as fellows 
of the Royal Society at the meeting held at Burlington 
House, London, on March 16 were: ; 

M. Born, Tait Professor of Natural Philosophy, Uni- 
versity of Edinburgh, distinguished for his work in many 
branches of mathematical physics and particularly for his 
contributions to quantum theory and its applications to 
physics and chemical! physics. 

A. J. Bradley, Assistant Director of Research in Crystal- 
lography, Cavendish Laboratory, Cambridge, distin- 
guished for his methods of applying x-ray crystallography 
to elucidate the structure of metals and particularly the 
gamma-phase and order and disorder in alloys. 

G. W. C. Kaye, superintendent, Department of Physics, 
National Physical distinguished for his 
pioneer work in x-ray measurements and for his studies on 
acoustics and physical constants; has rendered valuable 
service to the Radium Protection 
National Radium Commission. 


Laboratory, 


Committee and the 


E. J. Williams, Professor of Physics, University of 
Wales, Aberystwyth, distinguished for his researches on 
the passage of electric particles through matter, and on 
individual collision processes, which have provided evi- 
dence for the existence of the heavy electron. 
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Summer at University of Chicago 


Symposium on Liquids and Polyatomic Molecules. The 
Department of Physics and Chemistry will cooperate in 
offering a symposium on liquids and polyatomic molecules 
from June 21 to June 24. This will deal with recent de- 
velopments in regard to the structure and optical properties 
of liquids and of polyatomic molecules and will be con- 
ducted by Professors Franck, Harkins and Mulliken. 

Symposium on Cosmic Rays. A symposium on cosmic 
rays, conducted by Professor Arthur H. Compton, will be 
held from June 26 to June 28. 

Conference on the Calculus of Variations. Under the 
direction of the Department of Mathematics a conference 
on the calculus of variations will be held from June 27 to 
June 30, with a number of invited addresses by members 
of the Faculty of this and other universities on topics 
connected with the calculus of variations in which the 


speakers have been especially interested. 
* 


Priestley Lectures 


The Pennsylvania State College celebrated the two 
hundred and sixth anniversary of the birth of Joseph 
Priestley with the presentation of the thirteenth annual 
series of Priestley Lectures on March 20 to 24. The 
lectures, a series of five, were given this year by Dr. 
Kenneth S. Cole, Associate Professor of Physiology at the 
College of Physicians and Surgeons, Columbia University 

Joseph Priestley was born in England in 1733 and did 
most of his work in chemistry in that country. His chemical 
investigations were carried out quite aside from his regular 
duties as a clergyman. Although he did not begin his 
scientific work until after he had reached the age of thirty, 
his enthusiasm and talent for experimental devices overt 
came his lack of previous training and unfamiliarity with 
chemistry in general and he became one of the outstanding 
chemists of his time, He is known primarily for his most 
important discovery of oxygen in the air and for the 
preparation of the first pure oxygen. His first recognition 
came, however, in 1773 when he was awarded the Copley 
Medal for correlating the then known facts concerning 
carbon dioxide and recommending that a solution of this 
gas in water (our present “soda water’’) be carried by 
sailors on long voyages for the prevention of scurvy. 

loward his later years Priestley became quite unpopular 
in England because of his liberal views toward religion and 
the American and French Revolutions. Because of this 
ill-feeling he emigrated to America in 1794 and built a 
home at Northumberland, Pennsylvania, where he con- 
tinued his experimenting and writing until his death. His 
Northumberland home was purchased in 1919 by the 
alumni of the Department of Chemistry of The Pennsyl- 
vania State College and a museum built nearby as a 
memorial to Priestley’s contributions to early American 
chemistry. 

In 1926 the members of the faculty of the Department 
of Chemistry inaugurated a series of lectures known as 
the Priestley Lectures. Since 1931 the lectures have been 


18 


sponsored by the local chapter of Phi Lambda Upsilon 
national honorary chemical fraternity, with the coo; eration 
of the Department of Chemistry. The lectures, therefore 
now constitute a living memorial to Joseph Priestley. 

rhe lectures are given each year by an authority jn , 
field related to that in which Priestley’s work falls, | havei 
chemistry. The lecturers for the past five years have bee: 
Dr. Ross A. Gortner of the University of Minnesota 
Dr. Matthew A. Hunter of the Rensselaer Polytechnj 
Institute, Dr. Warren K. Lewis of the Massachusetts 
Institute of Technology, Dr. G. H. A. Clowes of Eli Lith 
and Company, and Dr. Robert F. Mehl of Carnegie 
Institute of Technology. 

Professor Cole, this year, discussed the inter-relation 
ships between physical chemistry and cell membranes 
His lectures dealt with the properties of the cell membranes 
which have been determined by the methods of physics 
and physical chemistry. The general topic was subdivicd 
into the following parts for the five lectures: Membran 
Problems; Membrane Permeability; Membrane Imperme- 
ability; Membrane Activity; and Membrane Structur 
Ihe lectures were given in the Auditorium of the Hon 
Economics Building at State College, Pennsylvan 
beginning at 7:30 P.M. on the evenings of March 20. 2! 
22, 23 and 24. 


* 


Sigma Pi Sigma 


Sigma Pi Sigma, the national physics honor societ) 
held its convention on April 6, 7 and 8 at Columbus, O! 
with the Ohio State University Chapter as host. Delegates 
from over thirty active chapters attended the business ar 
scientific sessions, which concluded with an inspect 
trip to the Perkins Observatory. Dr. R. C. Colwell, 
West Virginia University, the retiring president, deliver 
the presidential address on ‘Electromagnetic Waves 
Radio Signals."’ After being received into honorat 
membership in the society, Dr. Paul E. Klopsteg, of th 
Central Scientific Company, spoke at an open meeting 
“Archery, a Physicist’s Hobby.” At a panel discussiot 
Dr. M. N. States, assisted by Dr. W. P. Davey, Dr. W. I 
Forsythe and others presented various phases of problems 
dealing with the placement of graduates in physics. 

Che Convention authorized the expansion of the services 
of the Placement Board of the Society and its establis! 
ment upon a permanent basis. This board will serve as 4 
agency to bring to the attention of prospective employers 
of physicists the availability of certain members ol tt 
society and will function as a general clearing house !0! 
information regarding the placement of graduates 
physics. 

The report of the Executive Secretary to the convent 
indicated that the society has a total membership of 2435 
and an active-chapter membership of 650. 

Officers elected at the close of the convention wer 
President, Dr. F. C. Blake, Ohio State University; "1 
President, Dr. Joyce C. Stearns, University of Denve' 
Executive Secretary, Dr. Marsh W. White, The Pennsy* 
vania State College. 
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Purdue Summer Session 


x Werner Heisenberg of the University of Leip- 
Visiting Professor at Purdue University during 
er Session, from July 1 to July 22. Professor 
rg will lecture on problems of nuclear physics and 
vs. Those interested in attending should com- 
with the Department of Physics to obtain in- 


regarding available accommodations on the 
* 


International Meeting on Physics 


\n International Meeting on Physics will be held in 
h. Switzerland, from September 4 to September 16, 
the Physical Building of the Federal Institute of 

vy, 35 Gloriastrasse, Zurich 7, in connection with 

e Swiss National Exhibition. The program is as follows: 


Section I. Nuclear Physics 
September 4 to 9, 1939 
P. Scherrer, Federal Institute of Technology, 


cturers: 
hr, |. Chadwick, H. Euler, O. Hahn, W. Heisen- 
Joliot, R. A. Millikan, F. Rasetti. 


Section II. Physics of the Solid Body 
September 4 to 9, 1939 
P. Niggli, Federal Institute of Technology and 
ersity of Zurich. 
ecturers: 


in Arkel, P. Debye, R. A. Fowler, R. de L. Kronig, 
LL. Pauling, W. Pohl. 


Section III. Technical Physics 
September 11 and 12, 1939 


Fischer, Federal Institute of Technology, 


ecturers 
H. de Boer, L. Brillouin, E. Briiche, R. Houwink, 
ising, R. Nitsche, W. Schottky. 


Section IV. Television 
September 13 to 15, 1939 
At the same time: 
| International Television Meeting 
I. Fischer and F, Tank, Federal Institute of 
gy, Zurich. 
urers: 
lemy, A. D. Blumlein, K. Kiipfmiiller, R. 
Okolicsanyi, F. Schréter, V. K. Zworykin. 
Section V. High Frequency 
September 16, 1939 
_ At the same time: 
High-Frequency Meeting of the 
\ssociation Suisse des Electriciens”’ 
lank, Federal Institute of Technology, 


raine and T. L. Eckersly. 


iain lectures time will be provided in every 
rtain number of smaller papers of 30 minutes 
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each. Registration of these papers must be made before 
May 31, 1939, with the Secretary of the International 
Meeting on Physics, 41 Gloriastrasse, Zurich 7, Switzer- 
land. 


7 


Summer at the University of Pittsburgh 


The University of Pittsburgh announces a second Sum- 
mer Session on the Physics of Metals to be held from July 6 
to August 11. This summer session is a cooperative effort 
on the part of physicists and metallurgists to bring out 
recent progress in the physics of metals and to discuss 
possible applications in the field of metallurgy. The fol- 
lowing tentative schedule of courses has been arranged: 
Introductory Electronic Structure of Metals, JoHN C. 

SLATER 
Elementary Atomic Physics, M. F. MANNING 
Modern Theory of Solids, WILLIAM SHOCKLEY 
Electric and Magnetic Properties of Metals, FREDERICK 

SEITZ 
Introductory Statistical Mechanics, E. U. Connon 
Theory of Atomic Vibrations, J. C. SLATER 


A booklet will be ready shortly describing these courses 
in detail and may be obtained by writing to Dr. E. 
Hutchisson, Head of the Physics Department, University 
of Pittsburgh, Pittsburgh, Pennsylvania. 


* 


Westinghouse Research Fellowships 


Announcement of appointment of the second group of 
five Westinghouse Research Fellows has been made by the 
Westinghouse Electric and Manufacturing Company. The 
first group of five was appointed last year, and all of the 
group have been reappointed for a second year of funda- 
mental research work at the Westinghouse Research Lab- 
oratories in East Pittsburgh. 


The new appointees are: 


Dr. William A. Johnson, who will work on problems of 
diffusion of metals using artificial radioactive metals pro- 
duced in the Westinghouse ‘‘atom-smasher” as tracers. 
Dr. Johnson’s home is in Washington, D. C. He graduated 
from Lehigh University with highest honors in metallurgy 
in 1935, and has just completed his graduate study in the 
Metals Research Laboratory of Carnegie Institute of 
Technology under Professor R. F. Mehl. 

Dr. Alois Langer, who will work on chemical applications 
of artificial radioactive isotopes. Dr. Langer comes to the 
Laboratories from the School of Chemistry of the Uni- 
versity of Minnesota, where he has been engaged in re- 
search in analytical chemistry with Dr. I. M. Kolthoff, 
chief of the division of analytical chemistry. 

Dr. T. A. Read, who will study the processes governing 
internal friction in metals. Dr. Read comes from Columbia 
University, where he has been doing research in the energy 
dissipated in metals during vibration, in collaboration with 
Professor S. L. Quimby and with Professor C. Zener of the 
College of the City of New York. 
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Dr. Alden IT, Ryan, who will study 


havior of gases at ultra-high radiofrequencies, correspond 


the dielectric be 


ing to wave-lengths of about ten centimeters. He will also 
conduct experiments on new methods of producing radio 


waves of this type. Dr. Ryan comes from lowa State 


College in Ames, Iowa, where his graduate research has 


been devoted to study of dielectric constants at high 


frequencies 
Dr. G. L 


processes ol thermionic emission from oxide crystals, con- 


Tawney, who will work on the fundamental 


tinuing basic studies in this field which he has been carrying 
on at Massachusetts Institute of Technology, where he has 
worked in collaboration with Professor W. B. Nottingham. 

lhese five Fellows are the second group to be appointed 
under a plan announced by the Westinghouse Company in 


the fall of 1937, according to which ten young research 


Contributors 


men will be supported for fundamental research in field 


broadly related to the electrical industry. These men work 
on basic problems not directly connected with the imme 
diate needs of the industrial research program of the con 
pany. The first group of five, whose reappointment for , 
second year was announced at the same time, include three 
nuclear physicists, Drs. R. O. Haxby, W. E. Shoupp, and 
W. E. Stephens, who are working with Dr. W. H. Wells 
on the large atom smasher, together with Dr. J. A. Hipple, 
who has constructed a mass spectrograph for molecular 
structure studies, and Dr. Sidney Siegel, who is st udying 
the effect of the order-disorder transition on the elasti 
properties of single crystals of copper-gold alloys, 

Work of the Feliows is carried on under the genera 
supervision of Dr. E. U. Condon, Associate Director of the 
Westinghouse Research Laboratories. 


to This Issue 





R. P. Johnson received his B.A. from the University 
of Richmond in 1929, his M.A. in Physics from the Uni 
versity of Virginia in 1930 and his Ph.D. from the Massa- 
chusetts 1936. 1930 to 
1931 he by the Geophysical Research 
Corporation; from 1931 to 1933 he was an instructor in 
the 


Institute of Technology in From 


was employed 


since 


of the 


mathematics at University of Richmond and 


1936 he 


General Electric Company. 


has been with the Research Laboratory 


Captain S. F. Varian was a graduate student at Cali- 
fornia Polytechnic Institute. He was a commercial pilot 
for fifteen years and from 1929 to 1936 was a Captain with 
Pan American Airways, flying between Brownsville, Texas 
and the Panama Canal. He is now a research associate at 


Stanford University. 


Russell H. Varian received his B.A. degree from Stan- 


R. H. Varian S. F. Varian 


ford University in 1925, his M.A. in 1927. From 1929 | 


1932 he was with Farnsworth Television Corporation 


R. B. 
California Institute of Technology in 1935 and his Ph.D 


Nelson received his B.S. in physics from the 


in 1938 from the Massachusetts Institute of Technology. 
He is now employed by the RCA Manufacturing Company 
at Harrison, New Jersey. 


Miklos Hetényi obtained his engineering diploma from 
the Technical University of Budapest in 1931 and, after 
spending several years in the design of Danube bridges 
came to this Country under the Jeremiah Smith Fellowshi 
in 1934. He obtained his Ph.D. degree in Engineering 
Mechanics at the University of Michigan in 1936, where 
Horace H. Rackham Fellow for 1936-37 
Since 1937 he has been engaged at the Westinghouse 
Research Laboratories. 


he became a 


M. Hetényi R. B. Nelson 
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A High Frequency Oscillator and Amplifier 


RussELL H. VARIAN AND SIGURD F. VARIAN 
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(Received January 6, 1939) 


stream of cathode rays of constant current and 

s sent through a pair of grids between which is an 

g electric field, parallel to the stream and of such 
ngth as to change the speeds of the cathode rays by 


ypreciable but not too large fractions of their initial 


ved. After passing these grids the electrons with increased 


needs begin to overtake those with decreased speeds 


i 


} 
u 


| of them. This motion groups the electrons into 


bunches separated by relatively empty spaces. At any 


t beyond the grids, therefore, the cathode-ray current 


I. INTRODUCTION 


has long been realized that apparatus for 
transmitting and receiving electromagnetic 
illations at wave-lengths of the order of 10 cm 
shorter should lead to the accomplishment of 
irge 


~ 


number of useful objectives, because 
ive-lengths of this order would permit con- 
ration of the radiation into rather sharply 
efined beams. Such waves have been produced 


many ways since the time of Hertz, and can 


be produced readily by a spark transmitter and 


other ways; but the perfection of the three- 
trod 


le tube, and in particular the ease with 
| obtains great amplifications, has raised 


our standards so far, at least for low frequencies, 


this at 


Counter 
The fu 
kinds: 


brief a 


eate a hope for similar apparatus for 

frequencies. It is the purpose of the 
paper to describe a new type of elec- 
levice which has made progress in this 


isons for certain important features of 
e may be seen most readily after a 
ission of some of the difficulties en- 
by the triode at very high frequencies. 
mental difficulties have been of two 


} 


the resonant circuits and with the 


tube its 


VoLUM: 
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can be resolved into the original d.c. plus a nonsinusoidal 
a.c. A considerable fraction of its power can then be con- 
verted into power of high frequency oscillations by running 
the stream through a second pair of grids between which is 
an a.c. electric field such as to take energy away from the 
electrons in the bunches. These two a.c. fields are best 
obtained by making the grids form parts of the surfaces of 
resonators of the type in by 


described this Journal 


Hansen. 


As a resonant circuit the ordinary coil and 
condenser become unsatisfactory as the fre- 
quency is increased, not only because it becomes 
too small to be mechanically convenient but 
even more because the losses become too great 
and therefore the shunt impedance becomes ioo 
low. Another related and yet partially inde- 
pendent difficulty is the decline in frequency 
stability as the losses increase. These difficulties 
can be, and are, partially avoided by using such 
resonators as concentric lines; but because of 
limitations due to leads and other reasons the 
resonant circuit often remains a limiting factor. 
Improvement in resonant circuits is, therefore, 
of primary importance. 

Turning to difficulties with tubes, it is always 
found, with any given triode, that, as the operat- 
ing frequency is increased too far, the efficiency 
declines and eventually the decline becomes so 
rapid that in a rather short range of frequency 
the efficiency drops from a usable value to zero. 
It has been found that this difficulty is associ- 
ated with the transit time of the electrons. 
Whenever the time of transit from filament to 
plate becomes appreciable by comparison with 
the time of a quarter-cycle, the tube fails to act 
in the normal There 
methods of attack: to increase the plate voltage, 


way. are two obvious 
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and to decrease the physical size of the tube. The 
former of these methods cannot be carried far; 
but the latter method has made tubes now com- 
mercially available for work at fairly high fre- 
quencies. The difficulties associated with this 


Fic. 1. Cross section of a copper enclosure, useful as a 
resonator. A qualitative representation of the electric field 
corresponding to the mode of oscillation used is also given. 
lhe horizontal dotted lines represent grids which should 
be of copper or other good conductor, and which should 
be designed so as to resist the current flow as little as 
possible. 


method, however, are serious: first because of 
the mechanical troubles in producing the tube, 
and, much more important, because each de- 
crease in size decreases the power capacity. It 
would seem therefore that there is a practical 
limit to useful results along this line. Accord- 
ingly, it becomes necessary to try another plan: 
not to attempt the effect of 
transit time, but to turn it to constructive use. 


elimination of 


The two objectives in this work, therefore, are 
first, the application to this problem of a new 
type of resonator, and second, a constructive use 
of the transit time. 


Il. THe RESONATORS, AND EXCITATION BY 
CATHODE Rays 


Considering first the resonators, we apply here 
the results obtained by Hansen,' who had con- 
sidered the resonant properties of practically 
closed metal vessels, like those sometimes called 
hohlraums. In our laboratory, these have been 
called the Greek 
“rhumba,”” meaning rhythmic oscillation, and 
the 


‘‘rhumbatrons,” from word 


familiar termination “tron,’’ making the 
name mean the place where these rhythmic 
oscillations occur. Hansen obtained mathemat- 
ical results on frequencies and power consump- 


1W. W. Hansen, J. App. Phys. 9, 654 (1938), and 


much information imparted in private conversation, 
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tion for shapes of sufficient variety to make 
rough estimates possible for many other cases. 
kor such resonators there are no quantities 
corresponding. exactly to the familiar self-jp. 
ductance L, resistance R and capacitance C of 
an ordinary circuit; and there are several essen. 
tially different but almost equally good ways of 
defining quantities to correspond roughly to 
them. By using any of these ways, and letting 
Q=Lw/R as usual,* any simple non-re-entrant 
rhumbatron, such as a spherical one, has a Q of 
the order of 4/6, where \ is the wave-length in 
free space and 6 is the skin depth in the metal of 
which the rhumbatron is made. The quantity 
shunt 


corresponding to impedance, by any 


reasonable definition, is then expressible in 
absolute electromagnetic units as something of 
the order of Qc, where c is the velocity of light, 
3X10" cm/sec. Since Q is a pure number, the 
conversion to practical units affects only <, 
which becomes 30 ohms. If 4\=10 cm and the 
metal is copper, (A/6) = 5.88 X 10*; so such shapes 
give high values to Q and to the shunt im- 
pedance. For reasons to be explained below, the 
shapes actually used are re-entrant, as, for ex- 
ample, in Fig. 1; and then Q and the shunt 
impedance are somewhat reduced but their 
values are still very high. 

Hansen proposed to use such resonators, ex- 
to accelerate 
electrons by running them through the resonator 
in the general and of the 


strongest electric field. So a reversal of phase 


cited by any feasible means, 


region direction 
relations gives the opposite result, the electrons 
being decelerated and the energy they lose being 
transferred to field. Thus to excite our 
resonator we make holes or grids in two of its 


the 


opposing walls adjacent to the region of high 
electric field strength, and pass through these 
grids and this field a periodically changing 
stream of electrons, the changes being so timed 
that more electrons pass when the field is such 


* This may sound as if there are various definitions 0! ( 
possible, each leading to a different numerical value. What 
is intended, rather, is that this is an equation of condition 
which, taken along with the equation LCw* = 1 forces Us, 
on choosing some particular definition of, say, L to us 
with this definitions of C and R which will make t 
lumped constant “equivalent circuit’? have the same !t¢ 
quency and the same damping ratio as the actual rhumba- 
tron. This means that, whatever definition of L one star's 
with, Q is always defined as 2x times the ratio of the ener$) 
stored to the energy lost per cycle. 
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is to take energy from them than when the field problem stated at the conclusion of our intro- 
i; reversed. The reason why a sphere, for ex- duction; namely, that of finding a constructive 
imple, is not a suitable shape of rhumbatron is _ use for the transit time. The idea of constructive 
that its wave-length is only about 1.14 times use of transit time, in general, is not new, since 
its diameter; so even if an electron had a speed for example it appears in the magnetron and 
equal to that of light, it would never get through Barkhausen-Kurz oscillators. In the present 
the sphere before the field reversed. For best case, however, the method of using this time is 
operation, the path of the electron in the reso- different, being more adapted to the present 
nator must not be much longer than 6/2, where — problem. 
3=v c¢ and v=electron velocity, and it seems This method involves the use of an auxiliary 
best to have the path still shorter. We therefore rhumbatron, in a way shown diagrammatically 
se as a resonator a copper enclosure of shape’ in Fig. 2. Here the electrons from the cathode 
venerally similar to that of Fig. 1. are accelerated by the electrostatic field produced 
It will be noted that no use is made of high by the battery. In the form shown here, using 
frequency lead wires, nor is there any solid grids and straight beams, rather than holes and 
dielectric. Thus we obtain the full benefit of the focused beams, the first grid serves to straighten 
desirable qualities of the resonator. Also, it the lines of force of this field. Beyond it the 
should be specially noted that there is no radio- electrons pass through a field-free space to the 
frequency field outside the resonator, and there- second grid, and then through the high fre- 
fore no radiation, until we add such accessories quency field of the rhumbatron to the third grid. 
is are needed for producing the radiation desired The latter field is not strong enough to change the 
for practical uses. speed of an electron by more than a small or 
very moderate fraction of its speed at the first 
II]. Transit TIME AND BUNCHED CATHODE Rays Or second grid. So the current carried by the 
cathode rays not only has the same value at all 
points in space from the cathode to the second 
grid, but also nearly the same to the third grid. 
Beyond the third grid, however, the changes of 
speed have important effects. These are under- 


lhe requirement for intermittent cathode rays, 
just described, cannot be readily met by grid 
trol as in a triode for the reasons stated in 
ur introduction, in connection with transit 


tumes. Consequently it leads back to the second — gtood best by first considering an electron which 


passes the center of the high frequency field just 
as that field is changing from opposing to helping 

electrons. At the third grid, this electron has 
practically the same speed as at the second; but 
another electron which passed the center of the 
field a few electrical degrees earlier has had its 
speed reduced; and a third one, passing a few 

degrees later, is going faster. If there is plenty of 
field-free space beyond the third grid, these 
differences in speed cause the electrons ahead and 
behind the one of unchanged speed to draw 
nearer toit. Another electron of unchanged speed, 

a half-cycle earlier or later, has its neighbors 

draw away from it. Consequently at a suitable 

distance from the third grid the stream contains 
igrammatic representation of a rhumbatron, 

Fig. 1, a battery, and a cathode which emitsa bunches of electrons denser than the stream at 
‘rons indicated by the parallel dotted lines. the third grid and separated by regions less 
S a concentric line and a coupling loop for ‘ 
iofrequency power which maintains a radio- dense. 

| inside the rhumbatron. This field alter- 
rates and retards the electrons passing 


























This stream is of just the sort required for 
driving the oscillations in the rhumbatron dis- 
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cussed in Section II. Therefore, if this stream is 
sent through such a rhumbatron as that, and if 


that rhumbatron is oscillating in the proper 


EMITTER 


BUNCHER 


CATCHER 


i 
TIME 


Fic. 3. Sequence of diagrams showing, for 12 times 
equally spaced throughout one cycle, the positions of chosén 
typical electrons in the beam shown in Fig. 2. 
As to phase, the first vertical row of dots shows the 
positions of the electrons at a time when the field in the 
‘‘buncher”’ is zero and is coming to be in the direction to 
increase the speeds of the electrons. Horizontal lines show 
the positions of the cathode and the buncher rhumbatrons. 
A suitable place for the catcher is designated by the 
horizontal dotted lines. 


dots) 


phase to retard the electrons in the bunches, and 
if its oscillating field is strong enough to take 
more energy from the electrons than was given 
them by the rhumbatron of Fig. 2, the apparatus 
as a whole converts part of the d.c. power of the 
cathode rays into high frequency a.c. power. 

Such an apparatus we call a ‘“‘klystron,”” from 
the Greek verb ‘‘klyzo,”” expressing the breaking 
of waves on a beach. The rhumbatron in which 
the electrons are given their first modulation of 
speed is called the ‘“buncher,”’ and the one into 
which the bunches ‘‘break’’ is called the “‘catcher”’ 

though it must be noted that the electrons may 
go through the catcher, if both its surfaces are 
grids, and that what it “‘catches”’ is not primarily 
current but power. 

To show diagrammatically the character of 
these bunches, in Fig. 3 we have reduced the 
spacing of the grids in the buncher to zero, and 
the and drawn vertical 


likewise in catcher, 


columns of dots each representing a line of 
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electrons in the places they occupy at one instant, 
and different columns representing such lines 
occurring at different instants throughout the 
cycle. To catch the most power, the catcher must 
satisfy three requirements: it must be placed 
where the electrons are in bunches of a certain 
best form, to be defined later in a paper by D. L, 
Webster; it must oscillate in the right phase; and 
the strength of the oscillations must be such as to 
reduce the electrons in the centers of the bunches 
just to rest. Under these conditions, as Webster 
will prove, with ideally perfect grids and other 
parts, it should convert 58 percent of the power 
of the cathode rays going through it into high 
frequency power. Practically, of course, we 
cannot claim to reach this ideal efficiency, but 
the efficiency is good. 

As the klystron has been described here, with 
the buncher driven by some independent source 
of power, such as an antenna receiving a signal, 
it is evidently an amplifier; but to convert such 
an amplifier into an oscillator, we need only 
take a small part of the oscillatory power in 
the catcher and use it to excite the buncher. 
The phase relations must, of course, be correct 
and if not otherwise adjustable electrically the 


proper phasing may be simply obtained by 


























Fic. 4. Diagrammatic representation of radiofrequency 


amplifier and detector. 
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varving the mean speed of the electrons, by 
varying the P.D. between the cathode and the 
frst grid. And, of course, a klystron almost 
converted from an amplifier to an oscillator 
makes a good regenerative amplifier. 

Another point of interest in connection with 
the amplifier relates to the energy of the electrons 
after they have passed through the second 
A substantial majority of the 
electrons give energy to the catcher; therefore 
they pass the last grid at lower speeds than that 


rhumbatron. 


it which they passed the first grid. Also, it is 
unavoidable that a small minority of the electrons 
vo through the catcher during the wrong half of 
the cycle and acquire more speed. Thus any 
device which sorts electrons according to their 
speed may be placed after the last grid of the 
catcher and be used as a detector. Perhaps the 
simplest such is a grid nearly at cathode potential. 
When there is no signal input almost all the 
electrons go through it. Then if there is a signal 
o the buncher, the bunecher bunches the elec- 
trons, Which in turn builds up a field in the 
itcher, and this field slows down some of the 
electrons so that they can no longer go through 
the detector grid: this causes a decrease in the 
current to any collector beyond the grid. Such a 
device is shown in Fig. 4. Alternatively, the grid 
may be biased so that with no signal the electrons 
just fail to pass, and a signal then causes an 
crease in the collector current. 


lhus we have shown how amplification, oscil- 


ition, and detection may be achieved. Modula- 
may be obtained by varying the beam 
; also by varying the cathode potential 


ind so making the phase relations come nearer to 


Current 


rfarther from the ideal ones noted above; and 
re many alternative methods growing out 


[V. PRACTICAL DETAILS OF A 
TYPICAL KLYSTRON 


section we describe certain con- 
| and design details that might not be 
‘or the most part these relate to a 
early oscillator which worked at about 
which is shown in Fig. 5, but informa- 
ied from other apparatus has been 
Also, in Fig. 6 we have shown a 


quency 
| of various parts of interest. 
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First we may note the means needed to obtain 
the columnar beam of electrons assumed above. 
Electrons are drawn from a plane oxide-coated 
emitting surface A through a plane grid of 
tungsten wires B. Were it not for space charge, 
the field between these two planes would be 
straight and the electrons would all leave in the 











Fic. 5. Scale drawing of essential parts of a typical 
klystron oscillator. The vacuum envelope, means for 
tuning, means for varying the feedback and various 
unimportant details have been omitted. 


same direction. But there is space charge, and to 
compensate for this ring C is added in an 
endeavor to keep the field due to the electrodes 
plus that of the space charge approximately 
straight. 

After passing through grid B the electrons pass 
through three, and in many other cases four, 
grids, the first of which is D. These grids must 
have high electrical conductivity to avoid spoiling 
the resonators; they must have high heat con- 
ductivity to dissipate the energy derived from the 
electrons striking them; and they should have as 
little projected area as possible to keep the 
number of electrons intercepted at a minimum. 
In an effort to satisfy these requirements we have 
sometimes made the grids of honeycomb form as 
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shown in the detail. They are drilled and filed 
from copper ;’s"’ thick and are perhaps the most 
tedious parts to make in the whole apparatus. 


Fic. 6. Photograph of typical cathode, 
grid and rhumbatron. 


The 


termined by the requirement that the 


“lattice constant’’ of the grid is de- 
radio- 
frequency field along the electron path should 
usually change from zero to full value in a 
distance short compared to the distance traveled 
by an electron in one half-cycle. The meaning of 
“short”? will be discussed more precisely in a 
forthcoming paper by Hansen. 

The distance between two grids of the same 
rhumbatron is usually less than the distance 
traveled by an electron in one half-cycle. For the 
buncher this would be 8/2; but in the catcher 
some of the electrons are slowed nearly to rest, so 
that 6/4 is a closer approximation. In practice 
the latter spacing is used for both so that the 
rhumbatrons can be made nearly identical, so as 
to simplify tuning. Smaller grid spacings can be 
used, but with some impairment of the shunt 
impedance. 

The two rhumbatrons are connected, mainly 
for mechanical reasons, by a tube E whose length 
is determined by a compromise between the 
the 
shorter the tube, the less distance the electrons 


following opposing considerations. First, 
have in which to bunch, with the result that, for 
best bunching, they must be given a larger spread 
in velocity. This, on the other hand, tends to 
reduce the efficiency because: (a) more power is 
required to drive the buncher as the buncher 
voltage is increased; and (b) since no electrons 
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should be thrown backward in the catcher, only 
the slowest ones can be stopped, and so a larger 
spread in velocities means more loss of energy due 
to incomplete stoppage of the faster electrons, 
On the other hand, a large bunching distance 
increases the difficulty of getting a large fraction 
of the beam into the catcher. In practice the 
compromise chosen is not critical, and the appa- 
ratus works well over a large range of distance, 
The exact shape of rhumbatron best suited for 
the not known but the 
shape shown works well and is convenient to 


present purpose is 
manufacture. 

Since the Q of the resonators is high (about 
1000) it is not easy to build them sufficiently 
nearly identical to resonate nearly enough to the 
tuning. One of the 
simplest methods of tuning is that shown, in 


same frequency without 


which the rhumbatron is slotted (along the lines 
of current flow) and then deformed elastically by 
a screw. The exact mechanism for doing this is 
not shown as it is felt that different experi- 
menters will prefer different mechanisms. 

The same idea applies to the vacuum envelope 
and pumping system, which are not shown, and 
also to the mechanism for rotating the concentri 
line which carries loops projecting into the 
rhumbatron. The purpose of this latter is, of 
course, to feed power back from the catcher to 
the buncher. 
taken from the catcher for 
radiation or other uses by a concentric line; or, 


Power may be 
for radiation, the power may be taken out by 
simply boring a hole in the catcher. 

The energy supply or “‘plate voltage” has one 
terminal connected to the cathode parts A and C 
and the other to the assembly of rhumbatrons, 
grid B, etc. Choice of a suitable voltage 1s 
governed by the following considerations. As th 
voltage is raised, at constant power, the needed 
current decreases and the needed shunt impedance 
increases, perhaps beyond what can be had. On 
the other hand, as the voltage is decreased, 
troubles with spreading of the beam by space 
Clearly the proper 
depends on the power output desired and other 


charge increase. voltage 
factors; but here also conditions are not Ver) 
critical. For example, it is possible to use many 
different voltages, from 300 to 4000, on a single 
klystron oscillator. 
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Finally we have shown a device that is not 
essential but very convenient, namely, an oscil- 
lation indicator. This consists of a magnetic field 
G produced by a small permanent magnet (not 
shown) through which a narrow pencil of elec- 
‘rons from the catcher pass enroute to the 
uorescent screen 47. When the rhumbatrons are 
not oscillating the electrons all have the same 
velocity and make a single spot on the screen. 
When oscillations occur some of the electrons 
ie retarded and others accelerated by the field 
1 the catcher, and so the spot spreads into a line. 
[his device is exceedingly convenient because of 
Also, by 


jattern on the screen one may obtain qualitative 


ts quick responses. observing the 


information about the bunching. 


V. SUMMARY 


The most important requisites for a klystron 
implifier are as follows: 


1) Efficient resonators, such as rhumbatrons, 


with holes or grids permitting cathode rays to go 
through them along the lines of force of their 
electric fields; 

(2) A beam of cathode rays going through both 
resonators and carrying power enough to make 
the catcher oscillate more powerfully than the 
buncher. 

For a klystron oscillator, the requirements are 
those of the amplifier and in addition: 

(3) A coupling loop or some _ practically 
equivalent line feeding power from the catcher to 
the buncher; 

(4) A P.D. between the cathode and the 
rhumbatrons such that the bunches arriving at 
the catcher are not too far out of phase with the 
field within it. 

Note added in proof.—Below are three references closely 
related to the above article which have come to our at- 
tention since the article was written. One of these has 
appeared since this article was submitted. A. Arsenjewa- 
Heil and O. Heil, Zeits. f. Physik, 95, 752-762 (1935); 
E. Briiche and A. Recknagel, Zeits. f. Physik, 108, 459-482 


(1938); W. C. Hahn and G. F. Metcalf, Proc. I. R. E. 27, 
106-116 (1939), 
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IV. The Experimental Variation of the Coefficient of Friction with the Strength of the Lubricant 
Source for a Complete Journal Bearing 
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mental study of the lubrication of a complete 
ring has been extended, emphasis being placed 
e played by the position and strength of the 
urce upon the coefficient of friction. As pre- 
tically, the friction on the bearing has been 
ise with increasing source strength and to be 
zher when the source is at the crown than when 
egion of highest load) of the bearing. It has 

| that when the source is at the base, where it 

il rotation effects in supporting the load, the 
be successfully operated at much lower speeds 
ids than when it is at the crown without 


INTRODUCTION 
rather definite rules have for a long 
been accepted by industry for the 
tageous location of the source of 
full journal bearings, they appar- 
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encountering the usual thin film phenomena. Moreover, for 
this case, there exists a critical value of the source strength 
at which the journal and bearing are concentric and the 
coefficient of friction is given by the simple Petroff theory. 
The experimental results confirm qualitatively all the 
predictions of the hydrodynamic theory previously de- 
veloped. Quantitatively, experiment and theory also agree 
for small values of the source strength and high values of 
the Sommerfeld variable, and especially when the source is 
at the base of the bearing, the discrepancies becoming 
marked in regions where the approximations of the theory 
are inherently inadequate. 


ently are, for the most part, of an empirical 
nature, and little information! is available con- 


1G. Welter and W. Brash (General Discussion of 
Lubrication and Lubricants, The Institution of Mechanical 
Engineers, October, 1937) have found that when the usual 
lubrication of a journal bearing is supplemented by high 
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cerning the effect of the strength of the source 
the 


bearings. 


upon operating characteristics of such 
In the first paper? of this series hydrodynamical 
the effect 


of a variation of the pressure of the source upon 


considerations were utilized to show 
the coefficient of friction of a finite journal 


Source 
Journal 


Bearing 





bearing, and the existence of this effect was 
demonstrated qualitatively by the experiments 
described in the second paper.* 

Since the original analysis was restricted to 
cases where the lubricant 90° in 
the the 
hydrodynamic theory was then extended in Part 


source is set 


advance of maximum film thickness, 
III* so as to apply to the thick film lubrication 
of finite journal bearings with arbitrary positions 
of the lubricant source, and for bearing lengths 
equal to or greater than one-half the bearing 
perimeter. In applying the theory, however, 
each position of the source has to be treated 
separately, and the actual calculations were made 
only for the cases where the source is at the 
crown or base of the bearing. 


pressure lubrication supplied in the region of highest load 
the frictional resistance at starting is very much reduced. 
hey apparently made no attempt, however, to study the 
variation of the coefficient of friction with the strength of 
the lubricant source under conditions of uniform operation. 
2 M. Muskat and F, Morgan, J. App. Phys. 9, 393 (1938). 
’F. Morgan and M. Muskat, J. App. Phys. 9, 539 (1938). 
‘MM. Muskat and F. Morgan, J. App. Phys. 10, 46 (1939). 
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The theory clearly shows that the Sommerfelq 


variable, S=(r/c)*uN/P—where r is the journal 


‘radius; c, the radial clearance; yu, the lubrican; 


viscosity; N, the journal speed; and P, the load 
per unit of projected bearing area—does no 
uniquely determine the coefficient of friction, by 
that the latter also depends upon the dimension. 
less lubricant source strength, which may be 


represented by 
go= 3uQ0r/reW, 


where W is the total load, and Q is the actual 
rate of lubricant flow through the bearing. |; 
also indicates that when the lubricant source js 
in the loaded portion of the bearing a critical 
value of go exists for which the journal and 
values of th 
Sommerfeld variable. At this critical value, whic! 


bearing are concentric for all 


depends upon the bearing length, the coefficient 
of friction is given by the Petroff equatio: 
((r/c)f=2n*S), and the curve of f vs. wN/P is 
consequently a straight line which passes throug! 
the origin of coordinates. Moreover, all of th 
(r/o)f vs. S 


constant were shown to approach the Petrofi 


curves obtained by keeping 


line as an asymptote when S becomes large 
regardless of the position of the source. In cas 
the source is in the loaded portion of the bearing 
these curves also pass through the origin 0! 
coordinates. 

In the first experimental paper of this series 
that 


lubricant source pressure the curve of frictio 


it was also shown for low values of the 


coefficient vs. Sommerfeld variable in the thick 
film region lies on a straight line passing through 
the origin. In this paper these results will l 
extended with regard to the roles played ) 
both the strength and position of the lubricant 
source upon the friction coefficient curves, and 
the latter effects will be compared with thos 
predicted theoretically. 


EXPERIMENTAL ARRANGEMENT 


The experimental apparatus was practical! 
the same as that already described,* except thi! 
copper tubing was so placed about the bearing 
that lubricant could be introduced either in te 


ase 


loaded or unloaded portion, that is, at the! 
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or crown of the bearing. The case’ in which the 
which determines the polar axis, is at 
vn is illustrated by Fig. 1. In order to 
etter torque characteristics at extremely 


source 


speeds, a sliding gear transmission was 
introduced between the Graham variable speed 
transmission and the journal. Vibration was 
oreatly reduced by mounting the driving mecha- 
nism on rubber on a separate base, the journal 


being driven through two universal joints. The 


hearing length was 9.75 in. and its radius 1.75 
giving a ratio of length to perimeter of 0.88. 
[he ratio r/c of journal radius to radial clearance 
is 1470. 
The flow was determined by collecting the 


il, that had passed through the bearing, in a 
burette. Since the flow, after leaving the bearing, 
is due entirely to gravity, the possibility of 
some error was introduced because of local 
iations in various parts of the collecting 
pparatus. The consistency of the observed 
results, however, indicates that this error was 
ipprec iable. 
all of the measurements described in this 
paper Gulfgem Oil A (viscosity 12.6, centipoises 
| 80°F) was used as the lubricant. 
\s in the previous experiments, the tempera- 
ire Was measured by a single copper-constantan 
ermocouple set in the bearing surface, 180 
legrees from the source. The speed of rotation 
{the journal varied from about jg to 50 r.p.m. 
In the first measurements the source pressure 
s kept fixed while the Sommerfeld variable, S, 
ed by changing the speed of rotation of 
il. It soon became evident, however, 
rves of (r/c)f vs. 
vuld not 


S for fixed values of 
readily be obtained by that 
because when S was varied the 
y of the system changed also and this 
ed by a change of lubricant flow. In 
to secure curves for fixed go, S was 
a parameter, and data were taken 
e variation of the coefficient of friction 
n Of go for various values of S. From 
curves of (r/c)f vs. S for several 
were constructed. In order to keep 

n the experiments, the source was in the 
en at the highest point because the bearing 
id. The change has been made in order to 


vith the geometry previously chosen, which, 
vith that ordinarily found in the literature. 
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S practically constant, changes of viscosity due 
to changes in room temperature were compen- 
sated by changes in the speed of the journal. 
The coefficient of friction was then further 
corrected, for slight variations in S, by assuming 
that the Petroff formula was valid for a small 
range, so that (r/c)Af could be expressed by 
27 AS. 

As an alternative the source pressure, p;, could 
have been utilized instead of go in the develop- 
ment of the theory. The use of go was retained, 
however, so as to avoid the introduction of the 


3-5} 
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Fic. 2. The variation of the friction coefficient with the 
dimensionless flux strength of the lubricant source for a 
journal bearing fed at the base (loaded part) of the bearing. 
Dashed curve—theoretical. Torque measured on the 
bearing. S=Sommerfeld variable. 


lubricant, a factor 
depending upon the geometry of the particular 
system, which appears in the equations connect- 


radius of the source of 


ing the strength of flow and the source pressure. 


EXPERIMENTAL RESULTS 
The solid lines of Fig. 2 show the effect of the 


variation of go upon the coefficient of friction for 
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S=0.472, 0.207, 0.095, and 0.0232, the torque 
being measured on the bearing. The source was 
in the loaded portion of the bearing. Fig. 3 
likewise shows the results obtained, for S=0.399, 
0.207 and 0.095, when the 


the source was in 


Fic. 3. The variation of the friction coefficient with the 
dimensionless flux strength of the lubricant source for a 
journal bearing fed at the crown (unloaded part) of the 
bearing. Dashed curve—theoretical. Torque measured on 
the bearing. 


unloaded portion of the bearing, 180° from the 
previous position. 

The plotting of (r/c)f against go, which is 
essentially the product of the viscosity and the 
flow per unit load, was suggested originally by 
the theory, which showed that the flux strength, 
rather than the quantity of lubricant itself, is 
the 


fundamental variable.* This result might 


have been anticipated because the effectiveness 
of the source in bringing about a displacement 
of the journal relative to the bearing must 


* An interesting physical interpretation of go is obtained 


by introducing the Sommerfeld variable S and the volume 


of the clearance space V. It then may be shown that: 
go=3SQ0/VN. 


Thus go may be expressed as 3S multiplied by the number 
of complete changes of lubricant per revolution of the 


journal, even though not all the lubricant will be replaced 
uniformly. 


30 


obviously be dependent upon the load. Moreover 
go is a dimensionless quantity while Q is not. 
It is, however, interesting to compare one of 
these curves with the components obtained by 
plotting (r/c)f against wQ, rather than against 0, 
In Fig. 4, the curve for S=0.207 of Fig. 2 has 
been replotted in this manner. It will be evident 
at once that, although the deviation from the 
average curve is quite large in the case of many 
of the points in Figs. 2 and 3, considerable 
justification exists for the adoption of the floy 
per unit of load as the fundamental variable. 
In view of the fact that in some of the plots 
of Figs. 2 and 3 the data for the three different 
the more or 
separate themselves into three different curves, 


values of load do seem less t 
it might appear that some further correctio: 
term should be introduced in order to compensate 
adequately for the variation in load. However, 
since no such definite trend is observable in all 
cases, this variation must be largely due t 
experimental error. Moreover this assumptior 
is easily justified because the deviation of almost 
any point from the curves shown is well withi: 
the experimental error. Thus in the case for 
S=0.207 of Fig. 2, the maximum deviation of 
from the average curve is 
about one percent. This corresponds to about 


any single point 


one gram on the scale pan and is of the order of 
the uncertainty in the torque determinatio 
of S the percent 
deviation from the mean is larger. The increas 


For lower values average 
of friction due to lubricant flow is also larger, 
but the value of the frictior 
coefficient is much smaller the greater fluctua- 


tions in the data are still within the experimental 


since absolute 


errors. 

Conditions might appear to be more favorable 
for accurate observations in case the load were 
increased, since for the same value of S th 
speed of rotation and consequently the torqu 
would be increased accordingly. This, however, 
is not the case. Although an increased torque !s 
favorable to greater accuracy of torque readings, 
go is at the same time reduced because of the 
the total effect to le 
observed becomes smaller. 

The broken lines of Figs. 2 and 3 are te 
curves predicted by the theory for the core 


additional load, and 


sponding values of S. All are for w, the ratio ® 
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orqut 


yevel 


que > 
dings, 


the wing length to its perimeter, equal 
to 0.88 

Fivures 5 and 6, which show (r/c)f as a function 

various values of go, were constructed by 

the values of (r/c)f from the curves of 


Fiys. 2 and 3 for several values of go up to and 


including ().4. 


Due to crowding and overlapping, only the 


experimentally determined curves are shown. 


fhe solid circles represent readings that have 
been taken from Figs. 2 and 3 in the above 


lhe curves of Fig. 5 have been extrapolated 


the origin of coordinates even though some 
tact between the journal and bearing was 
bserved when go and S were very small. If the 
journal and bearing were perfectly round and 
smooth, theory predicts that the curves should 
the trend indicated. The curves of Fig. 6 
ere broken off at the lower end because addi- 
il points could not be obtained without 
wssible damage to the journal and bearing. 
[he open circles and crosses in Fig. 5 corre- 
nd to friction data that were taken ‘by 
rving wN/P, keeping go practically constant 


1 
+} 


he critical value at which the load, according 
the theory, should be just supported by the 
pplied lubricant pressure. Under such condi- 
is the journal and bearing should be concen- 
ind the friction curve should be a straight 
rough the origin of coordinates, of slope 
in the previous paper® the value of r/c 
d for this slope is 1470. It was not found 
at the very low speeds, to keep the 
flow rigorously at the required value. 
r, in this region the flow was not quite 
for clockwise as for counterclockwise 
These probably are the reasons why 
tions of the points from the straight 
eater at very low values of S. 


DISCUSSION OF RESULTS 


h the curves are largely self-explana- 


he observed results are, for the most 
ood agreement with those predicted 
‘ly, there are some inferences that may 
worthy of note and some differences 
e experimental results and the theo- 
lictions that should be pointed out. 
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Thus when the source is in the loaded part of 
the bearing, a critical value of the source strength 
(go=0.091 for a bearing in which #=0.88) at 
which the load is entirely supported by the 
pressure due to the source is predicted by the 
theory. Actually, however, this value did not 
appear to be sufficiently large in some cases. 
When S was steadily reduced, under these 
conditions, contact between the journal and 
bearing apparently occurred at about } r.p.m. 
This became evident by a sharp reduction in 
electrical resistance and a sudden large increase 
of torque during a certain part of each revolution. 
After the source pressure was increased, contact 
no longer occurred. At gg=0.324, rotation at +, 
r.p.m. could apparently continue indefinitely. 
At this speed the friction torque was too small 
to be measured. 

These results indicate that it is possible by a 
proper choice of the position and strength of the 
source to operate a journal bearing under thick 
film conditions of lubrication at any value of the 


$-8.———- 


== | 
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Fic. 4. The variation of the friction coefficient with nQ 
(viscosity X lubricant flow) for a journal bearing fed at the 
base (loaded part) of the bearing. S=0.207. 


Sommerfeld variable. This is in direct contra- 
diction to the results usually found in the 
literature where a rise of the coefficient of 
friction with decreasing S is almost always shown 
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at very low values of S. Moreover, under such 
conditions the coefficient of friction must also 
remain finite for the case in which the torque is 
measured on the journal because the difference 


he 


ip 


>) 


Fic. 5. The variation of the friction coefficient on the 
bearing with S for a journal bearing fed at the base of the 
bearing. go=dimensionless flux strength. Solid points 
represent readings taken from Fig. 2. 


between the coefficient of friction for the torque 
measured on the journal and the coefficient for 
the torque measured on the bearing cannot be 
greater than ¢/r. 

While the agreement 


between the observed 


results and the theoretical predictions in Figs. 2 
and 3 is, in general, satisfactory for low values 
of go and high values of S, the deviations become 
greater at large values of gp and small values of S. 
\loreover, the agreement, in case the source is 
in the loaded part of the bearing, is better than 
when it is in the the 


eccentricity is, in general, larger. These results 


unloaded part where 
immediately suggest that the second approxi- 


mation of the analytical theory is not sufficient 


in case go exceeds about 0.25 for an eccentricity 
of 0.5 or more. It is believed, however, that this 
is only a part of the actual explanation, and 
that irregularities and non-uniformity in the 
journal and bearing lead, in a _ considerable 
degree, to experimental errors which tend to 
exaggerate the deviations at large values of g 
and small values of S. Inaccuracies in the journal 
or bearing could also, no doubt, account for the 
excessively high value of flow that was found 
necessary to float the journal at extremely loy 
speeds. 


("e068 616 084 032 


$s 


0-40 


Fic. 6. The variation of the friction coefficient on the 
bearing with S for a journal bearing fed at the crown 0! 
the bearing. go=dimensionless flux strength. Solid points 
represent readings taken from Fig. 3. 


A comparison of the experimental and theo- 
retical curves of Figs. 2 and 5 for the source 0 
lubricant in the loaded portion with those 0! 
Figs. 3 and 6 for the source in the unloaded 
part of the bearing shows that the coefficient o 
friction for thick film lubrication is, in general, 
a little higher for the same values of S and 9 
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when the source is in the unloaded portion. At 
very low values of S, however, the position of 
the source has a much greater effect upon the 
coeficient of friction. In fact, tendencies toward 
seizure were always observed for S<0.05 when 
the source was in the unloaded portion of the 
bearing, although operation could be continued 
indefinitely for S=0.01 with the source at the 
hase. This, of course, was due to the fact that 
vhen the source is at the base of the bearing it 
helps to support the load and the journal 
eccentricity is consequently less than that ob- 
taining When the source is at the crown. 
\lthough the curves of Figs. 5 and 6 show the 
riation of the coefficient of friction with S and 
ith gy in a satisfactory manner, it is instructive 
o study the behavior of the torque as a function 
of the load. It is evident immediately that a 
plot of the reciprocal of S as abscissa and of the 
orque as ordinate will essentially be one of 
torque vs. load. Such curves may easily be 
plotted from observed data, but a simple observa- 
tion makes the plot possible from the customary 
curve of coefficient of friction vs. S. Thus, by 


j 
lot 


finition, we have 
y ve 
f= 


c crW 
where 7 is the torque and W the load. Also, 


i €&Ff c* 
W. 


S re uN 2r8uNl 


‘, for any point on the (r/c)f vs. S curve 


-y 


ordinate (r/c)f c 


S 7 2r?>LuN 


abscissa 


of the above ratio (ordinate /abscissa) 
S are thus equivalent, in che case of a 
bearing and for constant values of 
ind speed, to plots of torque vs. load. 
. 7 several curves that were obtained in 
ier are shown. Curves 1 and 5, respec- 
rrespond to the Sommerfeld curve for 
on the journal and the Harrison curve 
rque on the bearing. The horizontal 
line at T=(4n*/c)r’LuN is for the 
se where the torque is independent of 
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the load. Curves 2 and 4 were obtained from 
the curves for go=0.5 and 0.05 of Fig. 5. 

As 1/S has been plotted as the abscissa, cases 
of fairly complete lubrication may be expected 





-. 
2r>LuN 
Fic. 7. The variation of the friction torque (7) with the 

load (W) for full journal bearings. 1, Torque on the journal 
of a journal bearing of finite length according to the 
Sommerfeld theory. II, #=0.88; go=0.5. Torque on 
bearing. Source at base. III, Petroff line. IV, ®#=0.88; 
go=0.05. Torque on bearing. Source at base. V, Torque 
on the bearing of a journal bearing of infinite length 
according to the Harrison theory. c=radial clearance; 
r=journal radius; L=bearing length; »=lubricant vis- 
cosity; N=journal speed; #w=bearing length/bearing 
circumference; go=dimensionless flux strength. Curves II 
and IV were calculated from Fig. 5. 


to lie well to the left on the curves of Fig. 7. 
In fact for S=1, the curves shown do not differ 
appreciably from each other or from the Petroff 
curve. This is because rotational] effects predomi- 
nate in this region in all of the cases illustrated. 

In the thin film region, however, i.e., large 
values of load, the character of the torque vs. load 
curve depends entirely upon the curvature of 
the friction coefficient vs. S curve. Thus if the 
latter curve is everywhere concave upward with 
a positive friction axis intercept, or has a mini- 
mum of the type usually found with practical 
lubricant testing machines or of that given by 
the Sommerfeld theory, the torque will rise 
monotonically as the load is increased. If the 
friction coefficient curve is concave downward 
in the thin film region and passes through the 
origin, as the curves of Fig. 5 do for go>0.1, 
the torque will rise and approach a constant 
value as the load increases. On the other hand, 
if the friction coefficient curve passes through 
the origin but is concave upward in the thin 
film. region, as are the curves of Figs. 5 for 
go<0.1, the torque will fall to a constant value 
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as the load is increased. And in the special case 
where the tangent at the origin is zero, the 
torque will ultimately vanish as the load is 
indefinitely increased, as in the case of the 
Harrison curve. 

The line of constant torque corresponds to the 
case in which the journal and bearing are 
concentric for all values of S. This condition 
may be realized experimentally, when the source 
is in the loaded part, by adjusting the flux 
strength to the critical value for each load. 

Comparison between the observed and calcu- 
lated source pressures for given values of lubri- 
cant flow may be made most easily when the 
flux strength has the critical value, that is when 
the pressure due to the source is just sufficient 
to carry the load. In this case the journal and 
bearing are concentric, the eccentricity is zero, 
and the pressure above atmospheric is given‘ by: 


W(1+ctnh? wo/2) 
pi= 


8rr° 


1 


x Wy—2 log 7—4 % 


— 


1 n(i+e?"”®) 


where W is the total load ; wo, the ratio of bearing 
length to journal diameter ; 7, the journal radius; 
and 7, the ratio of source radius to journal radius. 
Thus for a load of 115.4 lb., and for wo, r, and 
F equal to 2.76, 1.75 inches, and 0.0937, respec- 
tively, the calculated value of the pressure is 
26 |b./in.2 The observed value in that case was 
actually 25 |b. /in.? 

Finally, attention should be called to the fact 
that in these experiments the rotational velocities 
and loads were small so that the temperature rise 
of the bearing due to shearing of the oil film was 
scarcely detectable. Air bubbles in the oil at the 
ends of the bearing were frequently the source 
of much concern, although they usually dis- 
appeared before any measurements were made. 
Several unsuccessful attempts were made to 
obtain some correlation between the appearance 
of these bubbles and the occurrence of a pressure 
gradient tending to produce a flow of lubricant 


toward the center of the bearing. While the 
theory indicates that the probability of the 
occurrence of unfavorable pressure gradients js 
much higher when the source is at the crow 
than when at the base of the bearing, no appre¢i. 
able difference could be detected experimentally. 
This, of course, does not necessarily mean that 
the film was complete in both cases. The possi. 
bility still exists that the greater eccentricity 
necessary to obtain the same load carrying 
capacity in the case of an incomplete film might 
compensate for any decrease in friction due t 
the limited extent of the lubricating film. |; 
seems more probable, however, that with the 
small loads of these experiments the film was 
complete at all times. 


SUMMARY 


Although a detailed quantitative agreement 
between theory and experiment cannot be 
claimed, there can be no question but that the 
experimental results are at least in qualitative 
agreement with the essential features of the 
theoretical predictions. 

These are that : (1) a monotonic increase of the 
friction coefficient occurs when the pressure of 
the lubricant source is increased; (2) all of the 
curves for various values of source strength 
approach the Petroff line as an asymptote at 
high values of S; (3) in case the lubricant source 
is in the loaded portion of the bearing the 
curves pass through the origin of coordinates and 
no thin film region exists; and (4) at a certain 
value of the source strength the friction curve is 
a straight line when the source is in the loaded 
portion of the bearing. 

Finally, the position of the source of lubricant 
has been found to have considerable effect upon 
the practical operation of a journal bearing, 
particularly at low values of S. 

The authors are indebted to Dr. Paul D. 
Foote, Executive Vice President of the Gull 
Research & Development Company, for permis: 
sion to publish this paper and to Mr. D. W. Reed 
for assistance in performing the experiments. 
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Luminescence During Intermittent Electron Bombardment 


R. B. Nevson,* R. P. JoHNSON AND W. B. NotrrinGHAMt 
Research Laboratory, General Electric Company, Schenectady, New York 


(Received January 3, 1939) 


amic luminescence characteristics of five samples 
cial phosphors excited intermittently by electron 
ent have been studied. The decay of artificial 
follows closely a simple exponential law, except 
| fall during the first three milliseconds after 
ent, the prominence of this component in- 
y with the current density, An attempt is made to 
ite the observed dependence of build-up rate on 
nsity with the phenomenon of current-saturation 
continuous bombardment. A sample of calcium 


|. PROBLEM AND EXPERIMENTAL METHOD 


The course of the luminescence of a phosphor 
toward a new equilibrium value, after a change 
the intensity of the exciting agent, is an 
important datum for any theory of fluorescence. 
lhe information yielded by decay curves and 
build-up curves supplements that obtained by 
studying how the luminescence under continuous 
excitation varies with the chemical and physical 
composition of the phosphor, with its tempera- 
ture, and with the nature, energy, and rate of 
rival of the exciting particles. The behavior of 
phosphors under intermittent electron bom- 
bardment, which is the subject of this investi- 
gation, has also considerable technical interest. 
Five samples of commercial phosphors, fur- 
nished by the RCA Manufacturing Company, 
have been studied: an artificial willemite, 
ZnO: Si0g- Mn, luminescing green ; two identical 
samples of ZnS activated with Ag, luminescing 
blue; a mixed sulphide, ZnS-CdS-Ag, lumi- 
nescing yellow; and a calcium tungstate, pre- 
sumably pure, luminescing blue. The powders are 
settled onto glass disks from acetone suspensions. 
Experiments on the luminescent and electrical 
of these same samples under continu- 
ron bombardment, in the tube used for 
study, have been reported.' 


propertic 
Ous elect 


the prese 


t, RCA Manufacturing Company, Harrison, 
{ absence from Massachusetts Institute of 
me of the circuits were devised by W. B. N.; 
ors are responsible for the measurements 
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tungstate follows changes in the intensity of excitation 
within 100 microseconds, the resolving time of the 
apparatus. Two samples of ZnS-Ag and one sample of 
ZnS-CdS-Ag decay more rapidly the higher the intensity 
of excitation. The initial rate for ZnS-Ag is greater than 
would be predicted for a simple bimolecular reaction, and 
cannot plausibly be explained by non-uniformity in the 
initial distribution of reactants in a bimolecular reaction. 
The ZnS-CdS-Ag decays initially according to a simple 
bimolecular law. 


The samples are mounted off the axis of the 
tube (Fig. 1), to prevent contamination by 
material evaporated from the cathode of the 
electron gun.” The electron beam is directed by 
magnets to the phosphor under study. The 
current in the beam is controlled by the bias on 
the first grid G,; and is measured by a galvanome- 
ter in series with the second anode, a layer of gold 
covering most of the inside of the tube. Measure- 
ments with the cage structure at the end of the 
tube have shown that the beam contains 99 
percent of the current reaching the second anode. 
Two Fernico wires, completely shielded, are 
imbedded in each mounting disk and ground 
flush with the glass surface which is later coated; 
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Fic. 1, Tube and circuit for studying luminescence of 
phosphors bombarded intermittently with electrons. 


2C, H. Bachman and C. W. Carnahan, Proc. I. R. E. 
26, 529 (1938). 
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Fic. 2. Voltage waves at various points in the circuit of 


Fig. 1. 


with these probes the potential of the phosphor 
with respect to ground is measured by an 
electrometer. For each voltage V42 of the cathode 
with respect to the second anode, the second 
anode is made positive with respect to ground by 
a voltage V 4g just sufficient to bring the phosphor 
to ground potential. The difference Va42— Vig is 
then the energy of the incident electrons in volts.* 

The light from the phosphor is focused by a 
lens on an aperture in front of a vacuum photo- 
cell, and the average value of the photoelectric 
current is measured by balancing out the voltage 
drop in a platinum-film resistor, a Compton 
electrometer serving as a null indicator. That the 
current is proportional to the light intensity, 
over the range to be studied, has been established 
in a separate experiment. 

These phosphors all follow changes in the 
intensity of the electron beam far too rapidly to 
allow consecutive readings of the light during a 
single cycle of build-up and decay. It is required, 


*W. B. Nottingham, J. App. Phys. 8, 762 (1937). 
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then, to turn the beam on and off periodically 
and to measure the average value of the light 


emitted during some short part of the total 
period. The total period, the fraction of the 
period during which the beam is on, and the 
phase difference between the turning on of the 
beam and the sampling of the light, should be 
separately variable. The block diagram (Fig, 1 
shows how these requirements are met. 

Immediately in front of the photo-cell aperture 
is a rotating Duralumin disk with a hole which 
admits light from the phosphor to the photo-cell 
during 1/100 of a revolution. For suppressing the 
electron beam periodically, a square wave of 
voltage from a multivibrator circuit is added to 
the steady grid bias on the gun. During the “on” 
interval this multivibrator voltage is zero, and 
the beam current is determined by the stead, 
grid bias; in the “‘off’’ interval the grid voltage is 
—80 v, much farther than is 
necessary to stop the current completely; the 
change from one state to the other occurs in 
about one microsecond. 


carried below 


The remainder of the circuit is for synchronizing 
the multivibrator oscillations with the rotating 
disk, and for shifting the phases. Half a revolution 
from the measuring photo-cell the hole passes 
another stop diaphragm and lets light from a 











Fic. 3. Stationary monitory oscillograph pattern corre 
sponding to the cycle of Fig. 2. (a) Beam goes on; (b) beam 
goes off ; (c) photo-cell open; (d) end of sweep circuit cycle. 
By inspection, the beam persists for 0.150 cycle, and the 
measuring photo-cell is opened 0.160 cycle after the beam 
is cut off. 
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o another photo-cell. The resultant 
oulse, after amplification, synchronizes a 
conventional linear sweep circuit with the disk. 
To the sawtooth voltage output of the sweep 
circuit is added a variable negative bias, so that 
the resultant voltage is initially negative and 
rises through zero at a time in the cycle de- 
termined by the negative bias. This biased wave 
is applied to an “impulser’’ circuit which has a 
rigger action, yielding a pulse of positive voltage 
when its input voltage rises through zero. This 
pulse turns the multivibrator on and so stops the 


lamp 
curren! 


electron beam. An identical impulser, whose 
sawtooth input voltage differs from that of the 
first only in the amount of added negative bias, 
turns the multivibrator off, starting the beam. 
Rectifier tubes prevent the negative pulses which 
originate in the impulsers at the end of the sweep 
circuit cycle from reaching the multivibrator. By 
changing the two biases the beam can be turned 
on and off at any two points in the sweep circuit 
cycle; the measuring photo-cell is open, of 
midway in this cycle. Fig. 2 shows the 
voltages at several places in the circuit of Fig. 1. 

The operation is monitored continuously with 
an RCA Type TMV-122-B cathode-ray oscillo- 
graph. All time intervals are referred directly to 
the period of a 1000-cycle tuning fork oscillator, 
which is checked occasionally by beating the 17th 
submultiple of its output against the 60-cycle line 
voltage. This fork synchronizes the horizontal 


course, 


time-axis of the oscillograph at some convenient 
submultiple of 1000 cycles. On the vertically 
deflecting plates is placed a fraction of the output 
of the sweep circuit, to which is added, through 
ondenser coupling, the voltage on the grid of the 
un. The disk speed is held at 1/10 the frequency 
the horizontal time-axis. The oscillograph 
lig. 3) is then stationary, and consists 

of ten parallel traces and a vertical return trace, 
i bumps showing where the beam is 
and off. By inspection of the pattern, 
intervals can easily be read with an 

than one percent of the total period. 

he shape of the beam current wave is 

ntly monitored, on a_cathode-ray 

i (not shown) in series with the second 


s showing these circuits in detail are available 
the authors. 
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The apparatus is run at speeds from ten r.p.s., 
limited by flicker of the monitor pattern, up to 
100 r.p.s., limited by the motor rating. Tests 
have shown that the resolving time is actually 
about one percent of the period of revolution— 
100 microseconds for a disk speed of 100 r.p.s. 


II. RESULTS WITH ARTIFICIAL WILLEMITE 


The general behavior of the willemite sample is 
seen in Fig. 4, where the phosphor was excited for 
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TIME IN MILLISECONDS 


Fic. 4. Typical semi-log build-up and decay curves for 
artificial willemite, showing effect of changing current 
density. The vertical spacing is arbitrary. 


40 milliseconds and allowed to decay for 40 
milliseconds. The decay is approximately expo- 
nential, at a rate independent, after the first two 
or three milliseconds, of the current density. The 
rise toward the equilibrium luminescence is more 
rapid at the higher current densities. These same 
tendencies appear in Fig. 5, where the excitation 
period is reduced to eight milliseconds and the 
decay is followed for 72 milliseconds. Fig. 6 
shows that the shape of the decay curve is 
independent of the energy of the exciting elec- 
trons. In Fig. 7, the decay during nine milli- 
seconds, after excitation for one millisecond, is 
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followed in detail, and it is evident that the rapid 
initial decay is accentuated at the higher current 
density (curves A and B).® Curve C is added to 
show again that the course of the decay does not 
with the electron 


vary curve D, a 


reproduction of part of curve A in Fig. 6, shows 


energy ; 


that it is independent of both the excitation 
period and the total decay period. 

For this phosphor, and for all the others, a 
change in the color of the light during decay was 
sought by placing between the phosphor and the 
photo-cell a filter whose cut-off lay approxi- 
mately midway in the emission spectrum. In all 
cases, the presence of a filter changed the curves 
only by a scale factor : the color remains constant 
during decay, and the photoelectric currents are a 
true measure of the rate of emission of quanta. 

That the decay curves are nearly exponential, 
and nearly invariant under changes in energy, 
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Fic. 5. Typical semi-log build-up and decay curves for 
artificial willemite, showing effect of changing current 
density. The vertical spacing is arbitrary. 

’ In this connection, see S. T. Martin, J. App. Phys. 10, 
116 (1939), 
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density and duration of the exciting curren; 
suggests that at least during the first 72 mijj. 
seconds the decay process is monomolecular, The 
slight curvature of the semi-logarithmic plots cap 
formally be explained by non-uniformity of the 
phosphor grains. Fonda,* using ultraviolet light 
for excitation, has found that the decay charac. 
teristics of synthetic willemite depend marked) 
on particle size. Biinger and Flechsig? report 
that single crystals of KCI-T1 phosphor show ay 


accurate simple exponential decay. If these 
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Fic. 6. Illustrating that the shape of the willemi 
characteristic is independent of energy of the inciden! 
electrons. The vertical spacing is arbitrary. 


curves are composites of simple exponential 
decay curves for the individual grains, the tim 
constants for various grains must range from 
about 7 to about 14 milliseconds. 

The initial rapid decay, whose prominence 
depends on the current density, still lacks eve! 


‘G. R. 
laboratory. ; +s 

7W. Biinger and W. Flechsig, Zeits. f. Physik 67, + 
(1931). 
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terpretation. Heating of the phosphor 
mbardment and cooling during decay 
ount for it; a rough calculation shows 
temperature rise during bombardment 
eed 5°C, and the day-to-day variations 


time mperature, which had no effect on the 


om : _— ° 
fr \cteristics, were greater than this. 


ay were strictly monomolecular, and 
bability of excitation per incident 
re independent of the current density, 
p should also follow an exponential 
same time constant as the decay. It 
lowever, that 


ence 
evet 


the luminescence of 
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willemite under continuous electron bombard- 
ment increases less rapidly than the density of 
the exciting current and approaches a saturation 
value. This phenomenon of current-saturation 
should have associated with it a dependence of 
the build-up rate on current density, such as 
appears in Fig. 4. The saturation effect is 
practically independent of the initial energy E of 
the bombarding electrons—when this energy is 
increased the light output at every current 
density is multiplied by a factor proportional to 
E*. It is consistent with these observations to 
suppose that the initial energy E determines only 
the depth of penetration of the electrons, this 
depth being proportional to E? (in keeping with 
the usual theory), and that the excitation per 
unit length of path of the electron depends only 
on the concentration of excitable ‘‘centers.”’ 
Assume that this concentration in the unexcited 
sample has the limited value No, and that each 
center which has been excited is unconcerned 
until it has reverted, with emission of light, to 
its ground state. Then the concentration N, of 
excited centers should change during bombard- 
ment according to the equation 


dN,/dt=kI(No— Ni) —ANi, 


where A is the natural decay rate (found by 
observing the decay after stopping the excitation), 
I is the current density, and & is a constant. 
Integrating this equation with N,(0)=0, and 
recalling that the light output is at every instant 
proportional to the concentration of excited 
centers, we find that the luminescence should 
build up during bombardment according to the 
law 


RIN» 


m 
L(t) =uN,i(t) = —[1—e- (ee) (1) 
kI+X 


and should approach an equilibrium value L. 


given by 
1 1 d 1 
a 
L x uNo ukNo 


Fig. 8 shows that the current-saturation effect, 
for this sample of willemite, follows Eq. (2) 
fairly closely. However, the build-up curves, 
such as those of Fig. 4, follow Eq. (1) only 
qualitatively. If a value of k is deduced from the 
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Fic. 8. (A) Current-saturation effect in artificial willem- 
ite under continuous electron bombardment ; (B) the same 
data, plotted to show reasonable agreement with Eq. (2). 


rate of build-up at high current density, then at 
low current density the rise is more rapid than is 
to be expected. 


III. Resutts with CaWQO, 


The light from the calcium tungstate sample 
rises to its equilibrium value within 100 micro- 
seconds (the resolving time of the apparatus, for 
the diaphragm sizes used) after the electron beam 
is turned on, and falls by a factor of at least 2000 
(too low for detection with our photo-cell) within 
100 microseconds after extinction of the beam. 


IV. RESULTS WITH SULPHIDES 


The two ZnS: Ag samples behave alike in every 
respect. The rise of luminescence after the beam 
is turned on is too rapid to be followed with any 
precision. The shape of the decay curve depends 
strikingly on the current density (Fig. 9), the 
initial decay being much more rapid at higher 
current density. A change in energy of the 
incident electrons, on the other hand, makes 
practically no change in the shape of the curve 
(Fig. 10). Every ordinate is multiplied by a scale 
factor—the same factor by which the light output 
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Fic. 9. Typical semi-log decay curves for Zn5:: 
showing dependence of shape on current density. '™ 
vertical spacing is arbitrary. 


under continuous excitation is multiplied whet 
the electron energy is changed. 
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Fic. 11. Data of Fig. 9(B) plotted to show initial agreement 
with an empirical hyperbolic decay law. 


as E**, in the case of these samples.' It appears, 
then, that the excitation per unit path length of 
the primary electron within the phosphor is 
constant, while the total rate of loss of energy 
varies with residual energy in such a way as to 
give a penetration depth proportional to E*°. 

The shape of the decay curves is of some 
interest. The simplest hypothesis for the sulphide 
phosphors, consistent with the fact that they are 
photoconducting, is that the excitation is an 
ionization process, and that light is emitted when 
an electron recombines with an ionized center.* 
If the electrons and ionized centers are equal in 
concentration and randomly distributed, the 
decay curve should then follow the same simple 
law as does the rate of a bimolecular gaseous 
reaction : 

k 
((k/Lo)'+t)? 


4 


where Lo is the initial Juminescence and k is a 
constant of the material. Actually, all our decay 
curves for ZnS-Ag can be fitted, over the first 
five milliseconds where the photo-currents are 


’ This mechanism has recently been discussed by F. 
Seitz, J. Chem. Phys. 6,'454 (1938). 
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Fic, 12. Typical semi-log characteristic for mixed sulphide, 
ZnS-CdS-Ag. 


most precisely measured,’ by the empirical 
equation L=a/(b+4?). Fig. 11 shows the appli- 
cability of Unfortunately the 
curves could not be traced out far enough to tell 


whether the hyperbolic law, which of course 


this formula. 


® The scattering of the observations after five milli- 
seconds, which is somewhat greater than the experimental 
error, is perhaps caused by changes in the temperature of 
the samples. Several observers have reported that the 
decay of sulphide phosphors, optically excited, shows 
marked temperature dependence at longer times. 
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cannot continue to hold, eventually yields place 
to a simple bimolecular law; neither could the 
initial decay be measured at current densities |oy 
enough to give useful precision to the intercept on 
the time axis. 

An abnormally rapid initial rate in a bimo. 
lecular reaction can formally be explained on the 
basis of an initial non-uniform distribution of the 
reacting components. In the present case, a nop. 
uniform the path of each 
primary electron is improbable for the reasons 


excitation along 
advanced above. It would conceivably be possible 
to invent an initial non-uniform distribution of 
excitation within a column surrounding the path 
of each primary electron which would yield an 
initial hyperbolic decay, and then to account for 
the dependence of the rate on current density by 
an overlapping of these columns. We are the 
more inclined to view any such hypothesis as 
ad hoc, since a number of investigators," studying 
the decay of long-lived sulphide phosphors after 
optical excitation, have found empirical laws 
which may be approximated by the form 


L=a/(b+1)*, 


with values of a ranging from 0.8 upward 
beyond 2. 

The mixed sulphide sample, ZnS-CdS-Ag, is 
like the ZnS- Ag in that the decay curves are not 
affected in shape by a change in energy of the 
electrons, but are made initially more steep b 
increasing the current density. The shape of the 
curves, however, is markedly different from the 
hyperbolic followed initially by the 
ZnS-Ag decay (Fig. 12). During the first four 
milliseconds the data is fitted reasonably well by 
the empirical formula 


course 


L=a/(b+1)?; 


the deviation at longer times is in the direction of 
an exponent larger than two. It should perhaps be 
remarked that a formula of this type, in which the 
exponent exceeds two, cannot be explained by 
non-uniformity in the initial distribution © 
reactants in a bimolecular decay. 


10 E, Becquerel, Comptes rendus 51, 921 (1860); \ 
Guntz, Ann. de Chimie 6, 5 (1926); R. Coustal, J. Chim 
Phys. 28, 345 (1931); W. L. Lewchin and W. W. Antono 
Romanovsky, Comptes rendus (Doklady) de l’Acad. ces 
Sciences, U.S. S. R. 17, 95 (1937); P. J. Mulder, J. Frank. 
Inst. 225, 527 (1938); W. H. Byler, J. Am. Chem. Soc 60 
632 (1938). 
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